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Preface 



The present doctoral thesis is the result of my research worlu in the DIANA research group 
under supervision of Giinther Hormann. 

As the title "Generalized regularity and solution concepts for differential equations" suggests, 
the main topic of my thesis is the investigation of generalized solution concepts for differential 
equations, in particular first order hyperbolic partial differential equations with real-valued, 
non-smooth coefficients and their characteristic system of ordinary differential equations. 

In the case of smooth coefficients classical distribution theory offers a convenient framework 
for solving such partial differential equations. If the coefficients are non-smooth (or even 
discontinuous) the well-known fact that the multiplication of distributions cannot be carried 
out unrestrictedly limits the scope of distributional techniques. 

Note that the product of two distributions can formally be written as the pullback of a tensor 
product of the two factors by the diagonal map 5 : x {x,x), so another viewpoint of this 
limitation is that the pullback of a general distribution by a C°° -function, as defined in |23t 
Theorem 8.2.4], exists only if the normal bundle of the C°°-function intersected with the wave 
front set of the distribution is empty. 

The relation between problems arising from non-smooth pullbacks and multiplication of dis- 
tributions, is emphasized by the fact that the pullback of the initial condition by the charac- 
teristic backward flow is a solution candidate for a homogeneous partial differential equation 
of first order. 

Generalized functions in the sense of Colombeau extend distribution theory in a way that it be- 
comes a differential algebra with a product that preserves the classical product • : C°° x C°° — > 
C°°. In addition [18', Proposition 1.2.8] states that the Colombeau algebra of generalized func- 
tions allows the definition of a pullback by any c-bounded generalized function. 

So it seems that the Colombeau algebra of generalized functions offers a suitable framework 
for dealing with the various problems caused by multiplication and pullbacks. 

By means of regularization it is easy to carry over any partial differential equation with non- 
smooth coefficient from distribution theory to Colombeau theory. In Colombeau theory there 
have been developed existence results (f.e. [18l Theorem 1.5. 2], [39]) that yield solutions for 
ordinary and partial differential equations beyond the scope of classical approaches. Never- 
theless this comes at the price of sacrificing regularity (in general a Colombeau solution may 
even lack a distributional shadow). It is prevailing in the Colombeau setting that the ques- 
tion of mere existence of solutions is much easier to answer than to determine their regularity 
properties (i.e. if a distributional shadow exists and how regular it is). 

This was my motivation for doing a microfocal analysis of the pullback of c-bounded Colombeau 
^supported by FWF-grant Y237-N13 
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generalized functions, since the solution of the (homogeneous) partial differential equation can 
be written as a pullback of the initial condition by the characteristic backward flow. The re- 
sults of these investigations have been collected in an article [19j and are presented in Chapter 
5 and Appendix B in a slightly adapted form. The first section of Chapter 2 is also based on 
this article. 

A further main topic of my thesis is to compare the Colombeau techniques for solving ordinary 
and partial differential equations to other generalized solution concepts, which has led to a 
joint article [20] with Dr. Giinther Hermann. Chapter 4 contains a slightly modified version 
of this article. 

Chapter 3 contains generalized solution concepts for ordinary differential equations, which 
are of relevance when studying the characteristic systems of ordinary differential equations. 

Finally I want to point out the important role the generalized graph, as introduced in Chapter 
2, plays throughout my thesis: In Chapter 3 it enables us to give a precise characterization 
of Colombeau generalized solutions of ordinary differential equations and its relation to other 
concepts as the Filippov generalized solutions. In Chapter 5 it serves as a technical tool, 
playing a crucial role in many of the proofs. Nevertheless I believe there are still further 
applications for the concept of the generalized graph that are yet to be discovered. 

I would like to thank Dr. Giinther Hermann for his excellent support. 
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1 Introduction 



1.1 Basic notation and overview 

Notation for sets: By 

Bp{x) := {yeW'\\x-y\< p} 

we denote the closed ball around x £ with radius p > 0. If is a compact subset of 
we define 

Bp{K) ■.= {yeW\3xeK:\x-y\< p}. 
We write S'^~^ for the unit sphere {x £ \ \x\ = 1}. 

Let denote an open subset of M". We use the notation K ^ il., if K is a compact subset of 
ri. The Lebesgue cr-algebra of J7 is denoted by and the Borel cr-algebra of Q is denoted 

Nets, subnets and clusterpoints: Let K Q. If {Xe)ee]o,i] is ^ ^^t with £ K M", then 
we denote 

(1.1.1) CP ((A,),) 

as the set of clusterpoints of the net {Xe)e- For sake of brevity we skip the index set ]0, 1] in 
our notation. 

A subnet of {Xe)e is defined by 
for some map r in the set 

(1.1.2) T := [t :]0, 1] ^0, 1] | lim r(e) = o} . 

This definition of a subnet is due to [371 p.70]. It holds that A G CP ({Xs)e) if and only if 
there exists some t £ T with A = lime_>o ^t{£)- 

Sometimes we will consider the subsequence of a net {Xs)e, i.e. a sequence (Ae^)jgN where 
Ej g]0, 1] tending to zero for j oo. 

According to [311 Definition 2.2] a net (As)^ in R is said to be slow scaled, if 
3e' g]0, 1] : Vt > 0, 3Ct > such that lA^I* < Cte^\ for ah e G]0,e'] 

holds. 
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Convex sets and functions: A set C C R"' is called convex, if for all xi,X2 G C it follows 
that AiXi + X2X2 G C, when Ai, A2 > and Ai + A2 = 1. 

Let A be a closed subset of M"', then the set 

X G I X = ^ AfcOfc : Ofc G A, Afe > 0, ^ Afc = 1, iV G N I 
k=l k=l } 

is called convex hull of ^ . It is the smallest convex subset containing A. 
A function / : C ^ M U {+00} is called convex, if 

/(Aixi + A2X2) < Ai/(xi) + A2/(x2), 

where xi, X2 G C and Ai, A2 > and Ai + A2 = 1. 

A function h : M" ^ M U {+00} with h ^ +00, which is convex, positively homogeneous and 
lower semi-continuous is called a convex supporting function. Due to [241 Theorem 2.2.8] or 
|46[ Theorem 13.2] it uniquely defines a non-empty, closed and convex set by 

(1.1.3) C7 := {x G M" I V-u; G M" : (x, w) < h{w)} 

and it holds that 

h{w) = sup(x, w), 

x&C 

so there is a one-to-one correspondence between the non-empty, convex and closed subsets of 
and their convex supporting functions. 

Some important examples: The supporting function of the space is defined by hi (0) = 
and hi{w) = +00 for if 0. If xq G M"", the supporting function of the point set {xq} is 
defined by h2{w) := {xo,w). The closed unit ball Bp{0) gives rise to the supporting function 
hsiw) := p\w\. 

If Ci and C2 are convex subsets of and hi, /12 their convex supporting functions, the set 

Ci + C2 = {z \ 3x e Ci,3y e C2 : z = X + y} 

is convex and its supporting is defined by h{w) := hi{w) + h2{w). It follows by the examples 
above, that the closed ball around xq with radius p has the supporting function w 1-^ {xo,w) + 
p\w\. 

If Ci,C2 are convex, closed subsets of M" with the convex supporting functions /ii,/i2i then 
we have Ci C C2 if and only if hi{w) < h2{w) for all w G M". For a proof see [46, Corollary 
13.1.1]. 

Let A be some (non-empty) closed subset of M". Then the function defined by h{w) := 
sup^g^(x, lo) is the supporting function of the convex hull ch{A), thus 

ch(^) = {x G M" I Vw; G : {x,w) < h{w)}. 

At this point we refer to the books |24j or [1^ for a detailed presentation of the theory of 
convex sets. See also the brief presentation in the beginning of Chapter 4.3 in [23j . 
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1.1 Basic notation and overview 



Function spaces and distributions: Let denote an open subset of M". 

The letter T will always be used for real number such that T > 0. We often write Jly to mean 
]0, T[xM" with closure = [0,T] x M". The space C°°{Q) consists of smooth functions on 
Q all whose derivatives have continuous extensions to For any s G M and 1 < p < oo we 
have the Sobolev space W'^'P(R") (such that W°'P = LP), in particular H^R"^) = T4^*'2(R"). 
Our notation for //^-norms and inner products will be ||.||^ and (., in particular, this reads 
||.||q and (., .)o for the standard notions. 

We will also make use of the variants of Sobolev and spaces of functions on an interval 
J C M with values in a Banach space E, for which we will employ a notation as in L^(J;E), 
for example. (For a compendious treatment of the basic constructions we refer to |")2. Sections 
24 and 39].) Furthermore, as usually the subscript 'loc' with such spaces will mean that upon 
multiplication by a smooth cutoff we have elements in the standard space. We occasionally 
write AC{J;E) instead of W^^^{J; E) to emphasize the property of absolute continuity. 

The subspace of distributions of order k on Q {k £ 'N, k > 0) will be denoted by V'^{Q). We 
identify T)'^{Q) with the space of complex Radon measures lu, on Q, i.e., = — + ^(77+ — 
Ty_), where iy± and r]± are positive Radon measures on fi, i.e., locally finite (regular) Borel 
measures. 

As an alternative regularity scale with real parameter s we will often refer to the Holder- 
Zygmund classes C|(M") (cf. [13 Section 8.6]). In case < s < 1 the corresponding space 
comprises the continuous bounded functions u such that there is C > with the property 
that for all X ^ y in we have 

\u{x) -u{y)\ ^ ^ 
\x — 2/1** ~ 

Notation and basic outline of the theory of Colombeau algebras: Let us recall the basic con- 
struction: A Colombeau (generalized) function on some open set il. C is defined as equiv- 
alence class {ue)e of ncts of smooth functions li^ S C°°{Q) and e G]0, 1] subjected to some 
asymptotic norm conditions (with respect to e) for their derivatives on compact sets. We 
have the following: 

(i) Moderate nets £m {^)- Me e C°°(J7)]0'i] g^^j^ ^j^^^^^ f^^, gjl K n and a G N", there 
exists p € M such that 

(1.1.4) sup ||5"ue(x)|| = 0{e~P) ase^O. 

x&K 

(h) Negligible nets 7^ (J^): {ue)e G C°°(l^)]°'^] such that, for all K (£ Q and ah g G M an 
estimate 

sup ||ue(x)|| = 0{e'') as e ^ 0. 

x&K 

holds. 

(iii) £m (^) is a differential algebra with operations defined at fixed £, M (0) is an ideal and 
Q{0,) := £m {^) /AA (0) is the special Colombeau algebra. 

(iv) If we replace the nets of smooth functions by nets of real numbers in (i)-(iii) we obtain 
the ring of generalized numbers M. 
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(v) There are embeddings, a : C°°(J7) ^ as subalgebra and 

i : G{0,) as linear space, commuting with partial derivatives. 

(vi) 0, G{^1) is a fine sheaf and Gd^) denotes the subalgebra of elements with com- 
pact support; using a cut-off in a neighborhood of the support, one can always obtain 
representing nets with supports contained in a joint compact set. 

Regular Colombeau functions: The subalgebra G°°{^) of regular Colombeau (generalized) 
functions consists of those elements in Gi^) possessing representatives such that the estimate 
(|1.1.4p holds for a certain m uniformly over all a E N". 

Rapidly decreasing Colombeau functions: A Colombeau function in ^(M") is called rapidly 
decreasing in the directions T C 5"^^, if it has a representative with the property that there 
exists a G No, such that for all p G Nq 

sup (1 + A2)P/2|n,(Aei)| = 0{e-^) as e ^ 

{A,5i)6K+xr 

holds. 

Fourier transform of a compactly supported Colombeau function: The Fourier transform of 
u £ Gc i^) is the Colombeau function in ^(M") defined by 



J'iu) : = 




where (ne)e is a representative of u with joint compact support in 0. 

^°°-wavefront set of a Colombeau function: If v G Gc {^), we define the set T,{v) c S""^^ to 
be the complement of those points having open neighborhoods F C S'"'^^ such that J-{v) is 
rapidly decreasing in the directions F. is a closed subset of Now let u £ Gi^). 

Then we define the cone of irregular directions at xq by 

•pectin), ip{xo)j^o 

Then the (generalized) wave front set of u is the set 

WF{u) := {{x,0 e X e 

We denote the projection onto the first component by 

singsupp(n) := {x G ft\{x,^) € WF{u)} 
and call this (generalized) singular support of u. 
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1.2 Set-valued Maps 



C-bounded Colombeau maps: Let $7i C and 1^2 ^ be open sets. By ^[$71,172] we 
denote the generalized maps F £ Q{0,i)"^ with the property that F is c-bounded on 0,i (into 
^2), i.e. it possesses a representative (-Fe)e with G C°° {Qi,Q2) satisfying the condition 

yK d ni,3K' (s ^2 : such that ^^^(i^) C K' for all e g]0, 1]. 

Our standard references for the foundations and some applications of Colombeau's nonlinear 
theory of generalized functions are [6l 171 1431 118j . We will employ the so-called special variant 
of Colombeau algebras, denoted by in [18], although here we shall simply use the letter Q 
instead. 

For microlocal analysis in context of Colombeau theory we refer to [111 BOl [26l \30\ [T6] . 

1.2 Set-valued Maps 

Here we introduce the notion of a set-valued map (also referred to as correspondence) and 
a few of their properties. This concept plays an important role in the field of mathematical 
economics and game theory. A main tool here is the generalization of Brower's theorem by 
Kakutani (see [361 ) for set-valued maps. Nevertheless we are not intending to give a complete 
overview of the topic and restrict ourselves to a few basic results, which will be needed in the 
existence proof for ordinary differential inclusions. 

The interested reader is referred to [SSj and [3]. 

Let n, ill be arbitrary non-empty subsets of and let 0,2 be a non-empty subset of . 

Throughout this section we define the projection maps of the product space of fli x ^l2 onto 
ill resp. 0,2 by 

TTi : Oi X fl2 ^ Oi, {x,y) X resp. 112 ■ Oi x 0,2 ^ 0,2, (x,y) 1— > y. 
Definition 1.2.1. We introduce the following notions: 

(i) Vq{0) denotes the set of non-empty subsets ofO, 

(ii) J^o{^) denotes the set of non-empty, closed subsets ofO, and 

(Hi) 1Cq{0) denotes the set of non-empty, closed and convex subsets ofO. 

Note that for M,N £ Vo{0) we have 

X-M :={y£0\3xeM :y = X-x}e Vo{0) 

for A G M, and 

M + N := {y e \ 3xi £ M,X2 € N : y = xi + X2} G Vo{0). 

It follows that Vo{0) is a linear space and J-oi^) is a linear subspace of Vo{0). Furthermore 
ICo{0) is a linear subspace of Tq{0). 
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For any M, N £ J-^oi^) we put 

p+{M,N) :=inf{pG [0,oo] : CB^(M)}, p_(M,iV) :=p^{N,M), 

then we define the (extended) Hausdorff metric by 

d{M,N) := max {p+{M,N),p^{M,N)). 

for all M,N £ ^o(^)- Note that in case ofVo{^) the map d : Vo{0,) x Vo{0,) is no longer a 
metric, since (i(5i(0), Si(0)°) = 0, but 5i(0)\Bi(0)° = 5"-^ / 0. 

We put \M\ := d{M, {0}) = sup^g^^ 

Proposition 1.2.2. The linear space {J-o(M^),d) is a complete metric space. Its linear sub- 
space (/Co(M"),(i |/Co(]R")x/Co(i:")) o.go-'i''"' complete. 

Proof. We refer to [38, Proposition 4.3.10] and [38, Proposition 4.3.11]. □ 

Remark 1.2.3. In the case of ICq{^}) we can write the Hausdorff metric in terms of convex 
supporting functions. Recall from the introduction that each M G ICo{Q) has a uniquely 
determined convex supporting function h : 0, {+00} , defined by 

h{w) = sup (a, w) 

and the set M satisfies 

M = {a G I tt; G M" : (a, < h{w)}. 

Then we can write p±{M, N) and the Hausdorff metric d{M, N) in terms of the convex sup- 
porting function by p+{M, N) = sup^^gj^n max {hN{w /\w\) — hM{w /\w\)., 0} and p-{M, N) = 
sup^gffin max {hM{w/\w\) — hN{w/\w\)),0}, thus 

d{M,N) = max {p+{M,N),p^{M,N)) = sup \hM{w/\w\) - hN{w/\w\)\. 
and we have \M\ := d{M,{0}) = sup^j^-^n \hM{w/\w\)\. 

Proposition 1.2.4. Let (Mj^gjo^ij be a net in J^(){Q) such that for some K i^Q, it holds that 
C K for all l g]0, 1]. If M £ J^Q{i}), then the following statements hold: 

• lim.^o p+(M„ M) ^ =^ M D limsup,_,o M,. 

• lim.^o p- (M„ M) ^ M Q liminf.^o M,. 

• lim.^o d{M„M) ^0 <^ M = lim inf ,^0 M, = \im sup,_o M, . 

Recah that limsup.^o^^ — fl.'elo,!] Uelo,.'] I'esp- liminf.^oM, := U.'e]o,i] fl^elo,.'] 
denote the set-theoretic limes superior resp. limes inferior. 

Proof Follows by \W, Proposition 4.2.2] and [Ml Theorem 3.3.11]. □ 
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Definition 1.2.5. A map defined by 



is called a set- valued map. If F maps into J'-o(^2) resp. A^o(f^2) uj^ call F closed- valued resp. 
closed- and convex- valued. 

If X CQi we use the notation Fx := Uiex 

Putting 

Grapli(F) := {{x,y) G Hi x ^2 | y G F^}, 

we say F is closed, i/Grapli(F) is closed in the relative topology of x Note that we 
use the upper-case spelling Grapli(F) to distinguish it from the classical graph grapli(/) of a 
continuous function f . 

F is called proper, if Graph(F) n fii x g x fi2 for all K <^^2- 

F is called upper semi-continuous at x, if for any W Q2 being an open neighborhood of Fx 
there exists some neighborhood X C. fti of x such that Fx C W holds. We say F is upper 
semi-continuous, if it is upper semi- continuous for all x G Qi. 

F is called locally bounded, if for all x there exists a neighborhood X C Oi such that Fx is 
bounded. 

F is called bounded, if the range of F denoted by Ran(F) := Fq-^ is a bounded subset of 0,2- 

Definition 1.2.6. Let f : ili ^ 0,2 be a (single-valued) function, then we say f is a selection 
of the set-valued map F, if f{x) G Fx for all x E fti. 

Definition 1.2.7. Let F : Qi ^ ^0(^2) be a set-valued map. Then the function 

Hp -.ni xW^ ^RU{+oo}, (x,w)^hxiw), 

where hx{w) = sup^^p^{a,w) is called the supporting function of the set-valued map F. It is 
the convex supporting function (cf. Introduction) of the convex hull of F^. In particular, if F 
is convex-valued it holds that h^ is the convex supporting function of F^, which implies 

= {a G O2 I Vu; G M"" : (a, w) < Hf{x, w)}. 

// F : Oi — > ^0(^2) is a set-valued map (allowing non-convex values), then we define the 
convex hull of F by 

ch{F) : J^i ^ /Co(J^2), ch{Fx). 

It holds that ch{F)a: = {a G O2 | Vu; G R"* : {a,w) < Hf{x,w)}. 

Example 1.2.8. Let F : M 1— /Co(M) be the set-valued map defined by 

( {0} a;G]-oo,0[ 
Fx := <^ [a,f3] x = 

( {1} x€]0,+oo[ , 

then F is upper semi-continuous if and only if a < and P > 1. Obviously F is bounded. Its 
supporting function H[x, w) is defined by 

'0 X g] - oo,0[ 
H{x,+1) := < (3 X = 

1 xG]0, -|-oo[ 
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and 

( xG]-oo,0[ 

H{x,-1) := < -a x = 

( 1 X g]0,+oo[ , 

which are upper semi-continuous in x if and only if a < and /? > 1. 

Lemma 1.2.9. Let F : J7i — ^ ^0(^2) be a closed set-valued map, then it holds that 

n U ^y=^^- 

e€]0,l] y€Be{x)nni 

Proof. We put M := flselo,!] lJyeB,{x)nni - ^2, then M D F^; foUows. We are going to 
show that M = Fx- We prove by contradiction: Assume there exists some a' G MjFx^ then 
it follows that (x, a') is contained in the open set (i^i x ^2)\ Graph(F). So there exists some 
neighborhood X C fii of a; and C of a' such that (X x W^) f) Graph(F) = 0. In 
particular a' ^ Fy for y G X, so \^y^Be{x)^o.\ — ^ small enough. It immediately 

follows that 

£e]0,l]j/6Be{a:)nQi 

a contradiction to a' G M, so the statement follows. □ 

Proposition 1.2.10. Let F : Oi — > ^0(^2) be a set-valued map, then the following statements 
are equivalent: 

(i) F is closed and locally bounded, 
(a) Graph(F) r\ K x <^ x ^2 for all K 

Proof. (^) {ii) : Let K <^Q,i, then Graph(F)niir x ^2 C iTx Fr- is closed and bounded {Fx 
is bounded since F is locally bounded), [ii) =^ (i) : Let X be some bounded subset of Qi, then 
there exists some compact set K with X C K. It follows that Fx = 7r2(Graph(F) nX x U2) Q 
7r2(Graph(F) Ci K x ^2) is bounded, thus F is locally bounded. Let (x^jag) G Graph(F) be 
a net converging to some {x, a) G x ^2, then there exists some companeighborhoodood 
K of X. Since Graph(F) D K x Q2 ^1 x ^2 and G -fC for e small enough, we obtain 
(x, a) G Graph(F) Ci K x Q2, thus F is closed. □ 

Proposition 1.2.11. Let F : Qi ^ "^0(^2) be a proper set-valued map. Then for all AC. 
closed, we have that Fa is a closed subset of 4^2 ■ 

Proof. Since F is proper we have that Graph(F) n Oi x i^T is compact for all K ^ Qi, thus 
7r2 |Graph{F) is proper and continuous, implying that Fa = 7r2(Graph(F) D A x O2) is closed, 
if ^ is a closed subset of Oi. □ 

Proposition 1.2.12. Let F : J7i — 7^0(^2) be a set-valued map, then F is upper semi- 
continuous if and only i/7ri(Graph(F) n {fli x A)) is closed in ili x il2 for all closed subsets 
A ofn2. 
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Proof. If F is upper semi-continuous, it holds that for any open neighborhood W F^, 
there exists some open neighborhood X of x, such that Fy C W. This is equivalent to the 
statement: For all open sets C we have that U{W) := {x e fli \ e F^ : a e W} is 
open. Thus Qi\U{W) = {x e Qi \ 3a e F^ : a e W} = 7ri(Graph(F) n x W)) is closed 
in showing that F being upper semi-continuous is equivalent to 7ri(Graph(F) n (Oi x A)) 
being closed, whenever A C is closed. □ 

Example 1.2.13. Let / : ili ^ be a (single-valued) function, then the set-valued map 
F : Qi /Co(i^2), X {fix)} is upper semi-continuous if and only if / is continuous. Note 
that Graph(F) = graph(/). Note that the equivalence in Proposition 1 1 . 2 . IT] translates to the 
classical equivalence: / is continuous if and only if f~^{A) is closed whenever A is closed. 
Since F^ = {f{x)} we have 7ri(Graph(F) n (Qi x A)) = {x e S^i | f{x) e A} = f-\A). 
Furthermore we observe that Fa ■= 7r2(Graph(F) n vl x 172) = /(^) is closed whenever 
^ C $7]^ is closed, if and only if / is a closed map. 

Theorem 1.2.14. Let F : Qi ^ J^o(i^2) be a set-valued map, then the following statements 
hold: 

(i) If F is closed and locally bounded, then F is upper semi- continuous and locally bounded, 
(a) If F is upper semi- continuous, then F is closed. 

(Hi) If F is convex-valued, then F is upper semi- continuous and locally bounded if and only 
if its supporting function H has the property that x i-^ Hix, w) is upper semi- continuous 
for all u> G M™. 

Proof, (i): Assume that F is not upper semi-continuous at x. Then there exists an open 
neighborhood W of F^, such that any neighborhood X of x contains some y with FyDW^ ^ 0. 
We prove by contradiction: Assume that F is closed and locally bounded. 

Picking the family of neighborhoods {Bi,{x) H f7i)£g]o,i] can find G B^{x) Pi Qi with 
Fxs n 7^ 0, thus we can find some net (a;^, ae)eg]o,i] with S Fx^ n W^. By Proposition 
11.2.101 we have that Uj/eBi(a;){y} x -Py is a compact subset of O2, implying that 

0/CP((x„a,),e]o,i]) ^ I U Mx^y)' 

and since — > x we obtain CP ((ae)£G]o,i]) ^ ^x- But G implies CP ((a£)ee]o,i]) ^ 
thus Fx n 5 CP ((ae)ee]o,i]) 7^ which contradicts F^ C W . So F cannot be closed and 
locally bounded. 

[a): If F is not closed in x ^2 there exist a net {(xejae)^} C Graph(F) converging to 
some (x, a) E ili x f72 such that a Fx. Let W be an open neighborhood of Fx such that 
a W . We prove by contradiction: Assume that F is upper semi-continuous at x, then there 
exists a neighborhood X oi x such that Fy <^W for all y £ X. 

It follows that Fx^ ^ W (for small e) which implies a = lim^^o £ yielding a contradic- 
tion to the initial choice of W. 
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[iii): Assume F is upper semi-continuous and locally bounded. By (ii) F closed. Let if G 
be arbitrary, then 

lim sup H{y,w) = inf sup H{y,w) = inf sup sup(a,ti;) 
y^x -Jelo,!] ye_Bj(a;)nni -^gIo,!] j,eB«(x)nni aeF^ 

and putting Ks := V}yeBi(x)r\^i{y} ^ -^y which is a descending family of sets, we can estimate 

inf sup sup(a,u;) < inf sup {a,w). 

Se]o,i] yeBs{x)nni a^Fy se]o,i] (z,a)eKs 

Lemma ll.2.91 implies that n<5e]o i] ~ i-^J" ^ ^^'^ Lemma I A. 41 yields 

inf sup {a,w) = sup {a,w) = sup{a,w) = H{x,w) < oo, 

which shows that x H (x , w) is upper semi-continuous for all w E M"'. 

Assume x ^ H{x,w) is upper semi-continuous: If {(xe,Oe)e} C \JyfzQ-^{y} x Fy is a net 
conyerging to some (x, o) £ Qi x Q2, then we have that {as,w) < H{xs,w) for all w S M" 
and e e]0, 1]. For w E M" fixed, it holds that 

{a,w) = lim(a£,u;) < lim sup H{xs,w) < H{x,w) < 00, 

for all w E M" because x H{x,w) is upper semi-continuous, thus a E F^. It follows that 
Uxeriii^} X -^a; is a closed subset of Oi x O2. Let K Q,i then 

sup \Fx\ = sup sup \a\ < sup sup H{x,w/\w\) < sup H{x^w) = H{xq,wq) < 00, 

xGK x&Ka^Fx xeKweR'" {x,w)£KxS"^-'^ 

for some {xo,wo) £ K x S*™^^. By ^24i Proposition] the convexity of the map w ^ H{x,w) 
implies its continuity, so (x, w) H{x, w) is an upper semi-continuous function and Lemma 
lA. 71 states that an upper semi-continuous function attains its supremum on a compact set. □ 

Example 1.2.15. Let u E be a distribution, then the wave front set of u defines a 

closed set-valued map by 

S : singsupp(n) /Co(S'„_i), x ^ ^.^{u), 

where 'Sx{u) denotes the cone of irregular directions of u at x. Since locally boundedness 
is implied by the compact range of S we obtain that T, is upper semi-continuous, i.e. if 
X E singsupp(n), then for any open neighborhood W ofT,x, there exists an open neighborhood 
X C r^i of X such that 5]jYpisjngsupp(ti) C W. 

In the same way the wave front set of a Colombeau function u E ^(r^i) defines an upper 
semi-continuous set-valued map : singsupp(n) — > ICo{Sn~i)- 

Definition 1.2.16. Let F : ^li ^ "^0(^2) be a set-valued map. 
We define the inverse graph by 

Graph(F)~^ = {{y,x) E Ran(F) x \ {x,y) E Graph(F)}. 



10 



1.2 Set-valued Maps 

We define the inverse map of F by 

F-^ : Ran(F) ^ Vo{Qi), y ^ 7r2(Graph(F)-i n {{y} x Q^)) e Voin^). 

We say F is jro(f^)-invertible i/ 7r2(Graph(F)-i n {{y} x Qi)) e J'oi^i) for all y G Ran(F). 
If F is J^Q-invertible it holds that F~^ maps Ran(i<') into ^o(^i)- 

Theorem 1.2.17. Let F : i— > ^0(^2) be a set-valued map, then the following statements 
hold: 

(i) It holds that 

7ri(Graph(F)~^ n (y x X)) = ^2(Graph(F) n {X x Y)) 
7r2(Graph(F)-^ H {Y x X)) = ^i(Graph(F) n {X x Y)) 

for all X oni,Y and Graph(F-i) = Graph(F)-i. 

(a) If F is J-Q-invertible, then F^^ is upper semi- continuous if and only if Fa is closed 
whenever A is closed. 

(Hi) If F is upper semi- continuous, then F is J^Q-invertible and F~^ satisfies that {F^^)a is 
closed whenever A is closed. 

(iv) If F is upper semi- continuous and locally bounded, then F is J^Q-invertible and F^^ is 
proper. 

(v) If F is proper, then F is J^Q-invertible and F~^ is upper semi- continuous and locally 
bounded. 

Proof, (i) Observe that 

7ri(Graph(F)-i n {Y x X)) = 7:i{{{y,x) e Q2 x ^1 \ {x,y) G Graph(F) n (X x Y)]) 

= 7r,(Graph(F) n (X x Y)), i,j G {l,2},i 

for all X C and y C Jl2- 

If {y,x) £ Graph(F~i) it holds that 

X G 7r2(Graph(F)-i n {{y} x Qi)), 

thus {y,x) G Graph(F)^-'^ implying {x,y) G Graph(F). 

If ix,y) G Graph(F) we have y G 7r2(Graph(F) n {x} x O2) = 7ri(Graph(F)"^ x {x}), 
thus X G 7r2(Graph(F)~^ n {y} x n2) = (F^^ and {y,x) G Graph(F)"^ 

(ii) Since F is ^o-hivertible it follows by (i) that 

7ri(Graph(F-^) n {Q2 x A)) = ^2(Graph(F) n A x Q^) = Fa 

and Proposition 11.2.12] yields that is upper semi-continuous if and only if Fa is closed 
whenever A is closed. 
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(iii) Pr op osition 1 1 ■ 2 ■ 1 2] yields that Fa = 7ri(Graph(F) n {Qi x ^4)) is closed, whenever A <^ Q2 
is closed. By (i) we conclude that 7r2(Graph(F)"^ n (A x Qi)) = 7ri(Graph(F) n {0,2 x A)) 
is closed whenever A is closed. We obtain {F~^)y = 7r2(Graph(F)^-'^ n {{y} x f^i)) G Tq{Qi) 
for all y £ Ran(^), so A is .Fo-invertible. Furthermore {F~^)a is closed whenever A C is 
closed. 

(iv) By (iii) we obtain that F is .Fo-invertible and since 7r2(Graph(F^"^)n(-firnf]2)) (s x 
for all K (S r^i by Proposition 1 1 . 2 . lU] we follow by (i) that 

7ri(Graph(F-i) n (1^2 x K)) = 7r2(Graph(F-i) n {K n ^2)) x ^2, 

so F^^ is proper. 

(v) If F is proper, then we have 7ri(Graph(F) n (rJi x K)) (g ili x 1^2 for all K (g ^2, thus 

7r2(Graph(F-^) n {K x S^s)) = vri(Graph(F) n (S^i x K)) 

is compact. In particular if we set K = {y}. So {F^^)y G J^q{Q.i) and Proposition 11.2.10] 
implies that F~^ is closed and locally bounded. □ 

Example 1.2.18. Let us reconsider Example ll.2.13t where the set-valued map F : x ^ 
{/(x)} was obtained from a continuous function / : Oi — > ^22- Then the inverse map F~^ 
is defined by the inverse images of /, i.e. {F-^)y = r^{{y}). Theorem [L21Z1 yields that 
F~^ is a set-valued map such that is closed whenever A is closed. In addition F~^ 

is upper semi-continuous if and only if / is a closed function, i.e. f{A) C 0,2 is closed for all 
A (ZOi closed. 

Theorem 1.2.19. Let F : Vti ^ Tq{VL2) and G : 0, ^ J^q{Qi) be set-valued maps, such that 
one of the following properties holds: 

• G is single-valued. 

• F is upper semi- continuous and locally bounded, and G is locally bounded. 

• Fa is closed whenever A C Q is closed. 

then we define the composition map 

{FoG) :Q^J^o{n2), x^{FoG)^. 
where {F o G)x '■= iF)G^ ■ It satisfies the following properties: 

(i) If Fa is closed whenever A <^ 0,i is closed and Gb is closed whenever i? C is closed, 
then {F o G)b is closed whenever B <^ Q is closed. 

(a) If F and G are upper semi- continuous, then (F o G) is upper semi- continuous. 

(iii) If F and G are locally bounded, then (F o G) is locally bounded. 
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Proof. Since Gx S ^0(^1 ) is non-empty, there exists some y G Gx and since Fy G J^q{Q2) we 
have that {F o G)x is a non-empty subset of 172- 

If G is single- valued it immediately follows that {F o G)x = [F)g^ G ^0(^2)- 

In the case where F is upper semi-continuous it follows by Theorem 11.2.141 (ii) that F is 
closed. Since G is locally bounded it follows that Gx is compact. By Proposition 11.2. lUI (ii) 
we have that 7r2(Graph(F) r\Gx x ^2) = Uj/gG ^ compact subset of VL2. This yields 

{FoG)x&To{^2). 

In the case where Fa is closed whenever j4 C is closed, we have that Gx S J-{){^i) implies 
that {FoG)x = {F)g^ eM^2). 

(i) {F o G)b = {F)gb is closed for all i? C O closed, since Gb is closed and Fa is closed for 
arbitrary j4 C Jli closed. 

(ii) Let W be some open neighborhood of {F o G)x, then W is an open neighborhood for all 
Fy with y S Gx- Due to the upper semi-continuity of F we can find neighborhoods Zy for all 
y & Gx, such that C for 2; E Z^. Put Z := UyeG^ ■^s/' then Z is an open neighborhood 
of Gx and since G is upper semi-continuous we can find an open neighborhood X of x such 
that Gw ^ Z for all w e X. It follows that UyeG» — ^ 2/ S X, thus F o B is upper 
semi-continuous. 

(iii) Let X be some bounded subset Q, then 

y (FoG).= u U i^.= U^^ 

x&X x£X y€Gx yeC 

where G := IJ^ex Since G is locally bounded it holds that C is a bounded subset of Oi. 
Due to the locally boundedness of F we conclude that UyeC -^y ^ bounded subset of 
thus F o G is locally bounded. □ 

Remark 1.2.20. Note that if the set-valued map F : 0,i ^ ^^0(1^2) is neither upper semi- 
continuous nor closed and G : Q ^ To{Qi) is a locally bounded set-valued map which is 
not single- valued, the composition {F o G) need not define a set-valued map 0, J'^o(^2)- 
Consider the following example: Let F be defined by 



Fx 



{0} X < 
{1 + x} x>0 



Observe that it is not upper semi-continuous at and F[o,+oo] = {0}U]1, 00]. Let Gx ■= [—1, 1] 
for all X G M, then {F o G)x = {0}U]1, 2] J='o{R) for any x G M. 

Corollary 1.2.21. Let F : ^li ^ .^0(1^2) be an upper semi- continuous, locally bounded set- 
valued map. If g : Q ^ Qi is a continuous map, then the composition F o g is an upper 
semi- continuous, locally bounded set-valued map $7 J^q{Q2)- 

In the case where F is convex-valued, we obtain that F o g is again convex-valued. If H is 
the supporting function of F, we have that {x,w) — > H{g{x),w) is the supporting function of 
Fog. 

Proof. Follows by Theorem 11.2.191 and Definition 11.2.71 □ 
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For the remaining part of this Section we assume that l^i is an open or closed subset of M"*, 
such that the set of ah Lebesgue measurable subsets of is complete with respect to the 
Lebesgue measure. 

Definition 1.2.22. Let F be a set-valued map from Qi into M™. Then we call F Lebesgue 
measurable, if the set 

{x G ill I n / 0} 
is a Lebesgue measurable set for all open sets C M™. 

We denote the set of Lebesgue measurable set-valued maps by J^q{M.'^)) . 

Theorem 1.2.23. Let F : i^i ^ jro(M™) be a set-valued map. Then the following properties 
are equivalent: 

(i) F is Lebesgue measurable. 

(ii) Graph(F) belongs to Ci^i) ® ^(M"^). 

(Hi) for all w G the map x ^ infagi;'^ |a — is Lebesgue measurable. 

Proof. See [3, Theorem 8.1.4]. □ 

Proposition 1.2.24. Let F : f^i — > Tq{W^) be a set-valued map. If F is measurable, then 
its support function H has the property that x i— > H{x,w) is measurable for all w G M". 

// F is convex-valued and locally bounded, then F is measurable if and only if its supporting 
function H has the property that x ^ H{x,w) is measurable for all w G M™". 



Proof. We refer to [3 Theorem 8.2.10]. □ 

Definition 1.2.25. Let K and F G /Co(IR™)), then we define the integral over 

K by 

j Fydy := \b \ yw : (6, w) < j H{y, w) dy 

K y K , 

If F £ Ai{J; ICo{M.'^)) where J is some subinterval of W, we use the notation Frdr := 
Jj^ F-r dr and F-j- dr := — Jj^ i<V dr for all s,t £ J with t < s. 

Remark 1.2.26. Note that the integral defined above is only available for set- valued maps 
with convex values. We do not investigate its relation to the Aumann integral (see jE, Def- 
inition 8.6.1.] or [381 Definition 17.1.1]), which is available (if the domain of integration is 
compact) for any F G M.{Vli; Tq{V12)) . 

Remark 1.2.27. We define now the space of (Lebesgue) integrable functions by 

C^{VLi-lCo{VL2)) = {x^F^£ Al(Oi;/Co(J^2)) \ 3C eM -.^K (£ J, j Fydy< C}. 

K 
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This space can be equipped with a pseudo-metric 

dc{F,G):= sup [d{Fy,Gy)dy= [ sup \HF{y,w/\w\) - HG{y,w/\w\)\dy 

K fti 

where Hp,Hg denote the supporting functions oi F,G G /Co(i^2))- We obtain an 

equivalence relation by putting F ~ G if dc{F,G) = 0. By a standard argument the space 
of equivalence classes L^{^li; ICo{W^)) := {Qi; ICo{W^)) / ~ is a metric space (with metric 
dc)- 

It is straight-forward now to define /Co(i^2)) and its local version L^^^(r2i; /Co(f^2)) for 

all p gN. 

Theorem 1.2.28. Let J be some subinterval o/M. Furthermore let F £ L^^^^J; JCo{W^)) and 
t & J fixed and Cq G 1Co{W^) . Then we have that 

s 

Go + j ^rdr 
t 

defines an upper semi- continuous, locally bounded and convex-valued set-valued map from J 
to M™" called the primitive of F. Furthermore for any c £ Gq there exists an absolutely 
continuous selection f : J ^ satisfying f{s) G Cq + F-rdr with f{t) = c. It has the 
property that 

max {r,s) 

|/(s)-/(r)|< / snp H{t,w/\w\) dr. 

min {r,s) 



Proof. Let H denote the supporting function of F and let Hq be the convex supporting function 
of the set Cq. We observe that for all b G Gq + F^ dr it holds that 



\b\ < sup 



ho{w/\w\)+ / H{T,w/\w\)dT 



< oo, 



thus 



s s 

sup |Co + / -FVdrl < sup sup \hQ{w/\w\) -\- / H{ 



T, w/\w\) drl < oo 



for any bounded M C f^i. Since s ^ hQ{w) + H{t, w) dr is upper semi-continuous for fixed 
w G M*^ it follows by Theorem OH that s ^ Go + // F^ dr is upper semi-continuous. 

It is sufficient to construct a continuous selection / of Frdr with f{t) = 0. Putting 
/c(s) := f{s) -|- c for some c £ Gq immediately yields a continuous selection of Cq + F^- dr 
with fc{t) = c. So without loss of generality we may assume Cq = {0}. 

Let us first consider some interval [t,T+] (s J for some T-|- G J with > t. For /c G N we 
put Xk := and t^j := t + jXk, so t^^ = t and t^^k = T+. We obtain a decomposition of 

the interval in k subintervals [tkj~i,tkj] for 1 < j < k. 
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We are going to construct a sequence of continuous functions {fk)keN such that 

s 

{fk{s),w) < J H{T,w)dT 
t 

holds for all w eW^ and A; G N. 

Put Bf^s '■= It -^rdT. Any bt^s € has the property that {bt,s-,w) < H{t,w) dr = 
- Jl H{t, w) dr and -{bs,t, w) > - // H{t, w) dr for all w e W^, thus 

{bt,s,w) < -{bs,t,w), weW, 

which immediately implies bg^r = —br,s- It follows directly that bs^s = and Bg^s = {0} and 
Bt^s = —Bs,t for any s,t e J. Furthermore we observe that 

r s s 

J H{t, w)dT + J H{t, w)dT = J H{t, w) dr 

t r t 

for all w G M"* which implies Bf^r + -Br-,s = -Bt,s- 

Put fk{s) = for all s G [t,t + Afe], then we define iteratively 

fk |[t+jAfe,t+(j+l)Afc] (s) = fk{t + j^k) + bt+jXk,s, i = 1, • • • , ^ 

for some btJ^jx^^g £ -Bt+(j-i)A^,s- It follows immediately that 

for s G [t + /Afe,f + (Z + l)Ajt], which implies that {fk)keN is an equi-bounded family. 
Assume s, r G [t, T+j then we have 

(1.2.5) 

h h 
{fk{s) - fk{r), w) = bt+^j_i)x^^t+jXk + bt+hXk,s - XI h+{j-i)Xk,t+j>^k " ^t+hXk,r, w) 

(max(;i,^2) \ 
-&t+(mmai,i2)-l)Afe,mm(s,r) + X] ^«+0-l)Afe,t+jAfc + ^t+max (Zi ,Z2)Afe,max (s,r) I , "U^) 
j=min(ii,/2) / 

andsmce -^t+(min(h,i2)-l)Afe,min(6,r) + Ej=min(li'!z2) ^*+jAfc,t+jAfe +^t+max (h,b)Afc ^ ^n^jn (^_5)_max (r,s) 

we have the estimate 

max (r,s) 



min (r,s) 

implying 



max (r,s) 

\h{s)- h{r)\< [ sup H{T,w/\w\)dT 

I (r-,s) 
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which show that {fk)keN is equi-continuous family of continuous functions. 

By the theorem of Arzela-AscoU we can find a convergent subsequence {fkj)jeN converging 
to some continuous function / G C([t,T+]). It has the property that 



s 

{fis),^) = lim {fkAs),w) < / H{T,w)dT, 



t 

so / is a continuous selection of the set-valued map s i— Ft dr for all s G [t, T~^]. Absolute 
continuity of / follows by 

max (r,s) 

\f{s)-fir)\ < \fis)-fkM\+\fkM-fkM+\fkA^)-fi^)\ I H{T,vj/\w\)dr. 

min (r,s) 

The proof for the backward direction [T_,i] d J for some & J with r_ < i is analogous. 
In the case of a non-compact interval J the extension of the selection to the whole interval is 
straight-forward. □ 
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2 The generalized graph 



In this section we introduce the concept of a generahzed graph for a c-bounded generahzed 
map. The generahzed graph extends the classical graph of a continuous map, in the sense 
that the generalized graph of the embedded map (which is a c-bounded generalized map) 
coincides with the classical graph. Furthermore the generalized graph is closed in the product 
topology of the domain and the image space of the generalized map. 

Throughout the entire Section we let i}i,0,2 be open subsets of M" resp. W^. We define the 
projection maps of the product space of f^i x onto f^i resp. by 

TTi -.0,1 X ^ r^i, {x, y) ^ X resp. 7r2 : l^i x ^ {x, y) ^ y- 



2.1 The generalized graph of a c-bounded Colombeau map 

Definition 2.1.1. Let F G ^[r2i,rJ2] be a c-bounded generalized map. Then the set 
Graph(F) := ^{x,y) £ Qi x Q2 \ ^ <i net {xe)e in ^1 '■ lirnXe = x G ^li 

(2.1.1) 

and y G CP {{F^(^i,-^{xe))e) for some map r £ T} 

is called the generalized graph of F . By an abuse of notation we use the symbol Graph(F) 
from Definition \1.2.5\ for the generalized graph of a Colombeau map. This is motivated by 
the observation made further below ( cf. Theorem \2.1.8\) that the generalized graph actually 
defines a set-valued map by x ^ 7r2(Graph(F) n {{x} x ^2))- 

Observe that the definition for the generalized graph does not depend on the representative 
of F: Let (Fg)^, (-Fe)^ be two representatives of F, then := F^ — Fgr is negligible and 

\PT{e){Xe) - FT{e){Xe)\ < SUp \nr(s){z)\ = 0{t{s)) 

as £ ^ (since t £ T), where X is some compact set containing the converging net x^ for 
small £. 

We use the notation Graph(F)^ := 7r2(Graph(F) n {x} x $72), thus 

Graph(F) = |J {x} x Graph(F)^. 

Definition 2.1.2. Let F £ 0,2] be a c-bounded generalized map. We say F is of natural 

type, if it has a representative such that 

is continuous on]0, 1] for all x G Qi fixed. Note that a Colombeau map which is the embedding 
of a distribution is of natural type. 
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Example 2.1.3. Consider the c-bounded Colombeau function F of natural type defined by 
the representative (-Fe)e with 

Fs := sin(x/e), 
then it is straight-forward to calculate its generalized graph 

Graph(F) = M x [-1,1]. 

Example 2.1.4. Consider the Colombeau function F := l{H), where H is the Heaviside 
function. If we set g := H * p, where p is the moUifier of the embedding l with J p{y) dy = 1 
and J y^p{y)dy = for all A; > 1, we have that -Fe(x) := g{x/e) defines a representative of 
F. For all nets (xe)^ tending to xq / 0, we obtain that F^{xi;) — > sign(xo) as e ^ 0. Let us 
consider the nets {xe)e tending to 0: Since we can find zero- nets (xe)^ such that {x^/e) tends 
to any point in R U {±00} and g is continuous and bounded, we have that 

Graph(F) := ({x G M | x < 0} x {0}) U ({0} x Ran(5)) U ({x G M | x > 0} x {1}) . 



Note that Ran(g) is a closed interval that contains [0, 1]. 

Definition 2.1.5. Let F G ^[r2i,ri2] be a c-bounded generalized map. F is said to be equi- 
continuous at xq, if there exists a representative (-Fe)e which is equi- continuous in xq: For all 
7 > 0, there exists some 5 > such that 

|F,(x)-F,(xo)| <7 

holds for all x G Bs{xq) and e g]0, 1]. The generalized map F is called equi-continuous on 
some subset X C f^i, if it has a representative which is equi-continuous in each point x £ X. 

Proposition 2.1.6. Let F G ^[r2i,ri2] be a c-bounded generalized map. If F is equi- 
continuous on Qi, then it follows that the generalized graph of F can be determined pointwise 
by 

(2.1.2) Graph(F) = {(x, y) G J^i x I y G CP ((F,(x)),) } . 

Proof. We can choose a equi-continuous representative {Fi;)^ and proceed along a standard 
argument: Let {xs)e, {x'^)e be two nets tending to x G ili. Then we have, using the local 
equi-continuity of (F^)^ on Oi, that for all 7 > we can find some 6 > 0, such that the 
distance is bounded by 

|F,(,)(x,) -F,(,)(x;)| < 7 

for all X£,x'j. G Bs/2{x) uniformly for all t £T. Since lim^^QXe = lim^^Qx'^ = x there exists 
e' g]0, 1] such that |xe — x'^\ < S holds for all e < e' . Thus \m\f.^Q{F^{^^^^{xs) — F^(£)(x^)) = 
for any r G T and it suffices to consider the constant net x^ := x in Definition 12.1.11 As 
T G T we have that (F^(j)(x))e is a subnet of (Fir(x)\ and 

CP((F,(,)(x)),) CCP((F,(x)),), 

which yields the statement ()2.1.2p . □ 
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2.1 The generalized graph of a c-bounded Colombeau map 

Lemma 2.1.7. Let F G ^[J7i,ri2] be a Colombeau generalized map. If {x,y) E (f^i x 
r22)\ Graph(F), then there exists neighborhoods Y Q2 of y, X Q Qi of x and e' G]0, 1] 
such that 

YnFe{X) =0 

for all e < e' . 

Proof. Let {x,y) £ (ili x r22)\ Graph(F). The proof proceeds by contradiction: Assume that 
for ah neighborhoods Y' ^ of y and X' C Oi of x and for ah e' > there exists some 
T-.e'^ t{X',Y',£') < e' with 

y'nF^(,,)(X')^0. 

By setting X' := Bf;{x) n f^i, Y' := Bf.(y) n Q.2 and e' := e we obtain a map r : e 1— > 
t{Bs{x) n ^i,B^{y) n 02,e). As T(e) < e for all e G]0, 1] we have t £ T. We choose a net 
Xe G Be{x) n satisfying G ^^(y) n 1^2 for all e g]0, 1]. 

Obviously {xe)e is a convergent net in f^i with lime_^o3;£ = x G f^i and the net 
converges to y. Since t £ T we obtain y G CP ((^^(^■((xe))^), so (x,y) G Graph(F) contra- 
dicting the assumption. □ 

Theorem 2.1.8 (Properties of the generalized graph). Let F G ^[r2i,il2] 0, c-bounded 
generalized map, then the generalized graph Graph(F) has the following properties: 

(i) It extends the classical notion of a graph (of a continuous function) in the following 
sense: If F is the embedding of a continuous function G := {Gi, ...,Gm) £ G{Qi,i}2), 
i.e. F = {l{Gi), L{Gm)), then the generalized graph of F coincides with the graph of 
the continuous function G. 

(a) It defines a closed and locally bounded set-valued map by 

r^i J^o{^2)i X ^ Graph(F)a;. 

Proof, (i) follows by Proposition I2.1.6] and the fact that the embedding of a continuous func- 
tion is an equi-continuous generalized function. 

(ii) Graph(F)a; 7^ for all x G r^i is trivial. 

We proceed showing that Graph(F) is closed in Cli x 0,2- Let (x,y) G (Oi x il2)\ Graph(F), 
then we know by Lemma 12.1.71 that there exists X Q 0,i,Y <^ ^}2 open neighborhoods of x 
resp. y and e' > 0, such that 

Fe{X)nY = 

for all £ < e' . Now for any point y € X and any net {ye)e in ^1 converging to y, there exists 
some e" > such that ye £ X for all e < e" . It follows that 

Graph(F) n (X x y) = 0, 

so (ill X ^2)/ Graph(F) is open in i7i x Q.2 and hence Graph(F) is closed in the product 
topology of VLi xVL2. 
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2 The generalized graph 



In the next step we proove that x Giaph{F)x is locally bounded. Let i? C be bounded, 
then there exists some K ^ Qi with B <^ K. Furthermore let Ki be some compact neighbor- 
hood for all points in K. Since F is c-bounded we can find some e' > and K2 (s ^2, such 
that 

Fe{Ki) C K2 

for all £ < e' . So Graph(F) n K x C Graph(F) n -fC x K2, which implies that Graph(F) n 
iiT X is a bounded subset of ili x ^22- So x 1— > Graph(F)j; is locally bounded. 

Putting K = {x} we obtain that Graph{F)x = 7r2(Graph(F) n ({x} x ^2)) is closed, since it 
is the continuous image of a bounded, closed set. So Graph(F)a; G J^q{Q2) for all x £ Qi. 

□ 

Corollary 2.1.9. Since the generalized graph can be seen as a closed and locally bounded set- 
valued map X ^ Graph(F)a;, it follows by Theorem \1. 2. 14\ that it is upper semi- continuous, 



i.e. for any open neighboorhood W of GTaph.{F)x, there exists an open neighboorhood X of 
X, such that Graph(F)x C W. 

Proposition 2.1.10. Let F £ Q[Qi,Q2] be a c-bounded generalized map of natural type, then 
the generalized graph Graph(i<') has the property that Graph(F)^ is connected for all x G M". 



Proof. We prove by contradiction: Assume that Graph(F)a; is not connected, then there exist 
V, W open subsets of M"^ with V := V D Graph(F)^, W := W D Graph(F)^ / such that 
VnW = and VUW = Graph(F)^.. Choose a £V,b gW, then we can find a net (xe)^, {x'^)s 
with Xs,x'^ X and t,t' G T such that j/e := F^(£)(x£) a and y'^ := F^n^^-^{x'^) — > h. Since 
y, W are neighboorhoods of a resp. b there exists some e G]0, 1] such that 

for all e G]0,e]. Observe that since F is of natural type we have that 

Qe-a^ ^ar(e)+{l-Q)T'(£)("a:e + (1 - a)x'^) 

is a continuous map with gd^) = Ue ^ind ^^(l) = y^. The intermediate value theorem then 
yields that 5e([0, 1]) is connected for fixed e. Since V r\W = ^ and for e g]0, e] we have that 
V r\ ^£([0, 1]), VFn5e([0, l]) are non-empty, open sets in the relative topology of Ran(5£) C VL2 
it follows that (V U W) n 5e([0, 1]) / ^£([0, 1]) (equality would contradict the connectedness 
of ^^([0, 1])), thus there exists some G [0, 1] with geiae) G U Wy. 

Putting ^{e) := aeT{e) + (1 — Qe)T'(e) and y^ ■= o-eXe -|- (1 — a^^x'^ we observe that /U G T 
and ye ^ x. Since g^ia,) = F^(^,){y,) G (F U Wf it follows that / CP {{F^ie){ye))e) ^ 
{V U wy n Graph(F)^. which is a contradiction to {VUW)n Graph(F)^. = Graph(F)a,. □ 

Remark 2.1.11. If F G ^[r2i,M] is of natural type, then Proposition 12.1.1^ yields that 
Graph(i<')i. is a convex set for all x G M". 

For the higher dimensional case (m > 1), the convexity of Graph(F)2. does not hold for 
arbitrary F G ^[J7i,M'"] of natural type. Consider the Colombeau map in ^[M,M^] defined 

by 

Fs{x) := (sin^/^(x/e),cos^/^(x/e)), 
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2.2 The generalized graph as a set-valued map 



then it is straight-forward to verify Graph(F)o = Uae[-i i] 0) U (0, a)), which is connected 
but not convex. 

If F £ ^[Oi,M] is not of natural type, then Graph(i^)i; is in general neither convex nor 
connected. Just consider the simple Colombeau map F € Q\M., M] defined by the representative 



Its generalized graph satisfies Graph(F)2; = {0} U {1} for all x G M. 

Lemma 2.1.12. Let F G ^[ili,02] be a c-bounded generalized map. Furthermore, let Ki (s 
r^i and K2 := 7r2(Graph(F) H {Ki x ^2))- IfY is a neighborhood of K2, then there exists a 
neighborhood X O fli of Ki and e' g]0, 1], such that 

Fe{X) C Y 

for all e < e' . 



Proof. First we proof the Lemma for Ki := {x} containing only a single point x G ili. We 
note that K2 is a compact subset of Q2 by Theorem 12. 1.81 Assume that the statement of the 
Lemma does not hold, that is for all neighborhoods X' of x and for all e' g]0, 1] there exists 
some r(e',X') < e' and Xg/ G X' with 

Then (by setting X' = Bi;{x) n ili and e' = e) we define r : e 1-^ T(e, B^{x) n ^i) and obtain 
some net Xg G B^{x) n f^i such that F^(£)(xe) y. Note that T(e) < e for all e g]0, 1] implies 
T G T. 

Since y is a neighborhood of K2, we have that n K2 = 0. The net (xg)^ converges to x 
and the net (^^(^■((xg))^ G has the property that the set of clusterpoints CP (xe))^) 
is contained in the closure of Y'^. By definition of the generalized graph Graph(F) we have 
that CP ((F^(£)(xe))£) C K2, which contradicts the fact that Y^r\K2 = %. 

Now we consider the general case, when Ki is an arbitrary compact set. Then y is a 
neighborhood of vr2(C/ H Ki x VL2) and it follows that y is a neighborhood of each set 
Ti2{{z} X Graph(F)2) for z G Ki. We can apply the first part of the proof for each point 
z and obtain (open) neighborhoods Xz of each point z such that F^{Xz) C Z for all e < e^, 
where Ez G]0, 1] depends on z. Since Ki is compact and {Xz)z^Kx is an open covering we can 
choose some finite subcovering {Xz^)\^i. Then X := IJfc=i ^z^. is a neighborhood of Ki such 
that feiX) C y holds for e < e' := min^^^^ e^^. □ 



2.2 The generalized graph as a set-valued map 



By Theorem l2.1.8l it turns out that we think of the generalized graph of a c-bounded Colombeau 
map F G ^[J7i,r22] as a set-valued map 



J^o(f^2), X ^ Graph(F) 
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2 The generalized graph 



Since the generahzed graph Graph(F) is closed in ili x and Graph{F)x is compact for all 
X E ill it follows by Proposition 11.2.10] and Theorem 1 1 . 2 . 1 4 1 that x i-^ Graph(i^)2^ is an upper 
semi-continuous and locally bounded set-valued map. 

It turns out that the convex hull of the generalized graph is an important tool to characterize 
the shadows of Golombeau solutions of ordinary differential equations, where the right-hand 
side A := [(^e)^] is a Golombeau map satisfying certain bounds. 

Theorem 2.2.1. Let F € ^[r2i,rJ2] and G G then the composition FoG G t/[r2,il2] 

defined by {F^ o Gfr)e has the property that 

Graph(F o G)^ C (Graph(F) o Graph(G))^ for all x £ Q, 

where the composition on the right-hand side denotes the composition of set-valued maps ( cf. 
Theorem \1.2.19\) . 



Proof. Note that (Graph(F)oGraph(G))2: is again an upper semi-continuous, locally bounded 
set- valued map due to Theorem ll.2.19[ 

Let a G Graph(F o G)x then there exists some t £ T and a net (xe)^ with — > x, such that 

a = liniF^(^)(G^(£)(xe)). 

Without loss of generality we can assume that = GT-(£)(xe) converges to some y G Graph(G)i.. 
It immediately follows that a G Graph(F)y and we have Graph(FoG)a; C UyeGraph(G)a; Graph(F)y 
(Graph(F) o Graph(G))^. " □ 

Proposition 2.2.2. Assume F G Q{0,i,Q2) is a c-bounded Golombeau map with representa- 
tive (i^e)e. Then the supporting function of the generalized graph Graph(F) defined by 

Hf{x,w) := sup {a,w) 

aGGraph{F)x 

can be obtained by 

(2.2.3) Hf{x,w) = \im sup {Ff,{y),w). 

'5^o?/eB5(x)nni,Ate]o,<5] 



Proof. Let w G be fixed. We put g{ii,y) := {Ffj_{y),w) for fixed u; G M" and Ms{x) := 
]0,5] X Bs(x) (which is a descending family of sets), then Lemma lA.51 yields 

hs{x,w) := sup g{fi,y) = lim g{fis,'Y,X5,^) 
for some converging net (/x^,^, X5^^)^g]o,i] in M^{x). 

For all /?, (5 > there exists a Ti{p,5) G]0, (5] such that \hs{x.,w) — g{lJ'y,STyy,s)\ < /o/2 for all 

Since the right-hand side of (|2.2.3|) is defined by 

h{x,uj) = lim hs{x,w) 
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2.2 The generalized graph as a set-valued map 



there exists for all p > some T2{p) £]0,p] such that \h{x,w) — hs{x,w)\ < p/2 for all 
S < T2{p)- Setting ip{e) := /i^^^^ .j.^^^)) ^-^(g) (which implies ?/;(e) < e since py^s ^ S, thus tp £ T) 
and Xs ■= ?/ri(e.r2(£)),T2(e) ^e obtain the estimate 

\h{x,w) - {F^(^^){xe),w)\ < \h{x,w) - hr2{e)(.x,w)\ + |/i^2(s)(x,u;) - g{'il){e),Xe))\ 

= \h{x,w)-hr2{e){x,w)\ + \hr.,(e){x,w)-g{pri{sMe))Me)^^Ti{s,T2{e))Me))\ ^ e/2 + e/2 =6. 

We conclude u)) = lim£^o(-^v(£)(^£)' Note that any clusterpoint a G CP (x£))e) C 

Graph(F)a; has the property that h{x, w) = {a, w) < Hp{x, w) for w £ M™ fixed. By repeating 
the steps of the proof for all w G M™ the bound 

h{x, w) < Hp{x, w) 

is obtained for all x G $7i and w G M™. 

It remains to prove Hpit^x) < h{t,x): Since Hf{x,w) = {aw,w) for fixed x G M", where 
ttw = lim£_>o -FT-(e)(2;e) for some map t £ T and x^ x. Without loss of generality we assume 
t{£) < S- It follows that 

Hf{x,w) = lim{F^^^^{xe),w). 

For all 5 > we can find some e'{6) £]0,6] such that \xe — x\ < 6 for e £]0,s'{6)] and 
T(e) <e< e'{6) < 6 for all e £]0,e'{6)]. It follows that 

HFix,w) = lim{F^r^){xe),w) = lim sup {F^M{xe),w) 

< lim sup {F^{y),w) = \\Tahs{x^w) = h{x,w). 

□ 

Example 2.2.3. Recall Example l2.1.4l where we considered the Colombeau embedding of the 
Heaviside function F = l{Q). If we put g := Q * p, where p is the mollifier of the embedding 
i with f p{y)dy = 1 and f y^p{y)dy = for all A; > 1, we have that := g{x/e) is 

a representative of F. Note that limy^_oo5(y) = and limy^+oo = 1- According to 
Proposition 12.2.21 the supporting function Hp can be obtained by 



Hf{x,+1) = lim sup g{y/p) = lim sup g[z) 

'5^Oj/e-B^(x0n!^i,AtG]O,5] '5^0 2e<5-i-Bi(x) 



f x > 

suPj/gM 5(y) a; = 
1 X < 



and 

Hf{x, -1) = lini sup -g{y/p) = - lim inf o(z) = 

'5^0'(/eB4(a;)nni,/ie]0,5] <5^0 2e<5 ^■B!,{x) 

Since F G ^[IR,M] (space dimension m = 1) it follows that Graph(F)2,. is convex for all x G 
thus the supporting function Hf generates the graph by 

Graph(F) = {a G M | Vw G M : (a, < Hf{x,w)] 




= ({x G M I X < 0} X {0}) U (^{0} X lm{g) ) U ({x G M | x > 0} x {!}) 
which is consistent with the result of Example 12.1.41 
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2 The generalized graph 



Theorem 2.2.4. Let F £ Q{ J x M")" be a Colomheau map with a representative {F^)^: such 
that there exists a positive [5 G -^ioc(^) supj^g]g,i \Ffr{t,y)\ < f3(t) for t G J. Then the 

generalized graph Graph(F(t, •)) exists for almost all t £ J fixed. Its supporting function 
defined by 

Hpit, x,w) = \im sup {F^{t,y),w), 

for almost all t £ J satisfies the following properties: 

(i) 1 1— > Hp{t,x,w) is Lebesgue measurable on J for all x,w £ R", 

(a) X Hp^t, x,w) is a upper semi- continuous for almost all t € J and w € M", 

(Hi) there exists a positive function j3 G L\^^{J) such that sup^jgjgn \Hp{t,x,w /\w\)\ < f3{t) 
for almost all t £ J and all w £ M". 

Proof. Observe that supj^g]jn [^^(t, y)| < f3{t) implies that x i— > Fg{t,x) defines a c-bounded 
Colombeau map for almost all t £ J, so the generalized graph Graph(i<'(t, •)) is defined for 
almost all t £ J. Thus property (ii) follows since x i— > Graph(F(t, x)) is a closed, locally 
bounded set-valued map for almost all t £ J, implying that its supporting function x i— > 
HF{t,x,w) is upper semi-continuos. 

Proposition I2.2.Z] immediatelv yields that the supporting function of Graph(F(t, •)) can be 
obtained by 

Hpit, x,w) = \im sup {F^{t,y),w). 

Property (iii) then follows from the estimate 

sup {F^{t,y),w/\w\) < sup \F^{t,y)\<f3{t) 

for almost all t £ J. 

Due to Lemma[A2]we have supyg5^(^)^^g]o,5] y)' = ^^^■y^o{^fj.5,^it, xs,'r),w) for some 
sequence (M(5,7) ^(5,7)7e]o,i] '^i^^ jj-s^^ G]0,5] and xs^-y £ Bs{x). 

It follows that supj^g^^(^)^^g]Q 5] {F^{t, y),w) is Lebesgue measurable as the pointwise limit of a 
net of Lebesgue measurable functions. Finally the same argument yields that t ^ HF{t, x, w) 
is Lebesgue measurable. □ 

Theorem 2.2.5 (Approximation property of the supporting function Hp). Let F £ Q[J x 

R^jM") be a Colombeau map as in Theorem \2.2.4\ then Hp denotes the supporting function of 
Graph(F(t, •)) for almost all t £ J. If {^e)e is a net of smooth functions converging uniformly 
on compact sets to some ^ £ C{J)'^ , then 

\im.sw£>{Fe{s,ie{s)),w) < Hp{s,C{s),w) 
for all w £ R™ and almost all s £ J. 
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2.2 The generalized graph as a set-valued map 



Proof. First of all we observe that 

is a family of functions in Ll^^i"^) and 

{F,{s,Us)),w) < sup \F,{s,y)\ < P{s) 
2/eK" 

for almost all s G J and e g]0, 1]. Thus 

limsup(F£(s,^£(s)),w;) 

exists for almost all s G J. 

Let K (E J he fixed, then there exists some Jl e T (without loss of generality we assume 
'jl{e) < e) such that sup^g^^ — ^{s)\ < e for all ^ g]0, ^(e)]. We have that 

sup {F^{s,^^{s)),w)\ 

VMe]0,M(e)] J ^ 

is a subnet of the net (sup^gjo^g] {F^,{s, ^^(s)), 'u;))^, which converges for almost all s G if, thus 



lini sup (F^(s,^^(s)),u;) = lini sup (F^(s, ^^(s)), iw). 

We conclude 



limsup(F£(s,^£(s)),w;) = lini sup (F^(s, ^^(s)), = lini sup {Fn{s,^^{s)),w) 
< lini sup (F^(s,y),u;) < lini sup {F^{s,y),w) = Hf{s,(,{s),w) 

for almost all s ^ K. □ 
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3 Generalized solutions of ordinary differential 
equations 



Throughout this chapter J will be some subinterval of M. Starting point of our overview of 
generalized solution concepts for ordinary differential equations is the initial value problem 

^(s) = a{s,^{s)) a.e., ^{t) = x 

where the right-hand side a{t,x) = {ai{t,x), . . . ,an(t,x)) is a map from J x M" to M". 
Occasionally we will indicate dependence on the initial conditions by writing for a 

solution ^ with ^{t) = x. 

Let us first consider the well-posedness: The left-hand side of the equation requires a possible 
solution ^ to be at least absolutely continuous in order to ensure differentiability almost 
everywhere on J. Furthermore it is apparent that the right-hand side of the equation has to 
be in Lj^^(J)"'. Nevertheless in general it is not sufficient to have a G L\^^{J x M")" in order 
to have the composition s ^ a{s,f{s)) in L\^^{J)'^ for / G j4C(J)". 

In the first section we investigate the classical Caratheodory theory, which requires the right- 
hand side a to be in the space LJ^^( J, with a{-,x) G L\^^{JT for aU x E M" 
and a(s, •) E C{J)'^ for almost all s E J (referred to as Caratheodory conditions), then the 
existence of a solution ^ E AC{J)'^ can be proved. 

Of course the continuity required in the space-variable of the right-hand side restricts the 
applicability of the Caratheodory theory. Nevertheless if we want to drop this requirement 
we cannot retain on the classical concept of a solution. 

A natural generalization of an ordinary differential equation is the differential inclusion 

C(s) E a.e., ^{t) = X, 

where At^x is a subset of M" for each (t, x) E J x M". The map A : {t, x) i— > At^x is called a set- 
valued map (cf. Section ri.2p . As in the case of Caratheodory theory we have to impose certain 
conditions (Filippov conditions) on the set- valued map A in order to obtain an existence result 
for the differential inclusion. These results will be presented in the second section. 

The concept of a differential inclusion can then be utilized to obtain a generalized solution 
concept for ordinary differential equations, when the right-hand side fails continuity in the 
space variable, but is still in LJ^^(J; L°°(M"'))". It is obtained by mapping the right-hand 
side a to some A in the space of set-valued functions satisfying the above-mentioned Filip- 
pov conditions in a way that At^x = {a{t,x)} whenever a fulfills the classical Caratheodory 
conditions. Then the absolutely continuous solution of the differential inclusion is called a 
Filippov generalized solution. 

Another approach to generalize the classical solution concept for ordinary differential equa- 
tions is the regularization of the right-hand side a by means of convolution with a net of 
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3 Generalized solutions of ordinary differential equations 

smooth mollifiers, such that the resulting net of smooth functions (Aj)^ has the property 
hm£^o(^£) v) = ('^1 v) fo'^ ^-ll ^ I''(M"). Then a net of solutions {^e)e can be obtained by 
solving the differential equation for fixed e G]0, 1] using the classical existence results. If the 
net of solutions {Ce)e possesses a subnet converging to some ^ G AC{J)"', i.e. there exists 
some T £ T such that Ct(£)(*) ~^ Ci^) uniformly on compact subsets of J, then the function ^ 
is a generalized solution in the sense that it satisfies the classical ordinary differential equation 
whenever the initial right-hand side a satisfies the Caratheodory conditions. 

It is remarkable that the latter approach is closely related to the former, since the generalized 
solution obtained by the regularization procedure satisfies a certain differential inclusion. 
The concept of a generalized graph (cf. Chapter [2]) plays a crucial role in understanding this 
relation. 

Another viewpoint of the regularization approach is that the net (A^),, can be considered 
as an element in the Colombeau algebra of generalized functions. We put A := [(^g)^] G 
g{J X M")". Note that in general A = {Ai, . . . , An) (i(ai), • • • , t(a„)), where t : P'(M") 
Q (W^) denotes the standard embedding, since we did not specify the net of smooth mollifier 
any further. Occasionally it may be convenient to use mollifiers with non-vanishing moments 
(f.e. positive or compactly supported mollifiers), which is not possible when using i to embed 
the right-hand side a. 

If the ordinary differential equation 



can be solved in the setting of Colombeau generalized functions by some ^ := € Q (J)", 

we call ^ a Colombeau solution with inital condition {t,x) £ J x M". 

It generalizes the classical solution concept whenever there exists some t G T such that 
('^r(e))e couvcrges to a function ( G ylC(J)" uniformly on compact sets of J. We say ^ has an 
absolutely continuous sub-shadow C- This notion is justified by the fact that C is an absolutely 
continuous shadow of ^, if {£,e)e converges to C uniformly on compact subsets of J. 

C is a generalized solution in the sense that it satisfies the classical ordinary differential 
equation whenever the initial right-hand side fulfills the Caratheodory conditions. 

3.1 Solutions of ordinary differential equations according to 
Caratheodory 

In this section we will review the classical Caratheodory theory for ordinary differential equa- 
tions 



with non-smooth right-hand side a : J x — > M" or equivalently the corresponding integral 
equation 



c{s) = A{s,as)) m = x 






s 



(3.1.2) 




t 
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3.1 Solutions of ordinary differential equations according to Caratheodory 



We follow closely the introduction from the book [15j. Throughout the entire section the 
right-hand side a = (ai, . . . , a„) satisfies the conditions 

(i) a{t,x) is Lebesgue measurable in t for all fixed x G M", 

(ii) a{t, x) is continuous in x for almost all t ^ J, and 

(iii) sup^jgign |a(t, x)| < f3{t) almost everywhere for some positive function /? € L^^^{J), 
referred to as Caratheodory conditions (CC). 

Note that the first two Caratheodory conditions ensure Lebesgue measurability of the com- 
position s a{s,f{s)) for all / G C{J)'^ (due to Lemma \A.9\\ . while the third condition is 
crucial in the existence proof. 

Theorem 3.1.1 (Existence theorem for ordinary differential equations). Let J be some subin- 
terval o/M and assume that a = (oi, . . . ,a„) satisfies (CC) on J x M". Let {t,x) G J x M^, 
then there exists an absolutely continuous solution = (,^1, . . . to the ordinary differential 
equation 

^(s) = a(s,^(s)) a.e., C{t) = x 



Proof. For sake of completeness we include this classical proof, following Filippov ([15, Ch. 
l,Thm. 1]). Without loss of generality we assume that t lies in the interior of J. 

Note that the Caratheodory condition (i) and (ii) imply (due to Lemma IA.9P that for any 
/ G C(J)"' the composition s 1— > a(s,/(s)) is in L\^^{J). Let us consider a compact subset 
of J in the forward direction, i.e. set := J f] [t,T^] for some T+ G J with > t. For 



any G N we put : = 



T+-t 



in order to obtain a decomposition of the interval in k 



subintervals by := [t + {I — l)Afc, t -|- ZA/c] for 1 < I < k, such that 

J+ := U J+. 

1=1 

We construct a function G AC{J^) inductively by setting 

max (t.s— Afe) 

(^kis):=x+ J a(r Afc,^fc(r))dr 
t 

for all s G J^. Whenever s G Jj^ for I > 2, the integral expression on the right-hand side of 
the above definition depends only on restricted to Jj^i. For s G J^^ we have ^fc(s) = x. 

Caratheodory condition (iii) implies 



Ms)-Ur)\ 



max (t,s— Afe) 



air + \k,(,k{r))dT 



max (t,r— Afe) 



max {r,s) 

< I /3(T)dr 

mill {r,s) 

T+ 



for s,r > t 



and Cfc(s) G Bp{x) for all s G J"*", where p := ^ /3(r) dr. Obviously {^k)k is a equi- 
continuous and equi-bounded family in C{J^)'^. The Theorem of Arzela-Ascoli yields an 

uniformly convergent subsequence ^fc, >° ^, where ^ G C(J''")" with ^(t) = x. 
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Observe that 



s— At. . 

'^3 



holds. So we finally obtain by the theorem of dominated convergence that 

niax{t,s-Xk.) max(t,s) 

^(s) =x + hm / a(r + Ajt . , . (r)) dr = x + hm / a{T, ^fe . (r - A^ . )) dr 
t t 

max (t,s) max (t.s) 

= x+ J \im^a{T,^kj{r - \kj))<ir = x+ J o(T,^(T))dT. 
t t 

Equivalently ^ satisfies the ordinary differential equation ^(s) = a(s, ^(s)), = 2; for almost 
all s G J+. The proof for backward problem on some J~ := Jn [T~, t] for some T~ G J with 
r~ < t is analogously. 

In the case where J is an open interval or unbounded, it is straight-forward to decompose the 
domain in countable many compact subintervals. Then the local solutions for each subinter- 
val are continuously glued together (by choosing appropriate initial conditions for the local 
solutions) to obtain a global solution. □ 

Definition 3.1.2. Let Q, be a non-empty subset of J x M", then we define 

E{n) := G AC{Jr I x)en: ^s) = a(s,e(s)) a.e., ^ = ^} 

denoted as ( Caratheodory) solution set with initial conditions in If we put Ef^x '■= '^({{t, x)}) 
it obviously holds that 

The solution set E](J7) is called forward (resp. backward) unique, if ^1,^2 € E(Q,) with ^i{t) = 
^2{t) implies = ^2{s) for all s >t (resp. s <t). 

In the case where J = R and Q' is a non-empty subset of M", we refer to a solution set 
E{R X fl') as autonomous, if ^ E E{R x il') implies + r) G E{R x Q') for any r G M. 

Theorem 3.1.3 (Properties of the solution set). Let be a non-empty subset o/M", then 
the solution set E{J x fl) has the following properties: 

(i) For every {t, x) E J x Q, there exists $, G E{J x Q) with ^{t) = x. 

(ii) If t £ J and ft is closed, then E{{t} x Q) is a closed subset of C(Jy''. 
(Hi) (0^eE:(Jxf2) equi- continuous family in C{J)"'. 

(iv) IfQm M", then E{{t} x 9.) is a compact subset of C{J)'^. 
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If the right-hand side is time-independent, we put J = M and the solution set H(M x Q.) is 
autonomous. 

Proof, (i) By Theorem 13.1.11 the set H(J x J]) is a non-empty subset of C{J)", such that for 
every {t,x) G J x Q there exists a ^ G H( J x fi) with ^(t) = x. 

(ii) First we prove that x i7) is a closed set in C{J)^ if Q is cfosed: Let {S,k)keN be 
sequence in x fi) uniformly converging (on compact sets) to some ^ G C{J)"'. In 
particular Xk = ^ki^) (,{t) =: x £ ^1, since il. is closed. It follows that 

^(s) = lim ^fe(s) = lim \xk+ a{T,(,k{T),w) dr 

V i y 

= x+ lim a(T,^fc(r))dr = x + / a(r,^(r))dr, 
y fc^oo J 
t t 

which implies G x i7) satisfying the initial condition ^(t) = x G il, so x i7) is a 

closed subset of C{J)"-. 

(iii) Using the integral equation ()3.1.2|) we have for all s,t € J that 

max {t,s) 

(3.1.3) \^(s)-m\< J 

min (tys) 

holds uniformly for all ^ G H(J x fi), thus H(J x fi) is an equi-continuous family of continuous 
functions. 

(iv) Since K (s M" the estimate 

max {t,s) 



sup |'^(s)| < sup |x| + sup 



/a(r, ^(r))dr < sup |x| + / /3(r) dr < oo 
x&K J 



holds for all s £ J, thus {(,)^eS{{t}xK) is an equi-bounded family. 

By the Theorem of Arzela-Ascoli it follows that : G S({t} x K)} is pre-compact in 
C{J)"'. The closedness of : ^ G H({t} x K)} in C{J)"' was already shown in (ii), thus 
E{{t} xK)^ C{J)-^. 

If the right-hand side of ()3.1.ip is time- independent, it is obvious that we can put J = R. For 
any ^ G H(M x Q) with ^(f) = a; G we can define = C(- — r) for some r G M. Upon change 
of variable in the integral on the right-hand side of ()3.1.2I) it follows directly that ^ again 
satisfies the ordinary differential equation with modified initial condition ^(t + r) = x G 0, 
thus ^G H(M X 17). □ 

Theorem 3.1.4 (Unique solution sets). If there exists some a G Lloc^J) such that the right- 
hand side a satisfies the forward uniqueness condition 

(3.1.4) (x - y, a{s, x) - a{s, y)) < a{s)\x - 
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resp. the backward uniqueness condition 

(3.1.5) {x - y, a(s, x) - a{s, y)) > -a{s)\x - 

for almost all s ^ J and all x,y G M", then the solution set S(J x M") is forward (resp. 
backward) unique. 

Proof. We prove only the forward uniqueness (as the proof of backward uniqueness is analo- 
gous): Let ^1,^2 G H(Jx 0) with ^i(t) = ^2(*) = x. Setting ((s) := |Ci('S) — ^2(s)P we observe 
C,{t) = 1^1 (i) — S,2{t)\^ = and using that .^1,^2 satisfy the ordinary differential equation, we 
derive that the estimate 

dsCis) = 2{Us)-Us)Ai{s)-Us))=2{Us)-Us),a{s,ii{s))-a{s,i2{s))) 
holds for almost all s £ J, thus Gronwall's inequality gives 



Cis) < C{t) exp ^2 j a{r) dr j = 



for all s > t. It follows that ({s) = |^i(s) — ^2{s)\'^ = for all s > t, so S(J x 17) is forward 
unique. □ 

Theorem 3.1.5 (Forward unique flow map). Let H( J x M") be a forward unique solution set. 
Put := [t, oo[nJ, then the map defined by 

Xt ■■ Jt X 

(s,x) ^ i{s]t,x) 

for some ^(•; i, x) £ H( J x M") with ^(t; t,x) = x is uniquely defined, called forward flow map. 
It has the properties that 

(i) Xt is a continuous map, 

(a) Xt{s,-) is proper and onto for all s G fixed, and 
(Hi) Xr{s,Xt{r,x)) = Xt{s,x) for s > r > t. 

If the right-hand side of 13.1.1\) is time- independent, we have the additional flow property 
Xt{s,x) = Xois — t,x) for all s >t, so property (Hi) can be written in the well-known form 

Xh {s, Xt2 {r, x)) = Xti+t2 {s + r, x) 
for s £ J^,r £ Jf^ with ti, t2 £ J- 

Proof. First of all we note that xt is uniquely defined, as the solution set H( J x M") is forward 
unique. 
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(i) Considering the continuity of xt, we first prove that xt{s,x) is separately continuous. The 
continuity in s is obvious by the definition of xt- 

Continuity in x we shall prove by contradiction: Assume Xtis,x) is not continuous in x at 
fixed s. Then there exists a sequence {xk)k£N converging to a;, 7 > and ko £ N such that 
\Xtis,x) — Xt{s,xi:)\ > 7 for all k > k^. Choose K d M" such that x^ £ K for all /c € N. 

For each Xk we can find a G '^{{i} ^ such that ^fc(t) = Xk and xt{s,Xk) = Ckis) for all 
s E J^. Furthermore there exists a ^ G '^{{t} >^ K) such that ^{t) = x and xt{s,x) = ^{s) 
for all s G . The compactness of x K)) in C{J)^ yields that, there exists a locally 

uniformly converging subsequence ^ki i S "({t} x K)). The forward uniqueness of 

H(J X M") gives that ^(s) = ^(s) for all s G J^"*", thus lim^^oo \xt{s,x) — Xt{s,Xki)\ = < 7 
contradicting the assumption. 

The equi-continuity of H( J x M") yields the estimate 

\Xt{s, x) - xt{s', x')\ < \xt{s, x) - xt{s, x')\ + \xt{s, x') - xt{s' , x')\ < \xt{s, x) - xt{s, x')\ 

+ sup \C{s)-as% 

5eH(JxR") 

which, as xt was already shown to be separately continuous, implies the joint continuity of 
Xt- 

(ii) There exists some ^ G Et,x with xt{s,x) = ^(s), thus 



\Xtis,x) -x\ = 



a(r,^(r))dr 



t 



< / /3(r)dr 



holds for all s > t. For any compact set K d M" we have that Xt ^({■^l ^ — ■— 
{x £W \ 3xo e K : \ X — xq\ < I3{t) dr}, implying the properness of x 1-^ Xt(*) 

To show that x 1— >■ xt{s.,x) is onto for all s G J fixed, we have to show = W^. We 

prove by contradiction: Assume there exists some y xt{s,^^), then by Theorem 13.1.31 (i) 
there exists a solution ^ G 2^,^ with ,^(s) = y. Putting x := (,{t) we observe that G '^t,x and 
forward-uniqueness of H(J x M") yields that ^(r) = xtl?';^;) for all r > t. In particular we 
have y = ^(s) = xt(s,a;) G xt{s,W^) which is a contradiction. 

(iii) For fixed x G M" we have that xt{s,x) = ^i{s),s G J^"^ for some ^1 G H(J x M") with 
^i(t) = X. Fixing r G J^"^ we may choose .^2 £ '^{J x 1^") with ^2('^) = Xt{r,x) and ^2(5) = 
Xr{s,Xtir,x)) for s > r. Since ^1,^2 G S(J x M") with ^i(r) = Xil^^^a;) and ^2('^) = Xt{r,x) it 
follows by the forward uniqueness of H(J x M") that ^i(s) = ^2(5) for all s > r, thus 

Xr{s,Xt{r,x)) = Xt{s,x) 

for s G J+ n J+ and r G J+. 

In the case where the right-hand side of l3.1.1l is time- independent, we put J = M and have that 
is autonomous due to Theorem 13.1.21 As xo{s,x) = (,i{s) for some ^1 G H(M"+^) 
with ^i(O) = X and xt{s,x) = ^2{s) for some^^2 £ H(M""''-'^) with C2{t) = x. Due to the 
autonomy of the solution set it follows that (,2{s) := ^2('S — t) is again in H(M"+-'^) . As 
^2(0 = Ciii) = X the forward uniqueness of H(R"+^) immediately yields (,2{s) = £,i{s) for all 
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s G J^, thus xt{s,x) = xo{s — t,x) for all s G Jj^ . Using this relation and replacing r by 
r + t — t2 and s hy s + r + t — t2 — ti we can rewrite property (iii) in order to obtain 

Xh (s, Xt2 {r, x)) = Xh+t2 (s + r, x) 
for ah r G J^+, s G with ti,t2 G M. □ 

Theorem 3.1.6 (Backward unique flow maps). Let H( J x M") be a backward unique solution 
set. Put Jf :=] — oo,t] n J, then the map defined by 

{s,x) I— > x) 

for some G H(J x M") with S^{t;t,x) = x is uniquely defined, called backward flow 

map. 

It has the properties that 

(i) is a continuous map, 
(ii) <&t(s, •) is proper for all s G Jf fixed, and 
(iii) ^r{s, ^tir, x)) = ^tis, x) for s < r < t. 

If the right-hand side of i3.1.1\) is time-independent, we have the additional property x) = 
$o('5 — t,x) for all s <t, so property (iii) can be written in the well-known form 

'^t^ (s, (r, x)) = ^ti+t2 (s + r, x) 
for s G Jf^,r £ with ti, t2 G M. 



Proof. Analogous to the proof of Theorem I3.1.51 □ 

Remark 3.1.7 (Discontinuous flow maps). Due to Theorem 13.1.11 it is of course possible 
(regardless of any uniqueness property of the solution set H( J x W^)) to define a "global" flow 
map 6^ : J X ^ M"-, (s, x) ^(s) where ^ G H( J x M") with ^(t) = x. But such a map 
is not uniquely defined and 0f(s,-) is in general neither continuous, nor onto M"" for fixed 
s,t £ J. Nevertheless at least (Of (•, a;))^;^^™ is equi-continuous family in C{J)^ and 0f(s,-) 
is a proper map, i.e. for all K d there exists some K' d M" such that the inverse image 
0^^({s} X K) is contained in K'. In the case where H(J x M") is forward (resp. backward) 
unique, we have that 

\j+= Xt (resp. Gf \j-= ^t) 

holds. In the case when the solution set is unique, the globally defined flow map @t is uniquely 
defined by the forward flow on and the backward flow on Jf , as we will see. 

Theorem 3.1.8 (Unique flow maps). Let H(JxM") be a unique solution set, then it generates 
a unique, global ( defined on the whole interval J) flow map by 



@t{s,x) : = 

for t £ J fixed. It satisfies the properties 



Xt{s,x) s>t 
^t{s,x) s<t 
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(i) Qt is a continuous map, 

(ii) &tis, •) is proper for all s £ J fixed, and 

(Hi) 0s(r, @t{s, x)) = 0t(r, x) for s,r,t £ J with the special case (setting r = t) 
@s{t, @tis, •)) = idR" for s,t e J. 

We can write property ( Hi ) in terms of the forward and backward flow as 



for all s,t £ J with s >t. 

In the case where the right-hand side of /13.1.1\) is time independent we again obtain Qt{s, x) = 
@o{s — t, x) for all s,t £ J and x £ W^. 

Proof. The uniqueness of Qt follows immediately by combining the uniqueness results from 
Theorem 13.1.51 and Theorem I3.1.61 It is globally defined since Jf |J = J. 

(i) The continuity of Qt is obvious by Definition, as Xtit, •) = ^tit, •) = idjRn, xt is continuous 
on Jf^ X R" and is continuous on Jf x M" due to Theorem 13.1.51 and 13 . 1 . 6l 

(ii) Qt{s,-) is proper follows immediately from the properness of Xtis,-) ^tis,) on the 
respective domains as proved in Theorem 13. 1.51 and 13.1.61 

(iii) The case where r<s<toit<r<s follows again directly from Theorem 13.1.51 
and 13.1.61 The r<s<t, t<r<s, r<t<s and s < t < r need some additional 
work, but the proof is analogous to the respective proofs in Theorem 13.1.51 and 13.1.61 To 
illustrate the proving technique one more time, we consider the case r < s < t: We have 
Qs{r,Qt{s,x) = ^s{r,Xt{s,x)) = (,i{r) for some £ H(J x W^) with ^i(s) = xtis,x) and 
0((r, = xt(r, x) = ^2{t) for some £ H(J x M") with ^2{t) = x. If r = s we have 
^i(s) = ^2{s), thus uniqueness of H(J x M"") implies that = ^2 on the whole interval J. 

The relation Qt{s, x) = Qo{s — t, x) for all s,t £ J and x £ M" follows directly from Theorem 



Remark 3.1.9 (C^ regularity of the flow). Let a satisfy the Caratheodory condition on 
J X M" and assume that a{s, •) £ C^(]R")'^ for almost ah s £ J. Then a the solution set H(J) 
is unique and the uniquely generated flow map Qt satisfies the properties 



^s{t,-) = {Xt{s,-)) 



1 



[3T5] and [3T6l 



□ 




and 




for all (s,x) G J X M". For a proof we refer to [25t Theorem 1.2.5.]. 
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3.2 Differential inclusions and integral inequalities 

In this section we introduce a natural generalization of the ordinary differential equation, 
namely the differential inclusion which is of the form 



where At^x is a non-empty, closed and convex subset of M" for all {t,x) €z J x M". 

It turns out that the notion of a set-valued map (cf. Section [TT2]) is suitable for the description 
of such generalized right-hand sides of the form {t,x) i— > At^x- 

As in the case of Caratheodory theory we have to impose certain conditions on the set-valued 
map A in order to obtain an existence result for the differential inclusions. 

Throughout this section we assume that the set-valued map A satisfies 

(i) t ^ Af^x is Lebesgue measurable on J for all fixed x G M", 

(ii) X ^ At^x is a upper semi-continuous for almost all t £ J, 

(iii) there exists a positive function /3 G L\^^{J) such that sup^jg^n \At^x\ < I3{t) almost 
everywhere on J, 

referred to as Filippov conditions (FC). 

It should be noted that sup^.g]gn \At^x\ < (3{t) implies that x ^ At^x is bounded for almost all 
t G J. Due to Theorem 11.2.141 condition (ii) is equivalent to have x At^x being a closed 
set-valued map for almost all t £ J. 

We should think of t i— > At^x as a set-valued map in Ll^^{J;ICo{W^)). This means that the 
Filippov conditions do not depend on the chosen representative of t ^ At^x for x £ M"" fixed, 
so without loss of generality we can assume that x i — > At x 

is upper semi-continuous (or 

equivalently - closed) for all t £ J. 

Equivalently we can formulate these conditions for the supporting H (according to Theorem 



(i) 1 1-^ H{t, X, w) is Lebesgue measurable on J for all x,w £ M", 

(ii) X ^ H{t,x,w) is upper semi-continuous for almost all t G J and w £ M", 

(iii) there exists a positive function /3 G L\oc{J) such that sup^g]gn \H{t, x,w/\w\)\ < (3{t) 
for almost alH G J and all w £ M". 

Furthermore the differential inclusion p.2.6p is equivalent to the following integral inequality 



(3.2.6) 



?(s) G As^^i^s) a.e., 



m 



X 



I1.2.14P by 



s 



(3.2.7) 




t 
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Theorem 3.2.1 (Existence theorem for differential inclusions). Let A be set-valued map 
satisfying (FC). Let {t,x) €z J x M" be the initial value. Then there exists a solution ^ G 
^C(J)" of the differential inclusion 



X. 



Proof. First of all we notice that for any / G C{J)^ it holds that the composition t ^tj(i) 
is a Lebesgue measurable set-valued map from J to M". By Definition 11.2.251 and Corollary 
11.2.211 the measurability of the set-valued map t i-^ ■^tj{t) is equivalent to the measurability 
its supporting function t ^ H{t, f{t),w) for all w G M", which follows by Lemma lA.91 since 
H satisfies the Filippov condition (i) and (ii). 

We construct a sequence of continuous functions {S,k)k£N which approximates the solution of 
the differential inclusion. 

Put S,o{s) = X for all s £ .J, then using Theorem 1 1 . 2 . 28 1 we iteratively construct a sequence of 
set- valued maps {Xk)kGn with : J /Co(IR") by 

s 

Xk+i{s) = {x} + y"^,,^,(,) dr 

t 

Furthermore Theorem 11.2.281 yields an absolutely continuous selection S,k+i of Xf^^i with 
^k+i{t) = X, so we obtain a family {(,k)k<^N of absolutely continuous functions with the property 

s 

(3.2.8) ms)-Ckir),w) < J H{T,Ck{T),w)dT 

r 

for s,r £ J and w G M". 

We observe that the (FC) (iii) yields the estimate 

max{r,s) max (r,s) 

Ms)-^kir)\< [ sup H{T,^kiT),w/\w\)dT< [ /3(r)dT for s,rG J. 
J u>eK" J 

min (r,s) min(r,s) 

Thus (Cfc)feeN is an equi-continuous family. The estimate 



s 

\^k{s)\ <\x\+ sup H{T,^kir),w/\w\)dT 



r 



yields the equi-boundedness of {Ck)keN- 

By the Theorem of Arzela Ascoli there exists an uniformly convergent (on compact subsets 
of J) subsequence ^fc. ■'— ?° where ^ G (7(7)" with ^(t) = x. 



"■J 

We have for all s, r G J, 



(3.2.9) \C{s) - ^(r)| < |^(.) - 6,(^)1 + 16, (5) - 6,(01 + |6,(r) - C(r)| 

max{r,s) max(r,s) 



min(r,s) min (r,s) 
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hence is absolutely continuous in J. 

By the upper semi-continuity oi x ^ H {t, x , w) , it follows that 

limsup H{t,Ck,{T),w) < H(t,C(t),w). 

The theorem of dominated convergence [21. Theorem 12.24] yields 

s s 

{^{s),w) <limsup{£,kj {s),w) = {x,w)+lim sup / H{T,£,kj{T),w) dr < {x,w)+ H{t,^{t),w) dr 

jr— »00 J— >00 J J 

t t 

Differentiation with respect to s yields 

{as),w) <H{s,as),w) 
for all w G M", thus S,{s) G ^s,^(s) foi' almost all s G J and ^(s) = x. □ 

Definition 3.2.2. Let he a non-empty subset of M" . Analogous to the case of classical 
ordinary differential equations, we define the inclusion solution set hy 

Ej{J X n) := {i G AC{JY I 3(t,x) G J X 17 : i{s) G a.e., ^{t) = x]. 

If we put r.t^x '■= '^liii't, x)}) it obviously holds that 

The solution set E^i^J x il) is called forward (resp. backward) unique, if S,i,S,2 £ 'Ej{J x J]) 
with Ci{t) = C2{t) implies ^i(s) = ^2(s) for all s > t (resp. s <t). 

In the case where J = M., we refer to a solution set H/(IR x il.) as autonomous, i/^ G x Q) 

implies + r) G H/(M x Q) for any r G M. 

Proposition 3.2.3 (Properties of the inclusion solution set). Let Q be a non-empty subset 
o/M", then the inclusion solution set S/(J x Q) has the following properties: 

(i) For every {t,x) G J x il, there exists G H/(J x Q) with ^{t) = x. 
(a) If t G J and Q is closed, then x Vl) is a closed subset of C{J)"'. 

(Hi) (Ogg=j{jxn) equi- continuous family in C{J)^. 

(iv) Ifn<£ M", then x ft) is a compact subset of C{J)'^. 

If the right-hand side of !i'J. 2. 6]) is time-independent, we can put J = M and the (Filippov) 
solution set H/(M x i}) is autonomous. 

Proof, (i) By Theorem 13.2.11 the set H/(J x il) is a non-empty subset of C{J)^, such that for 
every {t,x) G Q there exists a ^ G H/(J x fl) with ^(t) = x. 
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(ii) First we prove that x fi) is a closed set in C{J)^ if Q is closed: Let {^k)keN be 
sequence in x Q) uniformly converging (on compact sets) to some ^ G C{J)^. In 
particular xt = ik{t) ^{t) ='■ x ^Vt, since Vt is closed. It follows that 

{i{s),w) = lim = lim {xk,w) + / i?(T, ^fc(r), dr 

s 

<{x,w)+ / limsup //(r, ^fc(r), tt;) dr < (x, + / H{t,^{t),w) dr 

J k^oo J 
t t 

which implies ^ G x ri) satisfying the initial condition ^[t) = x G fi, so x il) is 

a closed subset of (7(7)"". 

(iii) For all t £ J, there exist some compact neighborhood Kf d J, such that by (j3.2.9p 



(3.2.10) |e(s)-^(r)|< 1 /3(r) 



max (r,s) 
mill (r,s) 



holds uniformly for all ^ G S/(J x il), thus ^/(J x (7) is a equi-continuous family, 
(iv) Since K d M" the estimate 



sup |C(s)| < sup |x| + sup 

ieEi{{t}xK) xeK ieSi({t}xK) 



sup H{t, ^,{t),'w/\w\) dr 



t 



max 



< sup |x| + / /?(t) dr < cxd 



min 

holds for all s G J, thus (0^eH/{{t}x_ft') is an equi-bounded family. 

By the Theorem of Arzela-Ascoli it follows that : ^ G H/({t} x K)} is pre-compact in 
C{J)"'. The closedness of : ^ G S/({t} x K)} in C(J)" was already shown in (ii), thus 
Ej{{t} xK)^ C{Jr. 

If the rightside of (j3.2.6|) is time- independent, it is obvious that we can put J = M. For 
any ^ G H/(M x f]) with ^(t) = x £ Q we can define = r) for some r G M. By 

integral substitution on the right-hand side of (|3.1.2p it follows directly that ^ again satisfies 
the ordinary differential equation with modified initial condition ^(t + r) = x G il, thus 
^GS7(Mxfi). □ 

Theorem 3.2.4 (Properties of the inclusion solution set). Let A satisfy conditions (FC) and 
assume there exists some positive a{t) G L\^^{J) such that the supporting function H of A 
satisfies 

(3.2.11) iJ(s,x,x - y) + H{s,y,y - x) < a{s)\x - 
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resp. 

(3.2.12) H{s, x,x-y) + H{s, y,y-x)> -a{s)\x - y\^) 

for all x,y € M" and almost all s G J, then the solution set 'Ej(J x M") is a forward (resp. a 
backward) unique transport set. 

Proof. We prove only the forward uniqueness (the proof of backward uniqueness is done 
analogously): Let 6, 6 e S/(J x M") with ^i(i) = ^(i) = x. Setting C(s) := |6(s) - 6(s)P 
we observe C{t) = \Ci{t) — '^2(i)P = and using the definition of Ei{J x M") and the forward 
uniqueness condition ()3.2.11|) . we obtain the estimate 

dsCis) = 2{^i{s)-Us),ii{s)-Us))<2{H{s,C,{s),^,{s)-Us))+H{s,Us),Us)-U 
< 2a(s)C(s) 

for almost all s G J, thus Gronwall's inequality gives 

C(s) < C(i)exp ^2 j a(r)drj = 0. 

It follows that C(s) = |6(s) - 6(s)P = for all s > t if ^i(t) = 6(t), so Ei{J x R") is a 
forward unique transport set. □ 

Theorem 3.2.5 (Forward unique (inclusion) flow map). Let 'Bj{J x M") be a forward unique 
solution set. Put := [t,oo[nJ, then the map defined by 

Xt ■■ Jt X 

{s,x) ^ i{s]t,x) 

for some £,{-;t,x) G S/(J x M") with S^{t;t,x) = x is uniquely defined, called (inclusion) 
forward flow map. 

It satisfies the following properties: 
(i) Xt is a continuous map, 

(a) Xt{s,-) is proper and onto for all s G fixed, and 
(Hi) Xr{s,Xt{r,x)) = Xt{s,x) for s>r>t. 

If the right-hand side of 13.2.6\) is time-independent, we have the additional flow property 
Xt{s,x) = Xois — t,x) for all s >t, so property (Hi) can be written in the well-known form 

Xh {s, Xt2 {r, x)) = Xti+t2 {s + r, x) 
for s G J^,r G J^^ with ti,t2 G M. 
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Proof. First of all we note that xt is uniquely defined, as the solution set H/( J x M") is forward 
unique. 

(i) Considering the continuity of xt, we first prove that xtis,x) is separately continuous. The 
continuity in s is obvious by the definition of xt- 

Continuity in x we shall prove by contradiction: Assume xt{s,x) is not continuous in x at 
fixed s. Then there exists a sequence {xk)km converging to a;, 7 > and ko G N such that 
\Xt{s, x) — Xt{s, xi^)\ > 7 for all k > k^. Choose d M" such that x^ £ K for all k £ N. 

For each x^ we can find a E Ef{{t} x K) such that ^fc(i) = Xk and xtis,Xk) = Cfc('S) for all 
s G . Furthermore there exists a ^ S Ej{{t} x K) such that ^(t) = x and xt{s,x) = (,{s) 
for all s E . The compactness of x K)) in C{.J)^ yields that, there exists a locally 

uniformly converging subsequence ^ G ^/({^} ^))- The forward uniqueness of 

H/(J X M") gives that ^(s) = ^(s) for all s G J^^, thus lim^^oo \xt{s,x) — Xt{s,Xki)\ = < 7 
contradicting the assumption. 

The equi-continuity of H/(J x M") yields the estimate 

\Xt{s, x) - xt{s', x')\ < \xt{s, x) - xt{s, x')\ + \xt{s, x') - xt{s' , x')\ < \xt{s, x) - xt{s, x')\ 

+ sup \^{s) - C{s% 

56Hj{JxR") 

which, as xt was already shown to be separately continuous, implies the joint continuity of 
Xt- 

(ii) There exists some ^ G Et^x with xtis,x) = ^(s), thus 



\Xt{s, x) -x\= sup {xt{s, x) -x,w) < 



t 



s 

sup {H{T,C{T),w/\w\)dT < //3(r)dT 



holds for all s > t. For any compact set K <s we have that Xt ^ ^) — ■— 

{x G I 3xo £ K : \x — xo\ < (3{t) dr}, implying the properness of x 1-^ Xt{s, x). 

To show that x xt{s,x) is onto for all s G J fixed, we have to show xtis,^"') = IR"- We 
prove by contradiction: Assume there exists some y xt{s,M.'^), then by Theorem 13.2.31 (i) 
there exists a solution ^ G Eg^y with ^(s) = y. Putting x := ^{t) we observe that ^ G H^^a, and 
forward-uniqueness of H/(J x M") yields that ^(r) = Xf(r, x) for all r > t. In particular we 
have y = ^(s) = xt{s,x) G XiCS)!^"") which is a contradiction. 

(in) For fixed x G M" we have that xt{s,x) = £,i{s),s G Jf^ for some ^1 G Ej{J x M"-) 
with ^i(t) = X. Fixing r G J^"^ we may choose ^2 £ ^/l^^ x I^") with C2{f) = Xt{f,x) 
and ^2('S) = XriSjXti'f'jX)) for s > r. Since ^1,^2 G x with ^i(r) = xt{r,x) and 

C2('^) = Xil'";^) it follows by the forward uniqueness of H/(J x M") that ^i(s) = ^2('5) for all 
s > r, thus 

Xr(s,Xt(?",a;)) = xt{s,x) 

for s G J+ n J+ and r £ . 

In the case where the right-hand side of l3.2.6l is time- independent, we put J = M and have that 

is autonomous due to Theorem 13.1.21 As xois,x) = for some ^1 G 
with Ci(0) = X and xt(s,x) = ^2(5) for some ^2 £ with £,2{t) = x. Due to the 



43 



3 Generalized solutions of ordinary differential equations 



autonomy of the solution set it follows that ^2{s) := ^2{s — t) is again in H/(M"^^) . As 
^2(0 = = X the forward uniqueness of immediately yields i2{s) = Ci(-s) for all 

s G , thus Xt{s,x) = Xq{s — t,x) for all s G . Using this relation and replacing r by 
r + t — t2 and s hy s + r + t — t2 — ti we can rewrite property (iii) in order to obtain 

Xh («> Xt2 {r, x)) = Xh+t2 {s + r, x) 

for all r G s e with ti, t2 G IK. □ 

Theorem 3.2.6 (Backward unique flow maps). Let H/(JxM") be a backward unique solution 
set. Put :=] — 00, t] n J, then the map defined by 

{s,x) ^ i{s]t,x) 

for some ^{■■,t,x) £ '^i{J x M") with S^{t;t,x) = x is uniquely defined, called (inclusion) 
backward flow map. 

It has the properties that 

(i) is a continuous map, 
(a) ^tis, •) is proper for all s G fixed, and 
(iii) ^r{s, ^t{r, x)) = ^t{s, x) for s < r < t. 

If the right-hand side of \3.2.6\ is time-independent, we have the additional property ^tis,x) = 
$o(s — t^x) for all s <t, so property (iii) can be written in the well-known form 

$j,(s,$t2(r,x)) = $t,+t2(s + r,x) 

for s e Jj^,r £ with ti, t2 G ^■ 

Proof. Analogous to the proof of Theorem 13.2.51 □ 

3.3 Generalized solution concepts for ordinary differential equations 

It is possible to generalize the classical solution concept for ordinary differential equations 
with discontinuous right-hand side a G Lj^^{J; L°^{W^))^, by replacing a with a suitable set- 
valued map A satisfying (FC) and solving the resulting differential inclusion by means of 
Theorem 13.2.11 

Generalization of the classical solution concept means that whenever a right-hand side a 
satisfies the Caratheodory conditions the graph of the set-valued map A and the graph of 
a G ^[^^(J; L°°(M"))" are equal. One way of mapping functions a G Ll^^{J; L°°(K'"))"- into 
set-valued maps satisfying (FC) is by means of the essential convex hull. 

Another method to generalize the classical solution concept, is to regularize the right-hand 
side a to obtain a net of smooth functions (^e)£g]o.i] with lim£_^o(^e) V') = foi' 
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if G 'D(W"'). Then a net of solutions (^,s)e is obtained by solving the family of ordinary 
differential equations 

(3.3.13) Us) = Ae{s, Us)), Ute) = Xe, 

where (i£,Xe)e is converging net with liuisr^Q^t^^Xe) = {t,x). 

Such nets of solutions {Ce)e have convergent subnets (^r(e))e with t £ T and ^^(g) ^ G 
C{ J)^ uniformly on compact subsets of J, when the net of smooth functions {A^)^ satisfies 
certain boundedness assumptions. The generalized solution ^ does in general not solve any 
differential equation, but it will turn out that ^ solves a differential inclusion. So obviously 
these two solution concepts are closely related. 

3.3.1 The Filippov solution concept 

If a G L\^^{J\L°°{W^))'^ we can consider its essential convex hull At^x for almost all t G J 
and all x G M". We show that A : (t, x) i— > At^x determines a set-valued map, satisfying the 
conditions (FC). Furthermore the essential convex hull At^x is equal to {a{t,x)} for almost 
all t G J and all x G M", if a satisfies the Caratheodory conditions. 

By Theorem 13.2.11 there exists a solution ^ of differential inclusion 

i{s) G -4s,^(s) a.e., ^(t) = x, 

which we call Filippov (generalized) solution of the ordinary differential equation ^ = a{s,^{s)) 
when (t, x) i-^ At^x denotes the essential convex hull of a. 

This type of generalized solution concept for ordinary differential equations with discontinuous 
right-hand side was presented in [H], [2t Chapter 2], and [25i Chapter 1.4]. 

Definition 3.3.1. Let a G -Lj'^^( J; and define the essential supporting function for 

a by 

Ha{t,x,w) := liinesssupygB^(^) {a{t,y),w). 

which is defined for almost all t £ J and all x,w £ M". 
The generated set-valued map A : (t, x) i— > At^x obtained by 

At,x ■■= {y G I {y,w) < Ha{t,x,w),w G M"}. 
is called essential convex hull of a. 

Proposition 3.3.2. The essential convex hull of any function a G Lj^^{J; (M."'))^ satisfies 
(FC). 

Proof. By Definition 13.3.11 it is obvious that w i— > Ha{t, x,w) is positively homogeneous and 
convex in w for almost all t G J and all x G M". 

We put 

h5{t,x) := esssupy^Bs{x){a{i,y),w) 
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By |2H Theorem 20.14] there exists a null set A^, such that 

hs{t,x) = sup {a{t,y),w). 

yCiBsix)/N 

Due to Lemma lA.51 we can find a sequence {xk)ken in Bs{x)/N such that 

hs{t,x,w) =lim sup {a{t,Xk),w), 

thus {t, x) I— > hs{t, X, w) is measurable as the pointwise limes superior of a family of measurable 
functions. 

Observe that hs^(t,x,w) < hs^it^XjW) for 62 > Si, so {hs{t,x,w))s(z]o^i] is a decreasing net, 
bounded for almost all t £ J. For y £ B^{x)/N it holds that hs{t,y,w) > {a{t,x),w), thus 
hs{t,y,w) is bounded from below for almost all t £ J. It follows that 

lim hg{t,x,w) = inf hg{t,x,w) 
S^o <5e]o,i] 

and since the pointwise infimum of a net of upper semi-continuous functions is again upper 
semi-continuous (see [211 Exercise 7.23 (a)]), we have that x Ha{t,x,w) is upper semi- 
continuous for almost all t G J and all w £ M". 

We have 

Hait,x,w/\w\) < hi{t,x,w/\w\) < esssupygiRn \a{t,y)\ £ ll^^iJ) 
since a £ L\^^{J] L°°(M"))". □ 

Theorem 3.3.3 (Approximation property of the essential convex hull). Let a £ L\^^{J] L°°{W^)Y 
and let p he a positive moUifier in 2?(IR") with J p{z) dz = 1. Put Pe(x) = p{x / e) . 

Then the approximating net as{s,x) = a{s,-) * pe has the property, that if (^e)e is a net of 
smooth functions converging to some ^ £ C{ jy^ uniformly on compact sets, then it holds that 

limsup(ae(s,Ce(s)),u;) < Ha{s,C{s),w) 

for almost all s £ J. 

Proof. Let i? > be such that suppp C Br{0). Then we observe that 

sup (a^(s,^^(s)),u;) = sup {a{s,y),w)p^{(,^{s)-y) dy < sup sup {a{s,y), w) < 00. 

for almost all s G J. 
Let K ^ J fixed. 

Then for ah e > there exists p'{e) £]0,e/{2R)] such that |^^(s) - ^(s)| < e/2 for all s G K 
and p £]0,p'{e)]. It follows that 

BR.,iU^)) C BR.,{as)) + B,/2{as)) = i?^.^+e/2(e(s)) C BMs)) 
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for all s G K and n G]0,/i'(e)], implying 



sup sup {a{s,y),'w) < sup {a{s,y),'w). 



Observe that |^sup^g]o,£] j ^ decreasing, bounded net for almost all s S J, 
thus converging. We conclude that any subnet of the net {sup^^^Q^^^{afj,{s,(,^{s)),w))£, has 
the same limit, thus 

lini sup (0^(5,^^(5)),?!;) = lini sup {a^{s,^^{s)),w). 

Finally we obtain 

limsup(ae(s,^e(s)),u;) = lim sup (o^(s, ^^(s)), -w) = lim sup (a^(s, ^^(s)), u;) 

< lini sup {a{s,y),w) = Ha{s,^{s),w) 

for almost all s € K and all w G M"-. □ 

Remark 3.3.4. Theorem 13.3.31 shows that a Filippov generalized solution for an ordinary 
differential equation with discontinuous right-hand side a G L\oc{Jj can be obtained 

by regularizing the right-hand side by convolution with a positive, smooth and compactly 
supported mollifier (instead of applying the existence theorem 13.2.11 for differential inclusions). 
Such an approach for constructing Filippov generalized solutions is pursued in [251 Chapter 
1.4]. 

Putting as{s,x) = a{s, •) * ps we have that any net of regular functions {(,e)e solving 
has the property that 

s 

limsup(^£(s),'w) = {x,w) + limsup / {ae{T, ^s{r)),w) dr 

e^O e^O J 

t 

s s 

< {x,w) + / limsup(ae(T,^£(r)),'w) dr < {x,w) + / i7a(T, lim sup ^e(r), dr 
J e^O J e^O 

t t 

for all w G R", which implies that any subnet {$,T(e))e (with t £ T) converging uniformly on 
compact subsets to some ^ G C{J)'^, satisfies the differential inclusion 

^(s) G ^{t)=x, 

where A denotes the essential convex hull. 

Theorem 3.3.5. Let a G Li,^( J; L°°(M"))" satisfy the forward uniqueness condition 
(3.3.14) {x - y, a(s, x) - a{s, y)) < a{s)\x - 
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resp. the backward uniqueness condition 

(3.3.15) {x - y, a(s, x) - a{s, y)) > -a{s)\x - y\^, 

for some positive a G L^^^{J), then the supporting function of the essential convex hull Hj sat- 
isfies the forward uniqueness condition \3.2.11\) resp. backward uniqueness condition 13.2. 12\) . 

Proof. We prove only the case where a satisfies the forward- uniqueness condition. 

{a{s, zi),x-y)-{a{s, Z2),x-y) = {a{s, zi), zi-Z2)+{a{s, Z2), Z2-zi)+2l3{s){\zi-x\ + \z2-y\) 

< \zi - Z2\'^a{s) + 46f3{s) 

for all zi G Bs{x), Z2 S Bs(y) and almost all s £ J. It follows that 

esssup^g5_5(^.)(a(s,2:), j;-y)+esssup^gB4(?;)(a('5>^)>y-^) < sup \zi-Z2\'^a{s)+46/3{s) 

(zi,Z2)&Bs{x)xBs{y) 

= \x- y\^a{s) + 26a{s) + 4(5/3(s), 

thus 

Ha{s, x, x-y)+Ha{s, y, y-x) = |im(ess sup2eB^(^.) (a(s, z), x-y)+ess sup^gB^{j/)(a(s, z),y-x)) 

< |x — y\^a{s) 

holds for almost all s G J. □ 
3.3.2 Colombeau solutions for ordinary differential equations 

Finally we consider generalized solutions for ordinary differential equation in the setting of 
Colombeau generalized functions. 

Let (t, x) G J X M" be the initial condition, then we study the ordinary differential equation 

where A := [(Ae)^] G Q[J x M",R"] such that there exists some representative satisfying the 
bound 

(3.3.16) sup |^e(t,x)| < (3{t), e g]0, 1], almost everywhere in t G J 

for some positive function (] G L^^^{J). It will turn out that the generalized solution ^ := 
[(Ce)e] then allows for picking a subnet (CT(e))e converging uniformly on compact subsets of J 
to some C G ^C( J)" with C{t) = x. 

Theorem 3.3.6 (Existence). Let A G G[J x M^jM"] be Colombeau map of natural type with 
a representative (^e)£ satisfying 113. 3.1b]) . 

If{t,x) G J X M", then there exists a c-bounded solution ^ G Q{J)^ to the initial value problem 

i{s) = A{s,i{s)), m = x. 
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Furthermore, there exists some {t,x) £ J x Q, C E AC (J)" such that for any representa- 
tive {Ce)e of ^, (t£,X£)e of (t,x) there exists subnets (tr(£)) '^T{e))e,{^T{e))e with T and 

lim£_^o(^r(e)) 2?T(e)) = {t-,x) and ^^-(^i;^ uniformly on compact sets of J. 

It satisfies the integral inequality 

s 

(3.3.17) {Cis),w) < {x, w)+ J Ha{t, ar),w) dr 

t 

for all w G M", where {x,w) HA^t, x,w) denotes the supporting function of the generalized 
graph Giaph{A{t, •)). 

Proof. Let (^e)e be a representative of A and let (tejXe)^ be a representative of {t,x). By- 
classical existence and uniqueness we obtain a solution for each e €]0, 1] such that 



s 

Us) = + j ^e(r,e.(T))dr 



holds. ()3.3.16p yields |Ce('S)| < l^el + j3{T)dT for all s G J, hence is bounded on 

compact subsets of J. 

Since (tejXj)^ is bounded we have that for all s G J there exists some C > such that 
|Ce('5)| < C for all £ g]0, 1], thus {S,e)e is an equi-bounded family. The estimate 

max (r,s) 

\Us)-Ur)\< I P{r)dT s,reJ,ee]0,l]. 

min (r,s) 

yields that the family {S,e)e is equi-continuous. 

The theorem of Arzela-Ascoli yields a subnet (CT(e))e with t £ T converging to some C G C{J)^ 
uniformly on compact sets. Without loss of generality we can assume that (ir(£)) 3;T-(e))e 
converges to some {t,x) G J x M" with ^(t) = x. 

We have for all s,r £ J, 

m - C(r)| < \C{S) - + \^rie){s) - Ue)i^)\ + \^rie){r) " C(r)| 

r r 

<IC(s)-e.(e)(^)l + / /3(T)dr + |C,(,)(r)-C(r)r^°y" /?(t) dr, 

s s 

hence C is absolutely continuous on J. 

Note that the equi-continuity of (^e)e implies its c-boundedness. Then the moderateness 
of (^e)e follows by the same argument as in the proof of [18j [Proposition 1.2.8], since the 
derivative ^^(s) may be written as the composition of the moderate net A^ and the c-bounded 
net {s,Us))- 
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By Theorem 12.2.41 we obtain that A induces a supporting function satisfying (FC). By the 
theorem of dominated convergence [211 Theorem 12.24] and Theorem 12.2.51 it follows that 



{Cis),w) = lim {Crie),w) < {X,W) +limsup / ^r(e)(T, ^r(e) (t)) dr 



t 



5 

< {x,w) + / limsup^^(e)(T,^^(£)(r))dr 



t 

s 

< {x,w) + J HA{r,C{T),w)dT 
t 

for all t(; G M". □ 
Remark 3.3.7. By Proposition 12. 1.6) the equi-continuity of {S,s)e yields that 

Graph(0= U{s}xCP((C.(s))e). 

Then ( £ AC (J)"' as in ()3.3.17p is a selection of the set- valued map s i-^ Graph(^)s, i.e. 
C(s) e Graph(^)s for all s e J. 

If (r, z) E Graph(^) we can find some t G T such that lim^^o '^T{e) (^) = 

Then the Theorem of Arzela-Ascoli yields that it is possible to choose a subnet of (Cr(e))e 
converging uniformly on compact sets to some C £ AC{J)^ with C(r) = z and ({t) = x. 

This shows that the Graph(^) is the union of graph(^) of all possible C G AC{J)^ as in 
Theorem 13.3.61 We can reformulate property (j3.3.17p to obtain 

s 

(3.3.18) Graph(Os {x} + j ch(Graph(^) o ({r} x Graph(Or))dr 

t 

where "o" denotes the composition of set- valued maps (cf. Theorem II. 2. 19p and the integral 
is according to Definition 11.2.251 

Example 3.3.8 (Hurd-Sattinger characteristic flow). Consider the ordinary differential equa- 
tion: 

^(s) = 0(-e(^)), m = x^ 

where 9 = [(6e)e] G Q {M) is defined by 0£(x) = Pe * H{—-) with = \p{-/e) for some 
p e 5(M) with / p{x) dx = 1. 



We put g{x) := p{z) dz and denote Ran{g) = [a_,a+] for some a_,a+ G M. Note that 
liuix^-oo g{x) = 1 and linix^oo g{x) = implies q_ < and a-f > 1. It holds that a_ = 
and a+ = 1 if /3 > 0. 

In our example we allow p to have non-vanishing moments, so Q is not necessarily the em- 
bedding of the Heaviside function. 

Nevertheless {Qe)e converges to the Heaviside function in Lj'Q^(M) and since sup3..gR x)| < 
sup^g]g |(7(a;/e)| < q+ it satisfies (I3.3.16p . 
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Due to Theorem 13.3.61 there exists a c-bounded Colombeau solution ^ = [(Ce)^]) such that 
there exists a converging subnet with hm^^o ^T(e) ~^ C uniformly on compact set for some 

As in Example 1 2 . 1 . 4 1 we obtain that the generalized graph of Q is 



Graph(G(--)) := ({x G M | x < 0} x {!}) U (^{0} x Ran(c/)j U ({x G M | x > 0} x {0}) . 
Using its support function (Example I1.2.8P the property (|3.3.17p of C translates to 

s s 

x + j lc(r)<o(7")dr < C(s) < x + j l(^(^)<o(r)dT. 



Then it is straight-forward to derive that 
(3.3.19) C(s;i,x) G < 



{min (x + s — t, 0)} x<0 

{0} x = 0,s>t 

[s-t,0] x = 0,s<t 

{x} X > 0. 



holds, which implies that Q is forward unique. It follows that 



lim sup |^,(s)-C(s)| =0 

s£K^^[t,+oo[ 



for all (£ M. 



Note that the support function of Graph(0(— •)) does not satisfy the forward uniqueness 
condition ()3.2.1ip if a_|_ > 1 or a_ < 0, nevertheless C, is forward unique due to ()3.3.19p 
regardless of the values a_, a+. 

By Remark 13.3.71 we obtain that the Colombeau solution ^ of the Hurd-Sattinger equation 
has the property that 



Graph(Os ^ {x} + j ch(Graph(e) o ({r} x Graph(e)r)) dr = {x} + j I 

t 

which in (13.3. 19p implies 



Graph(Or<o('r) dr 



Graph(e). C <^ 



{min (x + s — t, 0)} x<0 

{0} x = 0,s>t 

[s-t,0] x = 0,s<t 

{x} X > 0. 



Theorem 3.3.9 (Forward uniqueness of the distributional shadow). Let A G ^(fi-jn)" with 
representative (Ae)^ satisfying i3.3.16\) and in addition the forward uniqueness condition 

(3.3.20) {As{s,x) - A,,{s,y),x- y) < a{s)\x - y\^ + 7s,e'{s) 
resp. the backward uniqueness condition 

(3.3.21) {A,{s, x) - Ae'is, y),x-y)> -a{s)\x - - je,e'{s) 
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3 Generalized solutions of ordinary differential equations 

holds for some positive functions a,7£^e' G L^{[0,T]) such that Jq ■Js,e'{'T) dr = 0{e + e') as 
e,e'^0. 

Then the supporting function Ha of the generalized graph Graph(A(s, •)) satisfies the forward 
uniqueness property 

Ha{s,x,x- y) + HA{s,y,y- x) < a{s)\x - 

resp. 

Ha{s,x,x -y) + HA{s,y,y - x) > -a{s)\x - yp). 

Proof. By the proof of Theorem 12.2.21 we can write 

Ha{s,x,x- y) = lim(^^^(e)(a;e),a; - y) 

and 

HA{s,y,y- x) = lim{A^^(^^){ye),y - x) 

for fixed x,y £ and for some ti,T2 £ T and nets {xs)e, {ye)e with x^ x, y^ ^ y. For 
almost all s G J and all x,y G M" we have 

Ha{s,x,x- y) +HA{s,y,y- x) = \\M{A^^(^^){s,Xe),x - y) + {A^^(^^){s,ye),y - x)) 

e^O 

- }^Jy'\^ri{e){s,Xe),Xe - ye) - (^^3 (e) ^e) ' " Ve) + P{s){\x - Xe\ + \y - ye\)) 

< lini(a(s)|Xe - yef + (3{s){\x - X^] + \y - ye\) + 7ri(e).T2 (£)(*)) < Oi{t)\x - y\^. 

The proof for the backward uniqueness condition is analogous. □ 

Remark 3.3.10. By means of Theorems 13.3.91 it immediately follows that the 'sub-shadow' 
( of the Colombeau solution ^ = [(^£)£] € G (J)^ (as obtained by Theorem 13.3.60 . is for- 
ward unique resp. backward unique if the right-hand side A := [(Ae)e] S G[J x M",R"] has 
the property that the supporting function of Graph(A(i, •)) satisfies the forward uniqueness 
condition 13.3.201 resp. the backward uniqueness condition 13.3.211 

Example 3.3.11. Let a £ L^{R,L'^{W)f satisfy the classical forward uniqueness condition 
(3.3.22) {a{t,x) -a{t,y),x -y) <a{t)\x-y\^ for almost all {t,x),{t,y) e M"+^ 
for some positive 5 G L-'^(]R). 

Assume p £ 2?(M"+^) positive, f^np{z)dz = 1 and define pe^tjx) := e"'"-^^ p{t/e, x/e). Con- 
sider the Colombeau function A := [{Ae)e] £ e[M"+\M"] defined by 

As :=a* pe, 

then [(ylj)^] satisfies the condition (13.3. 16p for the existence theorem and the forward unique- 
ness condition (j3.3.20p . 
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3.3 Generalized solution concepts for ordinary differential equations 



Proof. Since a G L1(M, it holds that for almost all t G M, we have P{t) := 

sup^gjjn |a(t,2;)| where (3 is some positive function in L-'^(M). First, observe that 



(3.3.23) sup \Ae{t,z)\ = sup 



a{t -£T,z- ey)p{T, y) dy dr 



< j Pit- I j p{r, y) dy I dr < (3{t) for almost alH G M 

where := sup^gjo^ij /jj — £s) p(r, y) dy) dr is a positive function in L-'^(M) and 
set ip{s) := J^„p{s,y)dy. Hence A^{t,x) satisfies (I3.3.16p . Using (I3.3.22p and writing 
explicitly as convolution integral, a straight-forward calculation and estimation gives 

(3.3.24) {A,{t,x)-Ae>{t,y),x-y) 

< 2a{t)\x - y\^ + 2Ci5(t)(e + e'f + 2C2 sup \a{t, z)\{e + e') + r,^,,{t)l{l + \x - y\^) 



where Ci,C2 are positive constants and r^^i;r{t) := Xe{t) — Ae'(t) with 



A.(t) 



1 



e\x - y\ 
One easily verifies that 

(3.3.25) 



{a{s, X — ez) — a{t, x — £z),x — y)p{t — s/e, z)d.s dz. 



Xeit)dt = Oie), 



hence jo' r,,,,(r) dr < A,(r) dr + J^' A.^r) dr = 0(e + e'). 
The constants Ci and C2 are derived from the estimates 

\ez-e'z'\'^(f){z)(l){z')dzdz' < [ [ii£zf-2££'{z',z) + {e'z'f)(l){z)(l){z')dzdz' 



< 



{{ezf + ee'{{z'f + (zf) + {e' z'f)^{z)^{z') dzdz' 
2{e + e'f j z'^p{z) dz =: 2{e + e'fCi 



and 



\ez-e'z\cl){z)cl){z')dzdz' < j j {\ez\ + \e' z'\)(t){z)cp(z') dz dz 

< (e + e') j \z\p{z) dz =: (e + e')C2. 
Finally we can put ()3.3.24p and ()3.3.25p together to obtain 
{As{t,x) - A^>{t,y),x - y) 



< (25(t) + -7e,.'(t))k - y\ + 2Cia{t){e + e'Y + 2C2 sup |5(t, z)\{e + e') + -7e,e'(t) 
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3 Generalized solutions of ordinary differential equations 



Setting 

a{t) := {2a{t) + l sup r,,,,(t)) 

^ e,e'e]0,l] 

je,e'{t) ■■= 2Cia{t){e + e')^ + 2C2 sup \a{t,z)\{e + e') + l-r,^e' it)- 

we obtain that Ag{t, x) satisfies ()3.3.20p . So the requirements for Theorem l3.3.6l and l3.3.9l are 
fulfilled. □ 



54 



4 A comparative study of solution concepts for 
first order hyperbolic partial differential 
equations with non-smooth coefficients 

This chapter contains a sUghtly adapted version of the joint (with Giinther Hermann) article 

m- 

According to Hurd and Sattinger in [3l| the issue of a systematic investigation of hyperbohc 
partial differential equations with discontinuous coefficients as a research topic has been raised 
by Gelfand in 1959. Here, we attempt a comparative study of some of the theories on that 
subject which have been put forward since. More precisely, we focus on techniques and 
concepts that build either on the geometric picture of propagation along characteristics or on 
the functional analytic aspects of energy estimates. 

In order to produce a set-up which makes the various methods comparable at all, we had to 
stay with the special situation of a scalar partial differential equation with real coefficients. As 
a consequence, for example, we do not give full justice to theories whose strengths lie in the 
application to systems rather than to a single equation. A further limitation in our choices 
comes from the restriction to concepts, hypotheses and mathematical structures which (we 
were able to) directly relate to distribution theoretic or measure theoretic notions. 

To illustrate the basic problem in a simplified lower dimensional situation for a linear conser- 
vation law, we consider the following formal differential equation for a density function (or 
distribution, or generalized function) u depending on time t and spatial position x 

dtu{t, x) + dx{a{t, x)u{t, x)) = 0. 

Here, a is supposed to be a real function (or distribution, or generalized function) and the 
derivatives shall be interpreted in the distributional or weak sense. This requires either to 
clarify the meaning of the product a • u or to avoid the strict meaning of "being a solution" . 

An enormous progress has been made in research on nonlinear conservation laws (cf., e.g. 
|25[ H] and references therein) of the form 

dtu{t, x) + dx{g{u{t, x))) = 0, 

where 5 is a (sufficiently) smooth function and u is such that g{u) can be defined in a 
suitable Banach space of distributions. Note however, that this equation does not include 
linear operators of the form described above as long as the nonlinearity g does not include 
additional dependence on (t, x) as independent variables (i.e., is not of the more general form 
g{t, x,u{t, x))). Therefore the theories for linear equations described in the present paper are 
typically not mere corollaries of the nonlinear theories. Essentially for the same reason we 
have also not included methods based on Young measures (cf. [25\ Chapter V]). 
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4 A comparative study of solution concepts for Erst order hyperbolic partial differential 
equations with non-smooth coefEcients 

Further omissions in our current paper concern hyperbolic equations of second order. For 
advanced theories on these we refer to the energy method developed by Colombini-Lerner 
in [9j. An overview and illustration of non-solvability or non-uniqueness effects with wave 
equations and remedies using Gevrey classes can be found in ^47j . 

Of course, also the case of first-order equations formally "of principal type" with non-smooth 
complex coefficients is of great interest. It seems that the borderline between solvability and 
non-solvability is essentially around Lipschitz continuity of the coefficients (cf . [35[ [5^ [55] ) . 
Moreover, the question of uniqueness of solutions in the first-order case has been addressed 
at impressive depth in [TO] . 

Our descriptive tour with examples consists of two parts: Section 4.1 describes concepts 
and theories extending the classical method of characteristics, while Section 4.2 is devoted 
to theories built on energy estimates. All but two of the theories or results (namely, in 
Subsections 4.1.3 and 4.2.3.2) we discuss and summarize are not ours. However, we have put 
some effort into unifying the language and the set-up, took care to find as simple as possible 
examples which are still capable of distinguishing certain features, and have occasionally 
streamlined or refined the original or well-known paths in certain details. 

In more detail. Subsection 4.1.1 starts with Caratheodory's theory of generalized solutions to 
first-order systems of (nonlinear) ordinary differential equations and adds a more distribution 
theoretic view to it. In Subsection 4.1.2 we present the generalization in terms of Filippov 
flows and the application to transport equations according to Poupaud-Rascle. Subsection 
4.1.3 provides a further generalization of the characteristic flow as Colombeau generalized map 
with nice compatibility properties when compared to the Filippov flow. In Subsection 4.1.4 
we highlight some aspects or examples of semigroups of operators on Banach spaces stemming 
from underlying generalized characteristic flows on the space-time domain. We also describe a 
slightly exotic concept involving the measure theoretic adjustment of coefficients to prescribed 
characteristics for (1 + l)-dimensional equations according to Bouchut-James in Subsection 
4.1.5. 

Subsection 4.2.1 presents a derivation of energy estimates under very low regularity assump- 
tions on the coefficients and also discusses at some length the functional analytic machinery 
to produce a solution and a related weak solution concept for the Cauchy problem. Sub- 
section 4.2.2 then compares those three theories, namely by Hurd-Sattinger, Di Perna-Lions, 
and Lafon-Oberguggenberger, which are based on regularization techniques combined with 
energy estimates. Finally, Subsection 4.2.3 briefly describes two related results obtained by 
paradifferential calculus, the first concerning energy estimates and the solution of the Cauchy 
problem for a restricted class of operators, the second is a method to reduce equations to 
equivalent ones with improved regularity of the source term. 

As it turns out in summary, none of the solution concepts for the hyperbolic partial differential 
equation is contained in any of the others in a strict logical sense. However, there is one feature 
of the Colombeau theoretic approach: it is always possible to model the coefficients and initial 
data considered in any of the other theories (by suitable convolution regularization) in such 
a way that the corresponding Cauchy problem becomes uniquely solvable in Colombeau's 
generalized function algebra. In many cases the Colombeau generalized solution can be shown 
to have the appropriate distributional aspect in the sense of heuristically reasonable solution 
candidates. 
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4.1 Solution concepts based on the characteristic Row 

4.1 Solution concepts based on the characteristic flow 

In this section wc introduce solution concepts for first order partial differential equations, 
which are based on solving the system of ordinary differential equations for the characteristics 
and using the resulting characteristic flow to define a solution. 

To illustrate the basic notions we consider the following special case of the Cauchy problem 
in conservative form 

n 

Lu ■= dtu + J2 9x,Mt, x)u) = 0, u{0) = no e V'{W), 
k=l 

where the coefficients are real-valued bounded smooth functions. The associated system 
of ordinary differential equations for the characteristic curves reads 

ik{s) = ak{s,^{s)), ^k{t) = Xk {k = l,...,n). 

We use the notation ^(s;t,x) = {(,i{s;t,x), ^n(s; i, x)), where the variables after the semi- 
colon indicate the initial conditions x = {xi, ...,Xn) &t t. We define the smooth characteristic 
forward flow 

X: [0,T] xR"^R". (s, x) ^ ^(s; 0, x) 
Note that % satisfies the relation {dx denoting the Jacobian with respect to the x variables) 

dtxit, x) = dxxit, x) ■ a{t, x) V(t, x) G [0, T] x M", 

which follows upon differentiation of the characteristic differential equations and the initial 
data with respect to t and (k = 1, . . . ,n). Using this relation a straightforward calculation 
shows that the distributional solution u G C°°([0, T]; P'(M")) to 

Lu = 0, «(0) = no G P'(M") 

is given by 

{u{t), ^) := {uo, ^{x{t, .))) VV' G <t<T. 

If there is a further zero order term b-u in the differential operator L, then the above solution 
formula is modified by an additional factor involving b and x accordingly. 

In a physical interpretation the characteristic curves correspond to the trajectories of point 
particles. This provides an idea for introducing a generalized solution concept when the 
partial differential operator has non-smooth coefficients: As long as a continuous flow can be 
defined, the right-hand side in the above definition of u is still meaningful when we assume 
uo e P'°(M"). The distribution u defined in such a way belongs to AC{[0,T];V'^{W)) and 
will be called a measure solution. 

This approach is not limited to classical solutions of the characteristic system of ordinary 
differential equations, but can be extended to more general solution concepts in ODE theory 
(for example, solutions in the sense of Filippov). Although such a generalized solution will 
lose the property of solving the partial differential equation in a distributional sense it is a 
useful generalization with regard to the physical picture. 
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4 A comparative study of solution concepts for Erst order hyperbolic partial differential 
equations with non-smooth coefEcients 

4.1.1 Caratheodory theory 

Let r > and Qt =]0, T[x]R". Classical Caratheodory theory (cf. Section [3?T|) requires the 
coefficient a = (ai, a„) to satisfy 

(i) a{t,x) is continuous in x for almost all t G [0,T], 

(ii) a{t, x) is measurable in t for all fixed x £ M" and 

(iii) sup^gjgn |a(t, x)| < (3{t) almost everywhere for some positive function (3 G L^{[0,T]). 

Then the existence of an absolutely continuous characteristic curve = ■■■,(,n), which ful- 
fills the ODE almost everywhere, is guaranteed. Note that the first two Caratheodory condi- 
tions ensure Lebesgue measurability of the composition s ^ a(s, f{s)) for all / G AC([0, T])", 
while the third condition is crucial in the existence proof. 

A sufficient condition for forward uniqueness of the characteristic system is the existence of 
a positive a G L^([0, T]), such that ((., .) denoting the standard inner product on M") 

{a{t, x) - a{t, y),x - y) < a{t)\x - 

for almost all {t, x), {t,y) £ (cf- [H Theorem 3.2.2]). As well-known from classical ODE 
theory, forward uniqueness of the characteristic curves yields a continuous forward flow 

X : [0,T] X M" ^R". (s, x) ^ ^(s; 0, x) 

It is a proper map and for fixed time x(t, .) is onto. For the sake of simplicity we assume 
a G C([0,r] X M")" and b G C{[0,T] x R"). 

Let 

hh{t,x) :=exp ^- J 6(r, x)) dr 
then u G V' (^It) defined by 

(4.1-1) {u, v)-r,'(nT) ■= J ^^0' xi^, ■))hb{t, •))©'o(Mn) 





T 

dt 



(note that u can be regarded as element in AC{[0,T];'D'^{W^)), so the restriction u{0) is 
well-defined and equal to uq £ V "^(M")) solves the initial value problem 

n 

Lu := dtu + dx^. (o^ ■ u) + bu = 0, u(0) = uq 

k=l 

on r^r, where ■ u and b ■ u denotes the distributional product defined by 
• : Cinr) X p'°(17t) ^ v'^i^r) 
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4.1 Solution concepts based on the characteristic Row 



Applying L on we obtain 

n 

fe=i 

T 

/n 
{uq, {-dtif - ^ akdx^ip + h^p){t, xit, ■))hb{t, •))x)'o(]r„) dt. 

^==1 

Set (f){t,x) := ip{t,x{t,x)) and '^{t,x) := ^{t,x) ■ hh{t,x), then we have 

n 

dt(p{t,x) = dt(p{t,x{t,x)) = {dtip + '^ak{t,x)da,^ip){t,x{t,x)), 

k=i 

and 

= dt(pit, x)hb{t, x) + (j){t, x)dthbit, x) 

n 

= {dt^ + ^ak{t,x)d:ck^){t,x{t,x)) ■ hb{t,x) - (p{t,x{t,x))b{t,x{t,x))hbit,x) 



k=l 



{dttp + ^ ak{t, x)dxi^ip - bip){t, x{t, x)) ■ hb{t, x), 



k=l 

thus 

T 



{Lu,(p)j,>^^^^ = - J {uo,dtip{t,-))j,'oQ^n)dt = - J 5t(iio,V'(*,-))x>'o(Mn)dt = 0. 


for all if G 'D{CIt)- The initial condition n(0) = uq is satisfied, since x(0,x) = x, thus 
hb{0,x) = 1. 

Remark 4.1.1. In this sense, we can obtain a distributional solution for the Cauchy problem 

n 

Pv := dtv + ^ akdx^v + cv = 0, v{0) = vq, 
k=l 

whenever a G C([0, T] x M")" and c G V'{[0, T] x M"), such that -div(a) + c G C([0, T] x M") 
and vq G We simply set b := — div(a) + c and construct the solution as above. 

In other words, such a solution solves the equation in a generalized sense, relying on the 
definition of the action of Q := J2k=i ^k^k + c on a distribution of order by 

n 
k=\ 

In case where div(a) and c are both continuous, we can define the operator Q classically by 
using the product • : P'°(f^T) x C(J^t) ^ 2^'°(^^t) as above. 
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4 A comparative study of solution concepts for Erst order hyperbolic partial differential 
equations with non-smooth coefEcients 

4.1.2 Filippov generalized characteristic flow 

As we have seen in the previous subsection, forward unique characteristics give rise to a 
continuous forward flow. But in order to solve the characteristic differential equation in the 
sense of Caratheodory, we needed continuity of the coefficient a in the space variables for 
almost all t. In case of more general coefficients a £ Lj^jM, L°°(]R"))" we can employ the 
notion of Filippov characteristics, which replaces the ordinary system of differential equations 
by a system of differential inclusions (cf. I3.3.1|) . The generalized solutions are still absolutely 
continuous functions. Again, the forward-uniqueness condition on the coefficient a 

(4.1.2) {a{s,x)-a{s,y),x-y) <a{s)\x-y\'^, x,yeW 

for almost all s G M yields unique solutions in the Filippov generalized sense. Then the gen- 
erated Filippov flow is again continuous and will enable us to define measure- valued solutions 
of the PDF (cf. [45J), as before. 

In the Filippov solution concept the coefficient is replaced by a set- valued function A : (t, x) — > 
At^x where At^x are non-empty, closed and convex subsets of M". It has to satisfy some basic 
properties (FC) (as introduced in Section 13. 2p which imply by means of Theorem 13.2.11 the 
solvability of the resulting system of differential inclusions 

i{s) G a-e., ^{t) = X, 

with ^ € AC([0, ooQ". These basic conditions (cf. section [3?2]) are 

(i) s 1-^ Ag^x is Lebesgue measurable on R for all fixed x G M", 

(ii) X I— > Ag^x is upper semi-continuous for almost all s G M, 

(iii) there exists a positive function (3 G Lj'Q^(M) such that sup^j^jgn |^s,x| < P{s) almost 
everywhere. 

There are several ways to obtain such a set- valued function ^4 from a coefficient a G L^([0, T]; L°°(M"))"', 
such that the classical theory is extended in a compatible way, i.e. the set-valued function A 
should fulfill Atg^xo •= {d-itojXo)} whenever a is continuous at {tQ,xo) G [0,oo[xR"'. 

One way obtaining a set-valued function corresponding to a a G Lj'^^,(IR; L°°(]R"'))" is by means 
of the essential convex hull essch(a). According to Definition 13.3. II its supporting function is 

Ha{t,x,w) = limesssupj^gB^(^) {a{t,y),w) 

for almost all t G J and all x, ti; G M", so 

essch(a)t,^ = {a G M" | Vu; G : {a,w) < Ha{t,x,w)}. 

Measure solutions according to Poupaud-Rascle 

Let Ooo :=]0,oo[xM". We assume a G Lj^^(M+; L°°(E"))" to be a coefficient satisfying the 
forward uniqueness criterion ()3.1.4p . Let Lu := dtu + XlILi ^^Xiidiu) and ^ be the forward 
unique solution to 

(4.1.3) i{s) G essch(a)^^5(^), ^{t) = x. 
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4.1 Solution concepts based on the characteristic how 



The map 

X:M+xM"^M", {t,x) ^ ^{t;0,x) 
is the continuous FiUppov (forward) flow. 

Definition 4.1.2 (Solution concept according to Poupaud-Rascle). Let uq G 7Wfe(M)" be a 
bounded Borel measure, then the image measure at t £ [0,oo[ is 

(4.1.4) u{t)iB):= J lB{x{t,x))duoix), 

where B CR"' is some Borel set. Themapu: [0, oo[^ 7Wfc(M")) belongs to C{[0,oo[;Mb{^'^)) 
and is called a measure solution in the sense of Poupaud-Rascle of the initial value problem 

n 

(4.1.5) Lu := 9tn + ^ 5xj.(afc • n) = 0, ii(0) = uq. 

k=l 

Note that u defines a distribution of order in T)'{^loo) by 

oo 

{u, (p)v'{n^) ■■= J {uq, (p{t, xit, a;)))25'o(iRn) dt, yip e P'(Ooo). 



The solution concept of Poupaud-Rascle does not directly solve the partial differential equa- 
tion in a distributional sense, but it still reflects the physical picture of a "transport process" 
as imposed by the properties of the Filippov characteristics. Nevertheless, in the cited paper 
of Poupaud-Rascle ( [45] ) the authors present an a posteriori definition of the particular prod- 
uct a ■ u, which restore the validity of the PDE in a somewhat artificial way. We investigate 
this in the sequel in some detail. 

Definition 4.1.3 (A posteriori definition of a distributional product in the sense of Poupaud-Ras- 
cle). Let u € V iW"') be a distribution of order and a S LJq^([0, oo[, L°°(]R"'))"', satisfying 
the forward uniqueness condition \3.1.4^ , such that there exists a continuous Filippov flow x- 
Furthermore we assume that u is a generalized solution of the initial value problem as defined 
in Ili4.1.4\ ). Then we define the product a* u = [ok ■ u)k in ^'(]0, oo[xM")" by 

oo 

(a • u,v?)x)'(n^) := {uo, J dtx{t,x)ip{t,x{t,x))dt)^'o(^n), LpeV{Q.^). 



Remark 4.1.4. Note that the product a ■ u \s defined only for distributions u that are 
generalized solutions (according to Poupaud-Rascle) of the initial value problem ()4.1.5p with 
the coefficient a. The domain of the product map (a, ti) ^ a»u, as subspace of V "^(M") x 
V^iW^) has a complicated structure: Just note that the property to generate a continuous 
characteristic Filippov flow x is not conserved when the sign of the coefficient a changes, as 
we have seen for the coefficient a{x) = sign(a;). 

Example 4.1.5. Consider problem (I4.1.5P with the coefficient a{x) := — sign(x) subject to 
the initial condition uq = 1. Then the continuous Filippov flow is given by 

X{t, x) = -{t + x)_F(-x) + {x- t)+H{x). 
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4 A comparative study of solution concepts for first order hyperbolic partial differential 
equations with non-smooth coefficients 

We have x{t, 0) = t+ - = and 

dtxit, x) = -H{-t - x)H{-x) - H{x - t)H{x) for almost all t € [0, oo[. 

The generalized solution u is defined by {u, tp) := {uq, (f){t, x)) dt, where (f){t, x) := ip{t, x{t, x)). 
We have that 



4){t,x) := { 

thus 



' (p{t,x + t) X <0,0 <t < -X 
(^(t,0) t>\x\ 
ip{t,x — t) x>0,0<t<x, 



oo oo oo 

(u,(/?)^/o(Q^) := J {uo,(f){t,x))dt = J J (f){t,x)dxdt 

-oo 

oo CO —t oo CO 

= 2 J (p{t,0)tdt + j j if{t,x + t)dxdt + j j (f{t,x -t)dxdt 

-oo t 

oo oo oo oo 

= 2j (p(t,0)tdt + J J f{t,z)dzdt + J J (p{t,z)dzdt= {l + 2tS,ip{t,-))T:)'o(s.n) 



-oo 

This generalized solution gives rise to the following product 



oo 

{{-sign{x))»{l+2tS{x)),(p)D'(n^) := {1, J dt^{t,x)(p{t,^{t,x)) dt)j,>0Q^n) 



in V'{Cloo)- Evaluating the right-hand side we obtain 

oo / oo 



{1, J dtx{t,x)ip{t,x{t,x))dt) = J I- J H{-x)H{-x-tMt,x + t)dt 

-oo \ 

oo \ 

-/i^(.Mx-,Mt,x-,)* dx. 

/ 

Since H{—x)H{—x — t) = H{—x — t) and H{x)H{x — t) = H{x — t) for t > the latter gives 
upon substitution 

oo oo oo 

(l,y" dtx{t,x)(p{t,x{t,x))dt) = - J sign(2;) J ip{t,z)dtdz, 

-oo 

hence 

(— sign(x)) • (1 + 2tS{x)) = — sign(a;). 

However, we cannot define the product if — sign(x) is replaced by +sign(a;), since the Filippov 
characteristics ^{t;0,x) are no longer forward unique and thus do not generate a continuous 
Filippov flow X- 
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Example 4.1.6. We consider the same coefficient a{x) := — sign(x) as before, but now we 
set uo := S. We obtain the generahzed solution 

oo 

{u, 'P)iy'(n^) ■= (1 ® (p{t, xit, = J v{t, xit, 0)) dt 



This enables us to calculate the product 

oo 

{{-sign{x))»S{x),(p) =-((5,y" dtx{t,x)(p{t,^{t,x))dt). 



Putting il>{x) = /q°° dt^{t, x)(p{t, ^(t, x)) dt and observe that 

oo —X 

tp{x) := J dtx{t,x)(f{t,xit,x))dt = J ip{t,x + t)dt, if x < 



and 

oo X 

ipix) = J dtx{t, x)(p{t, xit, x)) dt = - J ip{t, x-t) dt, if a; > 0. 



At X = we obtain i/j{0) = lima;_>o_ i^{x) = lim2;_>o_,_ i^{x) = 0, so it follows that (—sign) • (5 = 
0. 

Example 4.1.7. Let a{t,x) := 2H(—x), so that the Filippov flow is given by 

X{t, x) = -{x + 2t)-H{-x) + xH{x). 
We have x{t, 0) = -2t_ = and 

dtxit, x) := 2H{-x - 2t)H{-x). 
Hence dtx{i, 0) = for almost all t G [0, oo[. If uq = 1 the generalized solution is 

oo oo 

(ti, <^)x,'o(n^) ■■= j j (t){t,x)Axdt, 

-oo 

where <p{t,x) = ip{t,x{t,x)). Since 

(t){t,x)\{^<_2t} = f{t,x + 2t) 

(t>{t,x) !{_2t<x<0} = V'(*,0) 

(t>{t,x) |{0<a;} = <P{t,x), 

we obtain 

oo / —2t oo \ 

{u,^p)x,'o{p.^) = J i J (pit,x + 2t)dx + 2t(p{t,0) + J (p{t,x)dx \ dt = {l+t6,(p{t,-))j,'o(^n 

Voo / 
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hence u = 1 + 2t6{x). Again we determine the product {2H{—x)) • (1 + 2t5{x)) by 



(2F(-x).(l + 2t5(x)),v7)©. 



H{-x)H{-x - 2t)ip{t, X + 2t) dt dx 



-oo 

oo oo 



2 j j H{-x-2t)ip{t,x+2t)dtdx = 2 j j H{-z)ip{t, z) dz dt = {l(g)2H{—),ip)Ty^ 



-oo 



-oo 



We obtain {2H{—x)) • (1 + 2t5{x)) = 2H{—x). Observe that together with the result in 
Example (|iX5]) (-sign(x)) • (1 + 2t5{x)) = {2H{-x) -!)•(! + 2t6{x)) we can conclude 
that either (—1) • (1 + 2t5{x)) is not defined or the product • is not distributive. In fact , 
it is not difficult to see that (—1) • (1 + 2t5{x)) cannot be defined in this way, neither can 
l«(l + 2t(5(x)). 

Example 4.1.8 (Generalization of Example I4.1.5p . Let ci > C2 be two constants, and a G 
[ci, C2]. Consider the a(t, x) := ciH{at — x) + C2H{x — at). We set ti(x) := if x < and 
t2{x) := ^^^^ for X > 0. The unique Filippov flow is given by 

cit + X x<0,t<ti{x) 
at X < 0,t > ti(x) 

x{t, x) = \ at X = 0, 

C2t + X X > 0,t < t2{x) 

^ at X > 0,t > t2{x) 

The generalized solution of the initial value problem Lu := dtu + dx{a ■ u) = 0, u(0) = uq G 
LJ^^(R), according to Poupaud-Rascle is given by 

T ti{x) 



(n,v?)^/o(Q^) =/ (-uo, x(i,-)))x,'o(K) dt = / / uo{x)ip{t,cit + x) dtdx 



-00 



T 



00 T 



+ 



no(x)93(t, at) dt dx 



uo{x)ip{t, C2t + x) dt dx 



uo{x)ip{t, at) dt dx 



>ti(x) 







t2{x) 



T -t(ci-a) 



T 00 



no(x)93(t, cit + x) dx dt 



no(x)(/?(t, C2t + x) dx dt 



-00 



-t{c2-a) 



T 



/ t{a-C2) 



iio(x)dx Lp{t,at)dt, 



hence 



u := uo(x — cit)H{at — x) + uo(x — C2t)H{x — at) + 



\-t(ci-a) 
/ t{a-C2) 



uo{x) dx 6{x — at). 



\-t(ci-a) 
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4.1.3 Semi-groups defined by cliaracteristic flows 

Let X be a Banach space and (X't)t<^[o,oo[ be a family of bounded operators on X. Consider 
the following conditions: 

(i) So = id 

(ii) Ss o St = S^+t for all s,t £ [0, oo[ and 

(iii) the orbit maps 

auo ■■ [0,oo[ X 

t ^ St(no) 

are continuous for every uq £ X. 

If (i) and (ii) are satisfied, then we call (S^) a semi-group acting on X. If in addition property 
(iii) holds, we say {'^t)te[o,oo[ is a semi-group of type Cq. 

We briefly investigate how the solution concepts discussed in Subsections 1.1 and 1.2 fit 
into the picture of semi-group theory when the coefficient a is time-independent. First we 
return to the classical Caratheodory case: Let a £ C(M")" and assume that a suffices the 
forward uniqueness condition (j3.1.4p . This implies that the characteristic flow x '■ ^oo I^" 
is continuous and x(t, •) is onto for fixed t £ [0, oo[. Furthermore we have x(S) x(^) ^)) = 
x{s + r, x) for all x £ M" and r,s £ [0, T] with s + r £ [0, T], since a is time independent. 

Consider the inital value problem Pu = dt + X]fc=i o,kdx^u = with initial condition n(0) = 
uq £ Co(M") (i.e. vanishes at infinity). It is easy to verify that 

St:C7o(M") ^ Co(M") 

defines Co semigroup on the Banach space Co(r2oo): Note that S^ is a bounded operator on 
Co(M") for each t £ [0,oo[, as x(i,IR") = IR", so 

||Sj(uo)||oo = sup \\uQ{x{t,x))\\ = sup ||no(x(t,x))|| = sup ||no(x)|| = ||no||oo. 

xeK" xeM" xeM" 

We have that \\Tit\\ = 1 for all t £ [0, cxd[. Condition (i) and (ii) follow directly from the flow 
properties of x- The continuity condition (iii), which is equivalent to 

lim ||St(tio) - uolloo = lim sup \\uQ{x{t,x)) - uq{x)\\ = 

holds, since {x{t^x))x(^Kn- is an equicontinuous family and uq vanishes at infinity. 

Remark 4.1.9. For a coefficient a in L°^{W^)'^ we can also define a semi-group on Cq{W^) 
by St(uo) := UQ{x{t,x)), where x is the generalized Filippov flow as introduced earlier. This 
is due to the fact, that the Filippov flow has almost the same properties as the Caratheodory 
flow. 
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It seems natural to understand the solution concepts as defined by ()4.1.ip and (14.1 .411 as 
action of the dual semigroup (SJ') on 

the Banach space of finite complex Radon measures, the dual space of Cq{W^) (cf. |12[ Chapter 
4], [441 Chapter 1.10] or [53, Chapter IX. 13] for the general setting). However, in general the 
dual semi-group is not of class Co (cf. [IZ, Example 1.31]). This is only guaranteed if we start 
from a Cq semi-group defined on a reflexive Banach space. 

Nevertheless the solution concepts in (|4.1.ip and (j4.1.4p still yield the semi group properties 
(i) and (ii) with weak-* continuity replacing the strong continuity property (iii). 

The situation is much easier with Hilbert spaces, of course. We conclude with an example 
involving a discontinuous coefficient. 

Example 4.1.10. Let a £ L°^(M) such that there exist co,ci > such that ci < a(x) < C2 
almost everywhere. We want to solve the initial value problem 

Pu = dtu + a{x)d:,u = 0, u(0) = -uq G -L^(M) 

for u G AC([0,r];L2(M)) nLi([0,r];Fi(M)). 

Let A{x) = Jq a{y)~^ dy, which is Lipschitz continuous and strictly increasing (thus globally 
invertible) and observe that x(t, x) = A~^{t + A(x)) defines the (forward) characteristic fiow 
that solves 

t 

Xit,x)=x + j a(x(r,x))dr. 


Let Q := —a{x)dx with domain D{Q) := H^iM). The resolvent of Q for Re(^) > is obtained 
from the equation 

(-Q + m)^ = /, /g^' 

Upon division by a we deduce 

(4.1.6) dxv + -v = -. 

a a 

Let us first consider uniqueness: Let w G //""^(M) satisfy 

(4.1.7) dxW + -w = Q. 

a 

Since w is absolutely continuous we have 

w{x) = C exp I — 2Re(^) / dz 
\ J a{z) 

\ -oo / 

for some constant C. But w G L^(]R) if and only if C = 0, thus = 0. 
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Existence: One easily verifies that 



v{x) = {R{f,)f){x):= / exp|-/i/a(z)"idz|Mdy= J exp {- f,{A{x) - A{y))) M 



is a solution of ()4.1.6|) in AC(M). Upon substitution y z = A{x) — A{y) in the right-most 
integral we obtain 



(4.1. 



iRifJ')f)i^) = J exp(-/xz)/(x(-^,a;))d2;, 





which is the Laplace transform of -,2;)). 

We denote the kernel of the integral operator by 

M{x,y) := H{x -y)exp{-fi{A{x) - A{y))a{y)-\ 

We briefly sketch the derivation of estimates for the operator powers Rin)'' for Re(/i) > 0: 
Note that R{id)^ is an iterated integral operator of the form 



k-l 



M(-,zi)/(zi)dzi (x) 



M{x, Zk)M{zk, Zk-i) ■ ■ ■ M{z2,zi)f{zi) dzk-i . . . dz2 dzi. 



To simplify notation let z = (zi, . . . , z^), d^z = dzi . . . dzk, h{z) := exp (— /u X^fLj^ zi), and 
g{x, z) := f{x{— Sf=i zi,x)). Using the flow property of x we obtain that 

R{liff{x)= j h{z)g{z,x)d^z 

[O.ooC^ 

holds, hence by by the integral Minkowski inequality 

2 X 1/2 



(4.1.9) \\R{^,ff\\L.< 



( 



\h{z)\\g{z,x)\d^z\ dx 



1/2 

< j \h{z)\ I / \g{z,x)\^dx\ d'^z. 

[0,00 [fe 



Since 



k 

|5(z,x)|2dx = j \f{x{-jZ^ux))\^dx = j \f{y)\ 



a{x{Y.i=iZhV)) 



dy<-||/lli., 

Co 
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we conclude 



<y^-||/llL- / \hiz)\d'^ 



L2- 



The Hille-Yosida theorem ( {441 Theorem 5.2]) yields that Q generates the Cq semigroup 

St : l2(M) ^ l2(M) 

^^0 ^ uo{x{-t,x)). 

The resolvent operator /i i-^ RifJ-) (defined for Re(/x) > 0) is the Laplace transform of the 
semigroup t — > as indicated in ()4.1.8p . 

Remark 4.1.11. Since L^(M") is reflexive the dual semigroup is Cq as well and has as its 
generator the adjoint operator Q* . 

Remark 4.1.12. If we assume additional regularity on the coefficient, e.g. a € C^(M) with 
CT > 0, in Example 14.1.101 then we obtain a Cq semigroup (Sj) acting on the Hilbert space 
H^{W) with < s < a. We may then use the fact that the (square of the) Sobolev norm \\v\\1 
is equivalent to the following expression (cf. [221 Equation (7.9.4)]) 

\v{x) - v{y)\^ 



\v{x)\'^ dx + Cs 



\n+2s 



■ dx dy, 



\x - y\ 

where the constant Cs depends only on the dimension n and s. From this it can be shown 
that we may have D{Q) = H^'^^{M.) as domain of Q (this also corresponds to the special case 
of the mapping properties stated in [51, Chapter 2]). Clearly, uniqueness in the characteristic 
equation (j4.1.7|) is still valid. A corresponding variant of the estimate (|4.1.9|) for the powers 
of the resolvent operator i2(/x)'^ on ff*(M) is obtained by the following calculation (with the 
notation h and g as in Example l4.1.1Up : 



miif f){x) - {R{„f f){y)[^ 



\x-y 



l+2s 



dx dy 



< 



£1 
Co 



Re(^) 



2k 



< 



£1 

CO 



2k 




IMz)| '"^"'^^"±'^^' dz 



l+2s 

x — y\ 2 



9{x,z) -g{y,z)\' 



\x-y 



l+2s 



dx dy 




To carry out the x and y integrations we use the substitutions x' = x{~ J2t=i ^) ^ v' 
x{-YA=iZi,y) to obtain 



\g{x,z) - g{y,z)\'' 



\x-y 



l+2s 



dx dy 



l/(x(-Eti^^^))-/(x(-Eti^/,y))P 



\x -y 



l+2s 



dx dy 



\f{x')-f{y')\ 



/\|2 



lx(Eti^^^'))-x(Eti^/,y')P+2^ 



a(x(EzLi zi,x')) 



o(x') 



a(x(Ef=i^^y')) 



a{y') 



dx' dy'. 
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R R R R 



Again by J^^^^^k \h{z)\ z — jie{fi)'' conclude 




R R 



i.e., < , ll/IU- 




4.1.4 Measurable coefficients witli prescribed cliaracteristics 

This subsection discusses a solution concept according to Bouchut- James (|5j|), which is set- 
tled in one space dimension and — from the distribution theoretic point of view — can be 
considered as exotic. The basic idea is to interpret the multiplication a ■ u occurring in the 
partial differential equation as a product of a (locally finite) Borel measure u and a function 
a from the set B°° of real bounded and Borel measurable functions. 

Multiplication of Radon measures by bounded Borel functions: We may identify locally fi- 
nite Borel measures on R with (positive) Radon-measures, that is the non-negative linear 
functionals on Cc(M) ([21, Remark 19.49]). Moreover, the space P'O(M) is the space of 
complex Radon- measures, which allows for a decomposition of any u ^ T> '^(M) in the form 
n = — i/_ + z(r/_(_ — ?/_), where z/„, 17+, 17+ are positive Radon- measures. 

The product of a bounded Borel function a £ i3°°(IR) with a positive Radon measure fi is 
defined to be the measure given by 




for all Borel sets B in M. Clearly, a /i is again a locally finite Borel measure. 

The product employed in [5] is the extension of to i3°°(M) x "D "^(R) in a bilinear way, i.e. 

o : ^°°(R) X P'°(R) ^ P'°(R) 

(a, u) a+ + a_ — (a_ + a+ i/_) 

+i(a+ ??+ + a_ r?_ - (a__ 7?+ + a+ ??_)). 



Consider the following sequence of maps: 



Cfe(R) ^ i3°°(R) A L^^(R) ^ p(o)'(R) 



69 



4 A comparative study of solution concepts for first order hyperbolic partial differential 
equations with non-smooth coefficients 

where li , L2 are the standard embeddings and A sends bounded Borel functions to the corre- 
sponding classes modulo functions vanishing almost everywhere in the sense of the Lebesgue 
measure. Although we may identify C5(IR) and with subspaces of 'D"'(M) this is not true 
of i3°°(M), since A is not injective. Note that t2 o A o z.^ is injective though. The following 
example illustrates some consequences of the non-injectivity of the map A for the properties 
of the product o. 

Example 4.1.13. Let a G M and Uaix) := l,x 7^ and UaiO) := a and u = 6 e I>'°(R). 
Note that A o Li{aa) = 1 as a distribution and the standard distributional product gives 
A o Li{aa) ■ u = 1 ■ u = 5 for all a G R. On the other hand S = a6. 

The product o will be used in the solution concept for transport equations on Q,t = [0, T] x 
M with coefEcient a in B°°{a^) and solution u G B°°([0, T]; P'O(R)), i.e., u is a family of 
distributions (ii(i))te[o,T] such that {u{t), ^)x)'q(r) ^ bounded Borel function on [0,T] for all 
(p G Cc(M). The extension of the product o to this space causes no difficulty. 

The solution concept according to Bouchiut and James: A key ingredient for the solution 
concept according to Bouchut- James, is to stick to a particular representative of the coefficient 
(in the L°° sense), by prescribing the value of the coefficient a at curves of discontinuity. We 
refer to the following requirements on the coefficient a G as Bouchut- James conditions 

: Assume there exists a decomposition Qt = C U V U S such that 

(i) is a discrete subset Qt, 

(ii) C is open, a is continuous on C, 

(iii) D is a one-dimensional C^'^-submanifold of Qt, i-c, for each (tQ,XQ) G T> there exists a 
neighborhood V of (to, xq) and a parametrization of the form 1 1-^ {t,^{t)) inVOV. 
Furthermore, a has limit values for each {t, x) e V from both sides in C \ P. These 
limits are denoted by a+{t,x) and a-{t,x). 

(iv) a{t,x) G [a-{t,x),a+{t,x)] for all {t,x) G V, 

(v) for any point {to,xo) G V with neighborhood V and local parametrization ^ as in (iii), 
we have ^{t) = a{t,^{t)). 

Condition (v) prescribes the values of the coefficient a{t, x) on the curves of discontinuity 
in such a way that the characteristic differential equation holds. In this sense, a coefficient 
satisfying (i) — (v) is a piecewise continuous bounded function, where the (non-intersecting) 
curves of discontinuity can be parametrized as regular curves. 

The Bouchut-James solution concept interprets hyperbolic Cauchy problems in (1 + 1) di- 
mension as 

(4.1.10) Pu:=dtu-\-aod,,u = 0, u{0) = uo e BViod^) 
and 

(4.1.11) Lu:=dtu + da:{aou) = 0, u{0) = uq e V'^ {R) . 
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Note that P (resp. L) is well-defined on the set B°°{[0, T]; BV ioc(M)) (resp. B°°{[0, T];V'^ 
Now the main results of Bouchut-James |5j are: 

Theorem 4.1.14. Theorem 3.4] Assume that a satisfies the Bouchut-James conditions 
(i)-(v). For any uq G ^MocW there exists u G Lip([0, T]; L/„^(R)) n e°°([0, T]; SVioc W) 
solving \4.1.10\) and such that for any xi < X2 we have for all t G [0, T] 

Var/(u(t, •)) < Varj(no), 

where I :=]xi,X2[ and J :=]xi — \\a\\oot,X2 + ||a||ooi[- If in addition the coefficient a satisfies 
the one-sided Lipschitz condition 

{a{t, x) — a{t, y),x — y) < a{t)\x — y\^ for almost all {t, x), {t, y) G f^r, 

where a G L^{[0,T]), then the solution u is unique. 

Theorem 4.1.15. |3 Theorem 3.6] Assume that a satisfies the Bouchut-James conditions 
(i)-(v). Then it follows that for any uq G P'^(M) there exists u G C([0, T]; P ''(R)) solving 
\4-l-ll^ - If o- satisfies in addition the one-sided Lipschitz condition 

{a{t, x) — a{t, y),x — y) < a{t)\x — for almost all (t, x), (t, y) G ^t, 

where a G L^([0,T]), then the solution u is unique. 



We compare the solution concept of Bouchut-James with the generalized solutions according 
to Poupaud-Rascle. 

Example 4.1.16. We come back to Example l4.1.81 where a(t,x) := ciH{at—x)+C2H[x—at) 
with C2 < ci and a G [c2,ci]. Let A G L-^([0, T]) such that ci > X{t) > C2. Consider a 
representative a in B°°{0,t) of the coefficient a given by 

{ci X < at 
X{t) X = at . 
C2 X > at 



We investigate wether the distribution u given in Example 14.1.81 solves (|4.1.1ip in the sense 
of Bouchut-James. 

Let us consider the case uq = 1. Then we obtain the solution u = l + t(ci — C2)S{x — at). Note 
that the requirement that a fulfills the Bouchut-James conditions (i)-(v) forces A(t) := a for 
all t G [0, T]. Thus we have 

aov = ciH{at — x) + a5{x — at)t{ci — C2) + C2H{x — at) 

and therefore dx(a o v) = at{ci — C2)5'{x — at) — (ci — C2)d{x — at). Since dtv = (ci — 
C2)5{x — at) — at{ci — C2)6'{x — at) we deduce that u solves the Cauchy problem in the sense 
of Bouchut-James. 
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Finally, we check whether the differential equation is fulfilled, if we employ the model product 
(cf. [ISl Chapter 7] or the introduction) instead of o. Let [a ■ v] denote the model product. 
We have that 

[a • = [{ciH{at — x) + C2H{x — at)) ■ u] 

= ciH{at — x) + C2H{x — at) + cit{ci — C2)[H{at — x) ■ 5{x — at)] 

= C2t{ci — C2)[H{x — at)-5{x — at)]+ciH{at — x) + C2H{x — at) + -{ci + C2){ci —C2)6{x — at), 
hence 

dx[a-v] = {c2 - ci)6{x - at) - ^{ci + C2){ci - C2)d'{x - at), 
dtv = (ci — C2)6{x — at) — at{ci — C2)6'{x — at). 

Therefore v solves the initial value problem 

dtv + dx[a ■ v] = 0, v{0) = l, 

if a = ^(ci + C2). 

Note that the coefficients H{—x) (with ci = 1,C2 = 0,a = 0), —H{x) (with ci = 0, C2 = 
— l,a = 0), and — sign(x) (when ci = 1,C2 = —I, a = 0) are included as special cases of the 
example presented here. In case the coefficient reads 

a(x) := H{—x) 

the unique solution in the sense of Bouchut-James is given by 

u=l + t5. 

It has been shown in [291 Theorem 5] that no distributional solution exists in this case when 
the model product is employed. 

4.2 Solutions from energy estimates 
4.2.1 Direct energy estimates 

We briefly review the standard techniques of energy estimates for the initial value problem 



Letg G [2,00]. We assume that / G L^{[{},T]-L'^{W)), a = (ai,...,an) G ^^([0, T]; l^i'''(M")) 
with real components, c G L^([0, T]; L'^(M")) and in addition 



n 



(4.2.12) 




in ]0,T[xR' 



n 



n(0) = no G L2(M"). 



(4.2.13) 



-div,(a)-c gLH[0,T];L-(M")). 
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Example derivation of an energy estimate 

We browse through the typical steps that lead to an estimate in the norm of -L°°([0, T] ; (M**)) 
for any 

u e AC{[0,T];L'^{W)) D L°°([0, T]; M^^'f (M")) 

with p £ [2, oo] such that ^ + | = ^ in terms of corresponding norms for ■u(O) and Pu. 
We write P = dt + Q with Q := Ylk=i ^kii, x)dxj. + c(i, x) and observe that 

P«eL^([0,r];L2(R")) 

holds since dfU G ^^([O, T]; L2(M")) and Qu G Li([0, T]; ^^(M")) (the latter follows from the 
facts that d:,.u{t, .) G andL«-i7 C when2/p+2/g = 1). Hence r Re(((Pn)(r), u(r))o) 
is defined and in L^([0,T]). Furthermore, the map t ^ \\u{t, •)||o is continuous. 

We put 

r 

Mr) := ||^div^(a(r,-)) -c(r,-) Hoc and A(r) := 2 ^ ds > (r G [0,T]). 



By assumption, /t G Li([0,r]) and A G ^C([0,r]). 

The standard integration by parts argument gives the Garding-type inequality 
(4.2.14) ^((Qu(r),u(r))o + (t^(T),Qu(T))o) =Re((Qu(r),u(r))o) > -h{T)\\u{T)\\l, 



and thus 



y"e-^('^)Re(((P^x)(r),^i(r))o) dr 



T T 

= ^ y" e-^w|;||t.(r)||2dr + j e-^('^)Re((Qu)(r), u(r))o dr 



> le-'^^-'^Hrnl - \\\uml - J (h{r) - ^) e-^W||Mr)||g dr. 



Therefore 



T 

{r)\\l < ||n(0)||2 + 2|e-^W||(Pn)(r)||o||n(r)||odr 



T 

< ||^(0)||2 + 2sup (e-W||^(^)||o) /e-W||(p^)(^)||odr, 
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where we may take the supremum over r G [0, t] on the left-hand side and thus replace r by 
t on the right-hand upper bound. A simple algebraic manipulation then gives 



sup ||e-^W/2n(r)||o- / e-^^'ymPu){r)\\odr < ||^x(0)||o+ / e-^'^^'ymPu){r)\\o dr 



Upon removing the squares and multiplying by exp(A(t)/2) we obtain the following basic 
inequality. 



Energy estimate: 



(4.2.15) sup ||n(r)||o < exp( / da) • ||n(0)||o + 2 exp( / h{a)da)- / ||(Pu)(r)||o dr 

refO.i] J J J 



t / t 



eMj h{a)da) ||7x(0)||o + 2 J \\{Pu){r)\\o dr 
V y 



We recall that the exponential factor depends explicitly on the coefficients a and c via h{r) = 
||idiv^(a(r, •)) - c(r, •)||oo. 

Note that this derivation of an energy estimate relied on the Garding inequality (I4.2.14|) . 

Example 4.2.1 (Failure of the Garding-inequality (j4.2.14p ). Let a g]1/2, 1[ and define 
a: M M by a{x) := 1 + when x < 1, and o(x) := 2 when x > 1. We have a G 
C-(M)\Lip(M). 

Let Q: H^{R) L'^{R) be the operator defined by {Qv){x) := a{x)v'{x) for all v G H^{R). 
Note that compared to the general form of the operator Q in the derivation of the energy 
estimate above we have here c = 0, a G C°°([0,r]; W'^''^{R)) but diva/2 - c = a'/2 L°°(M). 

Since Q is time independent, inequality ()4.2.14p with some h G L^([0,T]) (not necessarily of 
the form given above) would imply 

3C G M, V-y G C^{R) : Re{{Qv,v)o) > -C\\v\\l. 

We will show that there is no constant C G M such that the latter holds. Thus (j4.2.14|) cannot 
hold for Q (for any h G L''{[0,T])). 

Let p G C"'^(M) be symmetric, non-negative, with support in [—1, 1], ||p||o = 1, and such that 
p'{x) < when < x < 1. We define t'e(x) := e~^^'^p{x/e) (x G M, e > 0). Then clearly 
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Ve G C H\R) and \\ve\\o = 1 for all e > 0, but 

{QVe,Ve)o = / a{x)v'^{x)Ve{x) Ax 



L OO 

v'^{x)ve{x) dx + j x°^v'^{x)ve{x) dx + J v'^{x)vs{x) dx 



-oo 



=0 

1 

_a— 1 / n / 



= e°-' J z°p'{z)p{z)dz ^ -OO (e^O). 



We remark that even for a £ Cl{M.) \ Lip(M) the Garding inequality may fail as well: for 
example, with a{x) := —xlog\x\p{x) we have a G Cj(M) n W^^'^iW) for all q £ [l,oo[, but 

1 1 
{Qvs,Ve)o = -"^ j pie z)z log \ez\p{z)p{z) dz < 2|loge| j zp{ez) p {z)p{z) dz ^ -oo, 





smce 



lime^o Jo zp{ez)p'{z)p{z) dz = p(0) zp'{z)p{z) dz < 0. 



Remark 4.2.2. (i) Let Q* denote the formal adjoint of Q with respect to the inner 
product (on x-space) . Due to our regularity assumptions on a and c we have for any ip £ 
(since a is real) 

n 

Q*ip = '^{-ajdc.^ip) + (c- dwx{a))ip, 
i=i 

where the new coefficients —a, respectively c — diva;(a), in place of a, respectively c, satisfy 
the exact same regularity assumptions, including the condition 



idiv,(-a) - (c-div,.(a)) = _c e L\[0,T]; (R^)). 

Thus the basic energy estimate (|4.2.15p applies to zizdt + Q* as well. In particular, the function 
h in the exponential factor occurring in the energy estimates is the same for Q and Q* . 

(ii) Although the method of derivation discussed above relied on a Garding-type inequality, 
it seems that in essence energy estimates are, in a vague sense, a necessary condition for 
a hyperbolic equation to hold in any meaningful context of "suitable Banach spaces of dis- 
tributions". In other words, whenever a hyperbolic differential equation can be interpreted 
directly in terms of such Banach spaces it allows to draw consequences on combinations of cor- 
responding norms of any solution. For example, if the operator Q above generates a strongly 
continuous evolution system on some Banach space, then basic norm estimates for solutions 
follow from general principles of that theory (cf. [44 ^ [l9] ) . 



On the other hand, energy estimates are widely used to establish existence of solutions to 
(j4.2.12p by duality and an application of the Hahn-Banach theorem. We recall the basic steps 
of such method in the following. 
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Existence proof based on the energy estimate 

Let Rt ■■= {it,x) E R"+^ \t<T}. By abuse of notation we denote the trivial extension of 
a function v G C^{Rt) by zero for t > T again by v. Then £ := {/ e C°°([0,r] x M'") | 
3v G C^{Rt) with / = {-dtv + Q*v) |[o,t]xR"}- For < t < T and t; G C^iRr) we use the 
notation w{t) := v{T — t) and g{t) := {—dtV + Q*v){t). Then we have 

{dt + Q*{T-t))w{t) = g{T-t) 
w{0) = 

and an apphcation of (|4.2.15p (with Q* in place of Q; cf. Remark I4.2.2l fi) above) yields 

(T \ T T 

[ h{a)da] [ \\{-dt + Q*v){T -r)\\odr = [ H^WIIodr. 

J J J 

0/0 

We may deduce that for / G Li([0, T]; L2(]R")) and v G C^{Rt) 

T T 

j (/(r),t;(r))odr + (7Xo,t'(0))o < j ||/(r)||o|b(r)||o dr + ||^xo||o ||^^(0)||o 



T 

<C sup \\w{r)\\o<CCh [ ||5(r)||odr, 

rGfO.Tl J 







where C depends on / and uq. Therefore the assignment g = {—dtv + Q*v) |[o,t]xR"'^ 

/^^(/(r), f (r))o dr+(uo5 '^(0))o defines a conjugate-linear functional : £ ^ C on the subspace 
£ of L1([0,T];L2(M")) such that \v{g)\ < 

supo<r<T ll5'(^)l|o- Hahn-Banach extension of u 
yields a conjug ate-linear functional v' : L'^{[0,T]; L^{W')) C with the same norm estimate. 

Since L\[0,T]; L^{W)y ^ L°°([0, T]; there is tt G L°°([0, T]; L2(]R")) such that 

u'{g) = {u,g) for all g G ^^([O, T]; L2(IR")). When applied to g = {-dtv + Q*v) |[o,t]xR" 
with V G C^{Rt) we obtain 

T 

(4.2.16) / {u{t), -dtv{t) + {Q*v){t))^ dt = {u, {{-dt + Q*)v) l[o,T]xRn> 





T 

= J {fit),v{t,.))odt + {uo,v{0))o. 



l\/lodel discussion of tlie weak solution concept 

Case of smootli symbol: If the coefficients of Q (and thus of Q*) are C°° then the above 
identity implies that u is a distributional solution to the partial differential equation Pu = f 
in ]0,r[xM". In fact, with ip G Cc°°(]0, r[xM") in place of v we have 

T 

{{dt + Q)u, ^) = {u, i-dt + Q*)ip) = j (fit), <A(t, .))o dt = (/, if). 
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Moreover, since Qu G L°^{[0,T];H-^{W)) the differential equation implies that 

dtu = f-Qu £L\[Q,T]-H-\W)) 

and thus u S AC([0, T]; In particular, it makes sense to speak of the initial value 
u(0) E D'(M"'). Integrating by parts on the left-hand side of (|4.2.16p (now reading (j4.2.16p 
from right to left, and duality brackets in appropriate dual pairs of spaces) yields for any 

T T 

j {f{t),v{t))odt+{uQ,v{Q))Q = j {dtu{t) + Qu{t),v{t)l dt-{u{T),y{T))o + {u{0),v{0))o, 

{f{t),v{t,-))o =0 

hence u{0) = uq 

Of course, uniqueness of the solution as well as more precise regularity properties can be 
deduced in case of coefficients: For any s G M, / G Li([0, T]; F*(R")), and uq G i?'(IR") 
the solution u is unique in the space C{[0,T]; H^{R"-)) (cf. [22, Theorem 23.1.2]). 

Case of non-smooth symbol: The weaker regularity assumptions made above imply Q*v € 
Li([0,r];L2(M")) for ah v G C^{Rt)- We may thus define Qu G P'(]0,r[xM") by putting 

{Qu, if) := {u\Q*^)o yip G C^{]0, T[xM"). 

Then equation (j4.2.16p can be read as an equation in D'(]0, T[xM"'), namely 

{dtu + Qu, if) = if, if) G C7-(]0, r[xM"). 

Furthermore, we can again show that the inital datum is attained: Note that in Qu = 
Yli O'jdxjU each term can be interpreted as a multiplication of functions in L^{[Q,T];H^{W)) 
with distributions in L°°([0, T]; F-^(M")) (since u G L°°([0, T]; ^^(R"))) in the sense of the 
duality method (cf. [HI Chapter II, Section 5]). Applying Proposition 5.2 in [l3] to the 
spatial variables in the products then yields Qu G L^{[0, T]; W~^'^{W^)). Reasoning similarly 
as above, the differential equation then gives 

dtu = f-Qu e L\[0,T];W-^'HW')), 

which implies u G AC{[0,T]; VF~^'^(M")) and further also that n(0) = uq- 

Again higher regularity of u with respect to the time variable, namely u G C([0, T]; L^(M"')) 
can be shown by means of regularization and passage to the limit (e.g., similarly as in [U 
proof of Theorem 2.8]). 

4.2.2 Regularization and energy estimates 

Several advanced theories make use of regularization techniques or concepts at crucial steps 
in their construction of solutions. Some of these theories succeed by regularization and a 
careful passage to the limit via energy estimates (as with Hurd-Sattinger and Di Perna-Lions 
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theories presented below). Others even base their solution concept on a further generaliza- 
tion of the weak solution concept beyond distribution and measure spaces and still obtain 
existence of solutions essentially from asymptotic stability of energy estimates (cf. the Lafon- 
Oberguggenberger theory below). 

We introduce the following notation for partial differential operators that will be used in the 
sequel 

n 

(4.2.17) Pu := dtu + ^akdx^u + cu 

k=l 

n 

(4.2.18) Lu := dtu + '^d^^{aj u) + bu. 
Hurd-Sattinger theory 

We give a brief summary of the results from the first part in Hurd-Sattinger's classic paper 
|34j . We consider the Cauchy problem for the operator L on the closure of the domain 
n :=]0,cx)[xM". 

Definition 4.2.3. Let f G L|^^(J7) and aj (j = l,...,n) as well as b ^ Lf^^{Q). A weak 
solution in the sense of Hurd-Sattinger of the partial differential equation 

Lu = f on Q 

with initial condition uq £ L^^^ is a function u G Ly^^{Q.) such that for all (p £ C^(M"+^) we 
have 



(4.2.19) 




= //(*,.M,t..)d(M) + /„„w,(o..)d.. 

Note that if all coefficients are C°° functions then a solution in the above sense solves the 
partial differential equation on Q in the sense of distributions. 

Theorem 4.2.4. Let aj (j = 1,. . . ,n), b, and f belong to L'j^^{n) and uq £ Lf^^. Assume, 
in addition, that the following conditions are satisfied: 

(i) There exists ci > such that for almost all (t, x) £ Q: ak{t, x) < ci (k = 1 . . . ,n). 
(a) There exists a function fj. £ Ll^^{[0,oo[), fi>0, such that b{t,x) > —ii{t) for almost all 

{t,x) £ n. 

(Hi) For each k £ {l,...,n} there exists < fik £ Ll^^{[0,oo[) such that for almost all 
{t,x) £ Tl 

akit,x)-ak{t,xi,...,Xk-i,r,Xk+i...,Xn) ^ , , , „ ^ ^ 

> —IJ'kyt) for almost all r £ K. 

Xk-r 
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4.2 Solutions from energy estimates 
Then there exists a weak solution u G L'f^^{Q) to Lu = f with initial condition uq. 
Concerning the meaning of condition (iii) in Theorem 14.2.41 we mention two aspects: 

• In one space dimension we obtain > —fj,i{t), which resembles a one-sided 
Lipschitz continuity condition in the x variable (apart from the fact that need not 
be finite or defined for all t). In particular, it excludes jumps downward (seen when 
going from smaller to larger values in the x argument). 

• Heuristically — replacing difference quotients by partial derivatives — condition (iii) 
can be read as diva(t, x) > fJ-kii), thus giving an lower bound on the divergence 
of a. We observe that upon formally applying the Leibniz rule in the operator L we 
cast it in the form P as in (14.2. 17p with c = div a + b. In combination with condition 
(ii) of Theorem 14.2.41 we obtain that |divo — c = — (|diva + b) has an upper bound 
(uniformly in x), which can be considered a substitute for condition (j4.2.13p used in the 
derivation of direct energy estimates in Subsection 2.1. 

Remark 4.2.5. Hurd-Sattinger ([34J) also give a uniqueness result for first-order systems 
in case of a single space variable and 6 = 0. For scalar equations the hypotheses require 
condition (i) to be strengthened to boundedness from above and from below and condition 
(iii) to be replaced by a Lipschitz property with an upper bound instead; in particular, no 
jumps upward are possible. 

Example 4.2.6. For the operator L in one space dimension and coefficients o(x) = sign(x) 
and b = 0, the Poupaud-Rascle theory is not applicable (as mentioned in [351 Section 1, 
Example 2]), but Hurd-Sattinger theory ensures existence of weak solutions, if the initial 
value belongs to L'^^^. 



Di Perna-Lions theory 

The weak solution concept introduced by Di Perna-Lions in [T3] for the Cauchy problem for 
the operator P on a finite-time domain [0, T] x M" can be interpreted in the following way. 

Definition 4.2.7. Zei T > 0, 1 < p < oo, | + | = 1, / G L^([0, T]; LP(M")), at G 
Li([0,r];Lf^^(M")) (k = l,...,n), and c G ^^([O, T]; Lf^^(M")) such that 

div(a)-c GL\[0,T];Ll^{R-)). 

A function u G L°°([0, T]; LP(M")) is called a weak solution in the sense of Di Perna-Lions 
of the partial differential equation 

Pu = f on]0,T[xW 
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with initial value uq G LP{W^), if 

T ^ 
(4.2.20) j j u{t,x)(^-dt^{t,x)dx-^ak{t,x)d,^^^{t,x)^dxdt 

M" 

T 

~^ j j u{t,x)(^ — diva{t,x) + c{t,x)')(f{t,x)DCmxdt 

R" 

T 

~ J J fit:x)(f{t,x)dxdt + J uo{x)(f{0,x) dx 



holds for all ip G C°°{[0,T],W) with compact support in [0,r[xM". 

Clearly, in case of coefficients we obtain a distributional solution of the partial differential 
equation in ]0,r[xM". 

Theorem 4.2.8. Existence of a weak solution u G L°°([0, T]; ^^(M")) in the sense of and 
with assumptions as in Definition \4-2.7\ is guaranteed under the additional hypothesis 



-div(a)-cG Li([0,r];L°°(M")), if p > 1, 
P 

div(a),c G Li([0,r];L°°(M")), if p = 1. 

Remark 4.2.9. Uniqueness holds in general under the additional hypotheses that c, div(a) G 
Li([0, T]; (M")), and for j = 1, . . . , n also aj G L^{[0, T]; w/^f (M")) as weh as 

G L\[0, T];L'^{M.'')) + L\[0, T]; L°°(M")). 



1 + X 



Example 4.2.10 (Hurd-Sattinger applicable, but not Di Perna-Lions). Note that with a 
single spatial variable boundedness of div (a) = a' implies Lipschitz continuity. Hence, if 
a G (M) is not Lipschitz continuous but satisfies the one-sided Lipschitz condition in Hurd- 
Sattinger 's existence Theorem 14.2.41 (iii). then a weak solution in the sense of Hurd-Sattinger 
to the problem 

dtu + d^{au) = / G L^{^^), u |t=o= ^^o e L^' 



is guaranteed to exist, whereas the general statement of DiPerna-Lions' existence theory 
(Theorem 14.2.81 with p = q = 2) \s not applicable to the formally equivalent problem 

dtu + ad^u + a'u = f £ L^iR^), u \t=o= uq G L^{R). 

Example 4.2.11 (Di Perna-Lions applicable, but not Hurd-Sattinger). Let < cr < 1 and 

consider the identical coefficient functions ai = a2 G C'^,comp(l^^) (i-^-i comsupporteduported 
functions in C^(IR^)) given by 

ai{x, y) = a2{x, y) = - '^{x - y)\ x{x, y), 
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where x S I'(M^) such that x = 1 near (0, 0). Note that oi is not Lipschitz continuous, since 
for X > but X sufficiently small the difference quotient 



oi(x, 0) — ai(0, 0) X 



X a 



is unbounded as x — > 0. In particular, the latter observation shows that the Hurd-Sattinger 
existence theory is not applicable (condition (iii) in Theorem 14.2.41 is violated) to the Cauchy 
problem for the operator 

Lu = dtu + dx{aiu) + dy{a2u). 

On the other hand, we can show that with a = (ai, 02) the DiPerna-Lions existence theory is 
applicable to the Cauchy problem 

dtu + aidxU + a2dyU + {diYa)u = f e L^([0, T]; ^^(M^)), u \t=o= uq E 
To begin with, we observe that 

dxai{x,y) = dxa2{x,y) = - - -{x - y)"^ d^xix , y) 



{x - y) 



a 



yields 



dytti {x,y) = dya2 {x,y) = "^^^'f^l^ - - (x - y)+ dyx{x, y) 

{x-y). o 



diva(x,y) = - y)+ divx(x, y) e C^%mp(I^^)- 



Note that in the notation of Definition 14.2.71 and Theorem 14.2.81 we have c = diva E L°°(]R^) 
(and time-independent). Therefore, the basic assumptions for the solution concept to make 
sense as well as the hypotheses of the existence statement are clearly satisfied. 

As for uniqueness, we remark that all the conditions mentioned in Remark 14.2.91 are met if 
and only if cr > 1 /p. 

Remark 4.2.12. We mention that with coefficients as in the above example, the system of 
characteristic differential equations has forward-unique solutions, hence the Poupaud-Rascle 
solution concept for measures is also applicable. 



Lafon-Oberguggenberger theory 

The theory for symmetric hyperbolic systems presented in [31] by Lafon-Oberguggenberger 
allows for Colombeau generalized functions as coefficients as well as inital data and right-hand 
side. Thus we consider the following hyperbolic Cauchy problem in M""*"^ 

n 

(4.2.21) Pu = dtu^^ a^d^M + cu = f 

(4.2.22) u \t=o= uo, 
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where aj (j = l,...,n), c are real valued generalized functions in Q{W^^^) (in the sense 
that all representatives are real valued smooth functions), / G Q{W^^^), and initial value 
no G a(M"). 

The coefficients will be subject to some restriction on the allowed divergence in terms of 
e-dependence. A Colombeau function v G Q{M.'^) is said to be of logarithmic type if it has 
a representative (vg) with the following property: there are constants N G N, C > 0, and 
1 > r] > such that 

(J 

sup |we(y)| < A^log (— ) < e < ?7 . 

(This property then holds for any representative.) By a suitable modification of [421 Propo- 
sition 1.5] it is always possible to model any finite order distribution as coefficient with such 
properties (in the sense that the Colombeau coefficient is associated to the original distribu- 
tional coefficient). 

Theorem 4.2.13. Assume that aj and c are constant for large \x\ and that dxf.aj (k = 
1, . . . ,n) as well as c are of logarithmic type. Then given initial data uq G t/(M") and right- 
hand side f G Q{M.^'^^), the Cauchy problem |^.g.^j[ )-( (^.^.^^ ) has a unique solution u G 
g(R"+i). 



We also mention the following consistency result which shows that Colombeau theory 
includes the classically solvable cases: If we assume that the coefficients Uj and c are C°° then 
we have the following consistency with classical and distributional solutions (cf. [39j) 

• If / and no are functions then the generalized solution n G ^(M"'*'^) is equal (in Q) 
to the classical smooth solution. 

• If / G L^(IR; //*(M")) and no G i/'*(M") for some s G M, then the generalized solution 
u G ^(M"+^) is associated to the classical solution belonging to C{R; H'^{W^)). 

Example 4.2.14. Consider the (1 -|- l)-dimensional operator 

Lu = dfU + dx{H{—x)u). 

Since the coefficient (of the formal principal part) has a jump downward neither Hurd- 
Sattinger nor Di Perna-Lions theory is applicable. In fact, it has been shown in [291 Section 
2] that none of the distributional products from the coherent hierarchy (cf. [33] and the in- 
troductionary section) applied to H{—x) ■ u is capable of allowing for distributional solutions 
of the homogeneous Cauchy problem for arbitrary smooth initial data. 

Recall from Section 1 that measure solutions according to Bouchut-James exist for the corre- 
sponding Cauchy problem, if the Heaviside function (usually understood as a class of functions 
in L°°) is replaced by the particular Borel measurable representative with value at x = 0. 
For example, the initial value no = 1 then yields the measure solution n = 1 + t6{x) in the 
sense of Bouchut-James as seen in Example 14.1.161 

However, Colombeau generalized solutions are easily obtained — even for arbitrary general- 
ized initial data — if the coefficient H{—x) is regularized by convolution with a delta net of 
the form Peix) = log(l/e)p(x log(l/e)) (0 < e < 1), where p G C^(M) with J p = 1. Let a 
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denote the class of this regularization in the Colombeau algebra Q, then the operator L may 
now be written equivalently in the form 

Pu = dtu + adxU + a'u, 

where a' ^ 6 and u £ Q. Due to the logarithmic scale in the regularization the hypotheses 
of Theorem 14.2.131 are satisfied and the corresponding Cauchy problem is uniquely solvable. 
Moreover, for most interesting initial data (e.g. Dirac measures or L^^^) weak limits of the 
Colombeau solution u are known to exist and can be computed (cf. \29\ Section 6]). In 
particular, for the initial value uq = 1 we obtain the measure solution u = 1 + t6{x) as such 
a distributional shadow. 

Remark 4.2.15. (i) The basic results of Lafon-Oberguggenberger have been extended to the 
case of (scalar) pseudodifferential equations with generalized symbols in [57]. Special cases 
and very instructive examples can be found in and an application of Colombeau theory 
to the linear acoustics system is presented in [32]). 

(ii) Colombeau-theoretic approaches allow for a further flexibility even in interpreting distri- 
butional differential equations with smooth coefficients. For example, in [8] the concept of 
regularized derivatives is used, where partial differentiation is replaced by convolution with 
the corresponding derivative of a delta sequence. When acting on distributions this concept 
produces the usual differential operator actions in the limit. When considered as operators in 
Colombeau spaces, one can prove (cf. [HI Theorem 4.1]) that evolution equations with smooth 
coefficients all whose derivatives are bounded have unique generalized function solutions for 
initial data and right-hand side in generalized functions. In particular, famous examples like 
the Lewy equation become solvable and Zuily's non-uniqueness examples become uniquely 
solvable then. 

4.2.3 Paradifferential techniques 
Energy estimates 

Bony's paradifferential calculus has been successfully applied in nonlinear analysis and, in 
particular, to regularity theory for nonlinear partial differential equations. An ingredient in 
such approaches is often a refined regularity assessment of corresponding linearizations of the 
differential operators involved. A recent account of Metivier's methods and results of this 
type can be found in [4, Subsection 2.1.3], or with more details on microlocal properties in 

m- 

Let s E M and H^{W^) denote the Sobolev space i7*(R"') equipped with the weak topology. 
We consider a differential operator of the form 

n 

P^{x,t;dt,dx) ■■= dt + ^aj{v{x,t))dj, 
j=i 

where aj e C°°(R) (j = 1 . . . , n) and v E L°^([0, T];H'{W)) n C7([0, T];H^{W)) such that 
dtv G L°°([0,T];F^-i(M"))nC([0,r];i74-i(IR")). 
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4 A comparative study of solution concepts for Erst order hyperbolic partial differential 
equations with non-smooth coefEcients 

Remark 4.2.16. Not all hyperbolic first-order differential operators with coefficients of reg- 
ularity as above can be written in the special form of P^. In fact, this amounts to writing any 
given list wi, . . . , Wn of such functions as Wj = aj o v {j = 1 . . . , n) with aj G (M) and v as 
above. The latter is, in general, not possible, which can be seen from the following example: 
consider the Lipschitz continuous functions wi{t) = \t\ and W2{t) = t; if = ai o v and 
W2 = 0,2 o V with a Lipschitz continuous function v, then v is necessarily non-differentiable at 
0; on the other hand 

1 = w'2i0) = lim{a2{v{h)) - a2{v{0)))/h = lim a'2{^{h)){v{h) - v{0))/h, 

where ^(/i) lies between v{0) and v{h); hence a'2{^{h)) — > a'2{v{0)) and the second factor 
{v{h) — v{0))/h stays bounded, but is not convergent; in case 02(^(0)) = we obtain the 
contradiction 1 = 0, in case 02(^(0)) 7^ we have a contradiction to convergence of the 
difference quotient for W2- 

The key technique in analyzing the operator is to replace all terms aj{v)dj by Ta^(^y)dj, 
i.e., partial differentiation followed by the para-product operator (t,), and then employ 
estimates of the error terms as well as a paradifferential variant of Garding's inequality (cf. 
[H Appendix C.3-4]). This leads to the following result. 

Theorem 4.2.17 (gl Theorem 2.7]). 7/s > § + 1, then for any f € L°°{[0,T]; H^W)) D 
C{[0,T];H^{M.'')) anduQ G H'{W) the Cauchy problem 

PvU = /, U \t=0= Uo 

has a unique solution u G L'^{[Q,T]]H^{W)). Moreover, u belongs to C([0, T]; and 
there are constants K,'y,C > such that u satisfies the energy estimate 

t 

\\u{t)f^<Ke^'\\um\l + C I e^('-^^\PMr)\\ldT. 



Improvement of regularity in one-way wave equations 

We briefly recall some basic notions and properties concerning symbols with certain Holder 
regularity in x and smoothness in a la Taylor (cf. [51]). 

Definition 4.2.18. Letr > 0, < 6 < 1, andm G M. >1 continuous function p : R"-x]R" ^ C 
belongs to the symbol space C^S^^, if for every fixed x £ M" the map ^ ^ Pix,^,) is smooth 
and for all a G Nq there exists Ca > such that 

\d^p{x,0\ < Ca(l + |e|)™'l°l Vx,^ G 

and 

\\d^Pi.,o\\c: - ^"^^ + iei)"^-i"i+^' ve G M". 

Basic examples are, of course, provided by symbols of differential operators Yl o.^d'^ with 
coefficient functions Oq, G CI {\a\ < m) or any symbol of the form p{x,^) = a{x)h{x,(,), 
where a G and /i is a smooth symbol of order m. 
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4.2 Solutions from energy estimates 



Symbol smoothing: By a coupling of a Littlewood-Paley decomposition in ^-space with con- 
volution regularization in j;-space via a 5-dependent scale one obtains a decomposition of any 
symbol p £ C^S^^g in the form 

p = p^+p\ where p^ G S^s and p^ G C^S'^^''^. 

Observe that is C°° and of the same order whereas p^ has the same regularity as p but is 
of lower order. 



Mapping properties: Let < 6 < I and -{1 — 6)r < s < r. Then any symbol p G C^S^ 
defines a continuous linear operator p{x,D) : if^+'^(M") — > ff*(R"). 

Elliptic symbols: p G C^S^g is said to be ehiptic, if there are constants C, i? > such that 

>c(i + ieir vcGMMei >i2. 

One-way wave equations result typically from second-order partial differential equations by 
a pseudodifferential decoupling into two first-order equations (cf. |i50j Section IX. 1]). For 
example, this has become a standard technique in mathematical geophysics for the decoupling 
of modes in seismic wave propagation (cf. [i8\). The corresponding Cauchy problem with 
seismic source term / G C°°([0, T]; //*(M")) (with s G M) and initial value of the displacement 
no G i^^+^IR") is of the form 

(4.2.23) dtu + iQ{x,D)u = f 

(4.2.24) u\t=o = no, 

where Q has real- valued elliptic symbol q G C^S"^ with r > s. 
Lemma 4.2.19. If q £ C^S^^q is elliptic, then G S^g is also elliptic. 

Proof. By ellipticity of q and the symbol properties of q^ there are constants Ci,C2,R > 
such that 

ci{i + lei)™ < \q{x,o\ < \qHx,o\ + \q\x,o\ < \qKx,o\ + C2{i + leir-*^' 

holds for all x,^ G M" with |C| > -R > 0. Therefore 

\qHx,o\ > (Ci - ^2(1 + i^i)-^''')(i + 1^1)"' > c(i + lei)"^ vx,e G M", lei > R' 

for suitably chosen constants C and R' > 0. □ 

Let < 6 < 1. We have the decomposition q = q^ + q^, where G 5"^ ^ and q" G CS'J^^''. By 
Lemma [4. 2. 191 = q^{x,D) is elliptic and thus possesses a parametrix E'^ G S^j. 

We have 

{dt + iQ)E^f = {dt + iQ^ + iQ^)E^f = dtE^f + iQ^E^f + iQ^E^f = dtE^f + f + iR^f + iQ^E^f, 
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4 A comparative study of solution concepts for Erst order hyperbolic partial differential 
equations with non-smooth coefEcients 

where i?" is a regularizing operator. Therefore 

(dt + iQ)iu - E^f) = -dtE^f - iR^f - iQ^E^f =: /, 

where the regularity of the right-hand side / can be deduced from the following facts 

dtE^f £ C°°([0,T];i7"+i(M")), R^f e C7°°([0, T]; i/°^(M")), Q^E^f £ C'^{[0,T]; H'+^'^R'')). 

Hence f£ C°°([0, T]; ('^'^'i)(]R")). 

If we put w = u+E^f and wq := uo+-E''*/(0), then the original Cauchy problem (|4.2.23ti4.2T^ 
is reduced to solving the Cauchy problem 

dtw + iQ{x,D)w = f, w\t=o=wo, 

where the spatial regularity of the source term on the right-hand side has been raised by 
min((5r, 1). 

Remark 4.2.20. In case of a homogeneous (1 -|- l)-dimensional partial differential equation 
the precise Holder-regularity properties of classical as well as generalized solutions have been 
determined in [551 Section 3]. 



86 



5 Microlocal analysis of generalized pullbacks of 
Colombeau functions 



This chapter is based on the article [19j . 

The pullback of a general distribution by a C°°-function in classical distribution theory, as 
defined in [231 Theorem 8.2.4], exists if the normal bundle of the C°°-function intersected 
with the wave front set of the distribution is empty. These microlocal restrictions reflect also 
the well-known fact that in general distribution theory one cannot carry out multiplications 
unrestrictedly, since the product of two distributions can formally be written as the pullback 
of a tensor product of the two factors by the diagonal map 5 : x ^ {x,x). 

Generalized functions in the sense of Colombeau extend distribution theory in a way that it be- 
comes a differential algebra with a product that preserves the classical product • : C°° x C°° — > 
C°°. In addition Proposition 1.2.8] states that the Colombeau algebra of generalized func- 
tions allows the definition of a pullback by any c-bounded generalized function. The classical 
concept of a wave front set has been extended to generalized functions of Colombeau type in 

[IIlEllllo]. 

In this Chapter we investigate how the generalized wave front set of a Colombeau function 
u transforms under the pullback by a c-bounded map /. Our main theorem is a result 
corresponding to the classical Theorem [23^ Theorem 8.2.4]. 



5.1 Transformation of wave front sets 

In order to obtain a microlocal inclusion (corresponding to [23^ eq. 8.2.4]) for generalized 
pullbacks, we need the notion of a generalized normal bundle and find a method to transform 
wave front set under c-bounded generalized maps. 

In this section we consider a c-bounded generalized map / € ^[f^i, ^^2], where C M", Q2 ^ 
M"* are open sets and assume that F C x 5*""^ is a closed set. 

Although the pullback of a Colombeau function u by any c-bounded generalized map / is well- 
defined, we cannot derive a general microlocal inclusion for the pullback without requiring 
further properties for the generalized map /. We define an open subdomain D-'^ of $7i x S""~^, 
where the generalized map (x, rj) 1— > -^(^/^(x)?/ has certain properties which are needed to obtain 
a microlocal inclusion relation. This leads to the notion of a generalized normal bundle. 

Definition 5.1.1. Let f G ^[01,572] be a c-bounded generalized map, then we define the open 
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5 Microlocal analysis of generalized puUbacks of Colombeau functions 
set D-^ by 

Bf ■.={{x,ri) gQiX S""-^ I 3 neighborhood X x V C x S°'"^ of (x,??) 
and a positive net of slow scale ((Te)^, Bet, /3 g]0, oo[, 3e' G]0, 1] : 
inf \a^^df^{x)ri\ > a and 

{x,r])eXxV 

sup \aE'^df^{x)rj\ < P for all e < e'}, 

{x,r])&XxV-^ 

where V-^ := {ij G S'^~^ \ Jrjo G V : {r],r]o) = 0}. Then the generalized normal bundle of f is 
defined by 

:= {{y,r]) en2X S"^-' \ {x,y) G Graph(/), (x, r?) D^}. 
This is in correspondence to the classical normal bundle. 

Furthermore, we define the wave front unfavorable support of f with respect to a closed set 
r C ^2 X S""-^ by 

\J'{T) := {x&W'\ {x, y) G Graph(/), (x, r?) D^, (y, r?) G T}. 

Example 5.1.2. Consider the c-bounded generalized map defined by 

fe{x,y) = {x + Jey,X -JsV), 

where {'ys)s is a positive net converging to zero. 
Then the transposed Jacobian is 

which is constant with respect to the {x,y) variable, so we put := ^dfe{x,y). 

We observe that M/^rj = {rji + r]2,je{vi ~ V2))- Setting vi := (l/\/2, — I/a/S) and V2 '■= 
(1/V2, 1/^2) we observe 

laeM^vil = -\/2cr£7£, \asMeiy2\ = V2ae 

and {MeVi,M^U2) = 0. 

We can write any 77 G -S^ in terms of the orthonormal basis {i^k)k=i,2 by r? = {rj, i^i)i'i+{ri, 1^2)^2- 
Note that (77, i/i)^ + (77, 1/2)^ = 1 since \ri\ = 1. 

For any 77 G we have 

\a,M,ri\'^ = \a,M,{{rj,ui)ui + (77,1/2)1^2)^ 

= (?7, z/i) V^M^z^ip + (77, 1/2) VeM^z^aP 

= 2ali{v, vxfll + ^2?) = 2a,2(l - iv, i^ifil - 7s))- 

Now for any 770 G with 770 7^ ±1/1 we can find a neighborhood Vq such that there exists a 

6 g]0, 1] with the property that |(?7, i^i)! < 1 — 6 holds for all 77 G Vq. Choosing '■= 1 we 
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5.1 Transformation of wave front sets 



obtain 

inf \ae^dMx)rj\ = ^fla, inf ^1 " " tI) 



and 



sup \ae^d^e{x)r]\ = \f2ae sup \J\ - (??,z^i)2(l -72) 
(z,r))eM"xyQJ- ??esi 

so it follows immediately that D-^ D M2 >^ (S'Vjiz^i})- 

It remains to check whether M2 x {±1/1} belongs to D-^: Let y C 5^ be a neighborhood of i^i, 
then it follows that 

inf |a/d/,(x)r/| = ^Pla, inf y^l - (r?, - 7I) = V2ct,7„ 

(x,»?)eiR"xy j?ev 

and 

sup \a^dU{x)ri\ = V2ae sup \/l - (r/,i/i)2(l - 7I) = \/2cre, 

(x,77)eR"xy-i- jjev^ 

since the set V'^ by definition contains a vector i/' with (i^i, u') = 0. Any slow scaled net ((Te)e 
with the property that \/^asje is bounded away from zero, satisfies ae ^ 00 (since 75 tends 
to zero). It immediately follows that niM/^ x {±2^1}) = , thus = M x (S'^\{±z/}). 

Remark 5.1.3. Note that any non-empty closed set T <^ Q2 x S"^~^ can be considered as a 
set- valued map 

supp(r) :=7ri(r)^5'"-\ y^Ty, 

where Fy := 7r2(r n {{y} x S"^~^)). Then Theorem 11.2.141 yields that y 1— > Ty is upper semi- 
continuous, i.e. for all y G supp (T) and W some open neighborhood of Ty, there exists some 
open neighborhood Y C Q2 of y such that Ty = rynsupp(r) ^ where Fy = IJy^zy^y 
as in Definition 11.2.51 Due to Proposition 11.2.121 the upper semi-continuity of T yields that 
supp (F) is a closed subset of 1^2) so we obtain for all y supp(r), that there exists some 
neighborhood y of y such that Fy = 0. 

Finally it is straight-forward to prove that for all Iq <s ^2 and W some open neighborhood 
of Tyq, there exists an open neighborhood y C of Yq such that Fy C W. 

Lemma 5.1.4. The generalized normal bundle N-^ and the wave front unfavorable support 
U^(r) of f (with respect toT) are closed sets. //N-^ n F = 0, then U^(F) = 0. 

Proof. Let xq U^(F) and put Yq := Graph(/)a;(, , then since is open there exist open 
neighborhoods Xq of xq and Vq of Fy^ such that Xq x Vq ^ D-^. As pointed out in Remark 
I5.1.3l there exists a neighborhood y of Iqi such that Vq is still a neighborhood of Fy. The upper 
semi-continuity of Graph(/) provides a neighborhood X' of xq such that Graph(/) n {X' x 
^2) Q X' xY. Now we can choose a smaller neighborhood Xi C X' n Xq (such that X' n Xq 
is a neighborhood of Xi) of xq, such that XixVq D-^ and Graph(/) n {Xi x $12) ^ Xi x y. 
For all xi £ Xi and y £ with {xi,y) £ Graph(/), it holds that (xi,y) G Xi x Y and if 
(y, r]) G F it follows that rj G Fy C V and thus (xi, 77) G D-^. It follows that xi U^(F) for all 
xi G Xi, Xi n U^(F) = and thus U^(F)'= is an open set. □ 
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5 Microlocal analysis of generalized puUbacks of Colombeau functions 



A substantial step in the proof of the main theorem is the application of a generalized sta- 
tionary phase theorem (cf. Appendix B). In order to obtain a lower bound for the gradient 
of the occurring phase function we have to consider the map 



S' 



n—1 



T 



\T 



dfe{x)ri\ 



Lemma 5.1.5. Let W C. S"^ he a neighborhood of some r]o G S*" , then we have for any 
'r)i,r]2 &W that 



G := {r? G S"^'^ \3r]' : {r]' , rf) = 0}. 



Proof. Assume W C C5(ryo) := {?? G S'^~^ \ {rjo,r]) > 1 — 5} for some 6 G]0, 1[. If no such 
S g]0, 1[ exists Co(77o) C W implies W-^ = 5""^, so the statement is trivial. Observe that 



\vi + ml 



and since 



which implies 



( 



m + r]2 m- m 



{{m,vo) + {m,m)) > 



2(1-5) 



V2(l + (r?i,r/2)) 



>l-d 



vi + ml' \m - ml Ivi + mllm - ml 



m-V2 
\m-m\ 



□ 



So we introduce the following notation: 
Lemma 5.1.6. We define the Colombeau map 

, ^ / N ^dfs(x)r] 

on the domain T)^ , then 77 1— > Ms{x,rj) defines an equi- continuous Colombeau function for 
fixed X. 

Proof. The map (7 : ^ ||| is equi-continuous on W^/Bs{0) for any fixed 6 > since 





1 


^ \v\ 


1 


lei 1^1 


~W\ 


~W\ 



{--J + l^jllH-ISD 



Let {xo,rio) G Df, then we can find some neighborhood X x V C Q.i x S"" ^ of {xo,rjo) such 
that 

inf las'^dfs{x)r]l > a and sup \ae'^dfeix)r]\ < /3 

{x,r,)eXxV {x,r])eXxV-^ 

for some slow scaled net {(Te)^ and a,/3 g]0, oo[. Then we conclude that 

2 



\M,ix,rj) - M,{x,Ol <^,. . ^ , 

Vdfe{x)r]l 



dfeix)r]-^dfe{x)C 



\ae^dfeix)r]\ 
|C-r?|<2|e-7?| 



ae dfe{x){ri - < ^ ■ f 1 T^f t \ I 

^T^i{x,r,)exxv We^dfe{x}ril 



a 
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5.2 Generalized puUbacks of Colombeau functions 



uniformly for all x G Xq and S,,rj £ V. □ 

Remark 5.1.7. By the preceding Lemma we have obtained that [{t] Mi;{x,7]) Iq/)^] is an 
equi-continuous Colombeau map, thus it follows by Proposition 12.1.51 that 

Graph(M)a;^^ = |^ ^ S*""^ | 3 a net {xs)e in ^^i : lirnxg = x € 

and y € CP ((MT-(£-)(xe, ry))^) for some map r £ T} , 

which simplifies the determination of the generalized graph Graph(M) considerably. 

Definition 5.1.8. Let f G ^[l^i, ^^2] be a c-bounded generalized map and assume that / 0. 
We define the puUback of some closed set P C Jl^ x 5"^^ by 

rr := Graph(M) o ({.} x (P o Graph(/)).), 

where {•} x (PoGraph(/)), denotes the set-valued map defined by x ! — > l-^} ^ (P ^ Crraph(y))a;, 
using the notation from Definition \1.2.5\ and \1.2.16l 

//D-^ = we put if*T)^ = 5"*-^ for all xG^i. 

Remark 5.1.9. Obviously the set-valued map /*P is defined as the composition of upper 
semi-continuous and locally bounded set- valued maps, thus /*P is again an upper semi- 
continuous and locally bounded set-valued map Qi 1— > J^q[S^~^). According to Theorem 
11.2.191 the pullback is defined by 

(rr), = y Graph(M),,^ 

where P^; := UyeGraph{/)^ 'The upper semi-continuity of /*P yields that for any neighbor- 
hood W C S"""-^ of (/*P)x there exists some neighborhood X of x such that (/*P)x C W. 



5.2 Generalized pullbacks of Colombeau functions 

In this section we prove the main result which gives a microfocal inclusion relation for the 
generalized pullback of a Colombeau function. The proof of the theorem relies on a generalized 
stationary phase theorem, the details of which are discussed in the Appendix B. 

Definition 5.2.1. Let f £ 01^11,^2] be a c-bounded generalized map. Then we call f slow 
scaled in all derivatives on the open set Xq C Qi, if for all a € N" there exists slow scaled 
nets {ra^e)e such that 

(5.2.1) sup |5"/e(x)| < Cara,s as £ ^ 

holds, where Ca are constants. Furthermore we call f slow scaled in all derivatives at xq £VLi, 
if there exists a neighborhood Xq of xq such that \5.2.1\) holds for all a G . 

Define the sets 

:= {x £ \ f IS slowscaled in all derivatives at x} 
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and 

Kf{u) := Pi ^i(Graph(/) n supp(u - k) x ^2). 

Remark 5.2.2. The set 7ri(Graph(/) n Oi x supp(u — k)) is closed for all k ^ 0.2 due to 
Proposition 11.2.12] since Graph(/) is upper semi-continuous and supp(u — k) is closed. This 
implies that K'^(ti) is closed in the relative topology of Vli. 

Lemma 5.2.3. Let f £ Q[0,i,il.2] be a c-bounded generalized map. Furthermore let xq ^ 
Uf(r) U (S-^)^ and Yo := 7r2(Graph(/) n {xq} x Og)- U W is some open neighborhood of 
{f*T)xQ , then there exist neighborhoods X of xq, V of Tyq and Y of Yq with the following 
properties: 

XxVC D^(r), {f*r)x C W, fe{X) C y for ah e < e\ 
(/e) is slowscaling in all derivatives on X and Fy ^ V 

Furthermore there exists a positive constant c > 0, such that 

inf \MJx, ??) - CI > c for aU e < e' 

(x,r;,^)eXxVxVK= 

holds (we are using the notation from Lemma [5.1.6\) . 

Proof. From xq U^(F) it follows that {xq} x Fy^ is a compact subset of the open set 
Thus we can find some neie hborhood X' xV Q of {xq} x Fy^. Since Xq E S-' we can 
assume without loss of generality that / is slow scale in all derivatives on the compact set X' . 

Since is a neighborhood of {f*T)xQ we have that 

(/*F),„ = 7r2(Graph(M) n {xq} x Fy^ x S^~^) d W 

where Graph(M) denotes the generalized graph of the generalized map defined by (x, rj) 1— > 
M£(x,r]) on the open domain D'^. By Lemma 12.1.121 there exist neighborhoods X",V" of xq 
resp. Tyo such that 

(5.2.2) Me{X" X V") C W 

for all £ < e' . Let V := V D V" , then by Remark 15. 1.31 we can find some neighborhood Y of 
Yq such that Fy C V. By Lemma |2 . 1 . 1 2 1 there exists some neighborhood X'" of xq such that 
fe{X"') C Y for small e. Let X := X' D X" n X'" and Z := X x V x W. Then there exists 
(x£,y£,Ce) G Z (note that (x,r/, C) 1-^ |M£(x,r/) — C| is a continuous function for each e e]0, 1] 
and Z is a compact set) such that 

Ce := ^ inf |M,(x, ??) - ^\ = |M,(x„ tj,) - Ce\ 

holds for some net (xg, Ce)^. By (|5.2.2p we have that Ce > c > holds for all e < e' , where 
c is a constant. □ 
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5.2 Generalized puUbacks of Colombeau functions 
Lemma 5.2.4. Let v,w £ R" , then the inequality 

I 1 2 1 1 2 

\av - (1 - a)w\'^ > ^'^'J^'^^ (l - {{v/\v\,w/\w\)f) 
holds for all a G [0, 1]. 

Proof. Observe that 

\av — (1 — a)w\'^ = Q^|up — 2q(1 — a){v,w) + (1 — a)^|wp 

= a^\v + w\'^ — 2a{\w\'^ + {v,w)) + \ w\'^ 

and differentiation in a sliows that the expression has a local extremum at oq := {{v,w) + 
+ Since the second derivative with respect to a is ji; — if p > for v 7^ — w (the 
case V = —w is trivial), we have a local minimum at ao- Finally we obtain 

I ^1 2 1 ^ 1 2 

lav — (1 — a)w\'^ > \aov — (1 — ao)w\'^ = (l — (v/\v\,w/\w\)'^) , 

\v + w]"^ ' 

by checking the boundary case a G {0, 1} we verify that it is a global lower bound. □ 

Theorem 5.2.5. Let f G t/[r2i,ri2] o c-bounded generalized map with representative {fs)e- 
For u G: with representative (ue)e, we define the pullback f*u by 

f*u:={u,if,ix)))s+Mini). 

It satisfies the microlocal inclusion relation 

WF(/*'u) C /* WF(n) J (u^(WF(u)) x 5"-^) \J ((K-^(n) n (S^)") x 5"-^) . 

Remark 5.2.6. It is apparent that the microlocal inclusion relation for WF(/*ti) is splitted 
into three different parts: 

• /* WF(n): This part corresponds to its classical counterpart (cf. [ 231 eq. 8.2.4.]). 

• U^(WF(n)) X S'"^-'^: According to Lemma 15.1.41 this set is non-empty whenever N-^ n 
WF(u) is non-empty. In the classical case the pullback of a distribution u by a smooth 
map / is not defined if N-^ n WF(u) is non-empty, so there is no classical counterpart 
to this set. 

• (K^(n) n (S^)^) X 5"-^: {Kf{u) n (8-^)'=) is the set of ah points x G ^i, where the 
generalized map / is not of slow-scale in all derivatives, whenever the Colombeau map 
u is not constant in a neighborhood of Graph(/)a;. This set is required since we allow 
non-regular Colombeau maps /, which are only c-bounded. So the singularities of the 
map / may cause singularities of f*u. There is no classical correspondence, since in |23t 
Theorem 8.2.4] the map / is assumed to be smooth. 
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5 Microlocal analysis of generalized puUbacks of Colonibeau functions 



Proof of Theorem 1 5. 2. 51 According to [18, Proposition 1.2.8] the pullback f*u is a well- 
defined Colombeau function in Let (xo,Co) S x and Yq = {y £ \ ixo,y) G 
Grapti(/)(Oi)}. We set T := WF{u) and Ty := {rj \ {y,rj) £ T}. In order to prove the 
statement we show that if (xo,^o) /*rUU^(r) x Q (K-^(u) n (S^)'=) x 5"-^ it follows 
that (xo,eo) 0WF(rn). 

So let us assume that (xc-^o) /*rUU^(r) x 5""^ U (K^ (u) n (S^)^) x 5""^ Then we have 
Co (/*r)xo and xo U^(r) IJ K^{u) D {S-^y. This allows splitting the proof in two parts for 
the cases xq K-^(n) U U^(r) and xq (S^)^ U U^(r). 



In the first case (xq K^{u) U U^(r)): Since xq K'^(u) there exists a generalized constant 
A; G with representative {k^)^, such that Graph(/) n {xq} x supp(u — k) = 0. Hence for all 
y £ Yq it holds that y ^ supp(u — k). Since Yq and supp(u — k) are two disjoint closed sets, 
we can find a closed neighborhood Y of Yq such that Y n supp(u) = 0. Let x G ^^^(^2) with 
the property that x = 1 on some compact neighborhood Y' C Y° of Yq and supp(x) C Y. By 
Lemma 12.1.121 there exists a neighborhood X' of xq and e' g]0, 1] such that fe{X') C Y' for 
£ < e' . 

In order to show that ^0 (WF(/*ti))a;(j for xq K''^(n) we have to find a smooth function 
99 with support on a neighborhood of xq, such that J^{f*uip) is rapidly decreasing on some 
neighborhood of ^o- We choose 99 to be a smooth function with supp((/?) C X'. For all e < e' 
the identity f*u ■ (p = f*{xu) ■ ip holds, since the functions x and ip where chosen such that 
(x o /e) • (/J = 1 for e < e' . We have that x ' (^e ~ ^e) G (^^2) and it follows that f*u ■ 99 is 



The second case (xq U^(r)U(S'^)'^ and (/*r)a;(,) is more difficult to prove: Again we have 
to find a smooth cutoff function if with compact support containing xq, such that J-{f*Uip) is 
rapidly decreasing on some neighborhood of ^o- The main step will be the application of the 
generalized stationary phase theorem (cf. Appendix B), which requires a suitable lower bound 
on the norm of the gradient of the occurring phase function (it will turn out that xq U^(r) 
is essential) and slow-scaledness of the phase function in all derivatives in a neighborhood of 

Xq. 

First of all we note that (/j) is slow scaling in all derivatives in some neighborhood of xq, 
smce Xq € SA 

Now let W C S"""^ be an open neighborhood of {f*T)x^^ such that is a neighborhood of 
(this is possible since {f*'^)xo)^ then Lemma [5.2.31 implies that there exist neighborhoods 
X of xo, Y of Yq, and V of Fy^ such that {f*V)x Q W, fe{X) C Y for all e < e', Fy C V, 
and 

(5.2.3) inf mix, Til) -^i\>d>0 

holds for all e < e' , where d is some positive constant. Without loss of generality (/e) is slow 
scaling in all derivatives on X. 

Let X £ C^{fi2) with the property that x = 1 on some compact neighborhood Y' C Y° of 
lo and supp(x) ^ Y. By Lemma [2.1.121 there exists a neighborhood X' of xq and e' g]0, 1] 
such that fe{X') C Y' for e < e' . Without loss of generality we can assume that X' C X. 
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5.2 Generalized puUbacks of Colombeau functions 



In order to show that WF(/*n)2,.Q we are going to prove that there exists a smooth 
function ip with support on a neighborhood of xq, such that J-{f*uip) is rapidly decreasing 
on (which is a neighborhood of ^o); which imphes (xo,^o) WF(/*n). 

We choose to be a smooth function with supp(99) C X'. For ah e < e' the identity 
{f^)ue (p = f*{xu£) (p holds, since the functions x ^^id ^ where chosen such that (x°/e)"V = 1 
for e < ei. 

Set y := {ry G 5™"! | 3??i G A G M+ : 7/ = A • r/i} and W := e S"^-^ | 3^ G VF, A G M+ : 
^ = A • ^i}. Obviously !F(ux) is rapidly decreasing on V^. We have that 

\T{{ru\m)\ = \H{rxu)em)\ 



XUe{ri) I exp(i(/£(x),r?) - i{x,i))Lp{x) dx d-q 



XUe{r])Ie{tv) dv 



XU~e{ri)Ie{(.,ri) dr] 



V 



+ 



XUe{r])Iei(.,v) dr] 



where we have set 



r]) := j exp {i{fs{x),r]) - i{x, dx. 
We intend to apply the stationary phase theorem IB. 31 (cf. the Appendix B) with 

\^\ + \v\ 



ifeix), 



7] 



{x, 



Thus we have to find a bound for the gradient of the phase function = |^d/e(x)ar/i — (1 — a)^i | 

from below for all a G [0,1], rji = ri/\rj\, = (,/\(,\ and x £ X. By optimization in the pa- 
rameter a (according to Lemma l5.2.4|) and using the notation of Lemma 15.1.61 we can bound 
the expression 



inf 

xex,riiev,iieW,ae[o,i] 



^dfe{x)arji - (1 - a)Ci| 



from below by 



inf 



''dfeix)vi\ 



(x,rii,^i)GXxVxW'' \'^dfe{x)r]i 



inf 



^dfe{x)r]i\ 



> 



ix,r,i,(,i)eXxVxW'= 2 \'^dfe{x)r]i 

1 \^dfe{x)m\ 



inf 



\Meix,rji) - ^1 



inf 



{x,r)i,(,i)eXxVxW'' 4 |^4fe(rc)?7i| + 1 {x,rii,^i)eXxVxW' 



\Meix,m)-Ci\ 



> Ca. 



-1 



inf 



iM,(x,?7i)-eii, 



{x,rji,(i)&XxVxW'' 

where a,, is the slow scaling net from the Definition 15.1.11 By ()5.2.3p we have 

inf \Me(x,'ni) -^i\> d>0 

{x,r)i,^i)eXxVxW'^ 
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5 Microlocal analysis of generalized puUbacks of Colonibeau functions 



for e < e' , where d is a positive constant. Thus the gradient of the phase function is uniformly 
bounded from below by > Cd ■ o^^ for all (x, ^, r/) G X x F x W^. The stationary 

phase theorem IB. 31 yields 

(5.2.4) \m.r^)\ < c^e-^i + iei + ir?!)-" 

for all q € Nq. Note that we use C,Cp,Cq and C^,/ as generic constants. In the case 
where (x, ?7, ^) £ X x x the stationary phase theorem (now with the phase function 
exp(— i(x,0)) gives 



\^\'m,ri)\ = \c\ 

(5.2.5) 



exp {-i{x, 0) exp {i{fe{x),r]))ip{x) dx 



J2 \D^{exp{im,{x),v/\v\))^)\ 

\a\<k 

and by repeated use of the chain rule we obtain the estimate 

m 

\a\<kl3<a ^^^0 1=1 

E ^^(71,.. .,7/) n ™P \d^HUx),v/\v\)\ 

7l+...+7,=/3 l<j</^^^0 

where 7i + 72 + • • • + 7; denotes a partition of the multi-index /3 in exactly I multi-indices, 
that add up componentwise to (3. Using the notation \g\k := max|Q,|=fc \d^f{x)\, we can bound 
the expression by 



< Ck,2 E EEl^l',^+™=^ I n _«_^P_^_ \{feix),V/\v\\,l 



\a\<kl3<al=0 "' \l<i<l x<^X,rjeV'= 

sup \{fe{x),v/\v\)\j) 



m ( 

<Ck2Y yy\v\^ max 



\a\<k(3<al=0 \xeX,rjeV- 

<C'fe,3(l-h|r?|)''max max sup \{fe{x),ri/\ri\)\A . 

1=1 0<j<k~l+l \,^^x,v&V- J 

Since fs{x) is of slow scale in all derivatives at xq, it holds that for all j £ No there exists 
constants Cj and slow scaled nets r^j such that 

sup \{fe{x),Vl)\i <Cjre,j 
(x,77i)eXxV''= 

holds for small e. Summing up the estimates of the form (I5.2.5P give that there exists constants 
Cp such that 

(5.2.6) |j^(^,^)|<Cp(l + |r?|)-^'e-i(l + |e|)^ 
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holds for all p €No and (r?, ^ ^ We observe that 



\Hiru)sm)\ 



<C,e-i / \^,{r,)\{l + \C\ + \r,\)-idr, 



+ + / |xS^(r?)|(l + |r/|f d?? 



holds and using ()5.2.4p and ()5.2.6p leads to the upper bound 

C7,e-isup (l + |r?r'^x^(r?) /"(I + |e| + |r,|)'^-'? d?? 

ri€V J 

V 



Finally we set k := q — p — m and I := p + n — 1 and obtain 

C^e-^sup {I + Mr'^xueir]) [{l + lCl + lvD'-'dr^ 

V 



<Cp,„e-'-"(i + iei)~^ 

for all ^ G and Cp^q^rn,n some constant depending on g, m and n. It follows that 

□ 



5.3 Examples 

Example 5.3.1 (Multiplication of Colombeau functions). This example was presented in 
[301 Example 4.2] in order to show that an inclusion relation like in [23^ Theorem 8.2.10] for 
the wave front set of a product of distributions, cannot be extended to Colombeau function 
with wave front sets in unfavorable position. 

Consider the Colombeau functions u and v defined by 

Ue := e~^p{£'^{x + -^ey)) 

where 7^ is some net with lim£_>o7e = 0. Note that in [30, Example 4.2] 7^ := e^^'^. These 
Colombeau functions are both associated to 5{x) l(y) and the wave front sets WF(m) = 
WF(f) = {0} X M X {(±1,0)} are in an unfavorable position. 
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5 Microlocal analysis of generalized puUbacks of Colonibeau functions 



We are going to apply Theorem ()5.2.5p in order to calculate WF(u ■ v). First we rewrite 
u ■ V = f*i.{5), where fe{x,y) = (x + 7ey,x — 7ey) with 7^ some net tending to zero. In 
Example ([5X2]) we aheady showed that = R"^ x 5V(±l,=Fl)- From [29l Theorem 15] it 
follows that WF(i{5)) = {(0,0)} x S^. The wave front unfavorable support of / with respect 
to WF{l{5)) according to Definition 15.1.11 is 



Theorem 15.2.51 gives that 

WF(u-f) C (0,0) X S^, 

which is consistent with the result in [30\ Example 4.2]. 

Example 5.3.2 (Hurd-Sattinger). Let us consider the initial value problem 

dtu + ed^u + e'u = o 
5.3.7 * , 

u{0,x) =uo 

where G G : is defined by Ge(x) = p^^ * H{—-) with p^^ = ^p{-/je) where p G 5(M) 
and / p{x) dx = 1 and 7^ = l/log(l/e) is a net of slow scale. For the initial value we 
choose uq := i{5-sg) a delta like singularity at —sq (for a positive sq > 0). The Hurd- 
Sattinger example was first given in pi] (it was shown that it is not solvable in L^^^, when 
distributional products are employed). It was further investigated in [29j with methods from 
Colombeau theory. In [T7j the wave front set WF-y (with respect to the slow scale net 7) 
of the Colombeau solution was calculated. For sake of simplicity we do only consider the 
standard generalized wave front set WF, which is smaller since it neglects the singularities 
coming from the coefficient 0. 

We can write 

®e{x) = g{x/je), 

where g{x) := p{z) dz with lim^-^^-oo g{x) = and lima;_++oo g{x) = 1 (this implies that 
Ran((jr) = [a+,a_] for some a_ < 0,a+ > 1). 

We have already considered (cf. Example I3.3.8P the ordinary differential equation for the 
characteristic curves 




dsUs) 



X 



and we obtained 




X < 

X = 0,0 < i 
X > 0. 



If we set fe{t,x) := (t, ^^(0; t, x)) and / : 



[{fe)e] then 



u = f*{l0 no) • dxf{t,x) 
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is a solution of the initial value problem (15.3.70 . 

We observe that the derivatives of the characteristic flow satisfies 



e,(ge(0;t,x)) 



which follows from Remark 13.1.91 So dxfeit,x) is slow scaled in all derivatives, due to the 
fact that @e{x) = g{x/^e) where 7^ is a net of slow scale. We may conclude that WF(m) C 
WF(/*(1 CStio)), which enable us the apply Theorem 15 . 2 . 5 1 in order to estimate the wave front 
set of u. 



We put M = [iMe)e] G a[D^, 5^] defined by 

Ms{t,X,7]) -- 

as in Lemma l5. 1.61 Note that 

^dfe{t,x) = 



"dfe{t,x)r] 
''dfe{t,x)r]\ 



1 -e,(ee(0;t,x)) 



For r]± := (0, ±1) G (which are the only irregular directions coming from the wave front 
set of 1 (gi no), we obtain 



M,it,x,r]±] 



\dfe{t,x)^7J±\ 



VVi + e,(x)-2Vi + 0.(^)V 

Remarkably the result does only depend on the coefficient Ge(x). In Example 12.1.41 we 
obtained 

Graph((ee)e) = ({x G M | x < 0} x {1}) U (^{0} x Ran(c/)) U ({x G M | x > 0} x {0}) , 



so Theorem 12.2.11 yields 
Graph(M(-,r/±))t,^ C 



±1 



Graph(e) 



{(0,±1)} 

Uae[a_,a_|.]{( "F 



t,x 



for X < 
for X > 
}, for X = 0. 



Furthermore we observe that 
(WF(1 uo) o Graph(/))t,, 

WF(1 Mo)t,min (x-t.O) X <0 

^ < U.e[-t,o] WF(1 ^ no)t,. x = 
WF(1 uo)t,x = ^ > 



{(0,±1)} x<0,i = x-so 
{(0,±1)} x = 0,t>-so 
X > 
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5 Microlocal analysis of generalized puUbacks of Colonibeau functions 



and according to Definition 15.1.81 we obtain 



(r (WF(1 uo)))t,. = Graph(M) o ({(t, x)} x (WF(1 ® uo) o Graph(/))t,,) 

r {(^1,±1)} t < so,x = t- so 

= { U.e[._,.,]{(7&'7&^)} * = -o,:. = o 
I {(0,±1)} t>so,x = 

and Theorem [5X5] yields WF(u) C /*(WF(1 O uq)). 
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A Sets and measure theory 



In the Appendix A we have collected a few auxiliary measure- and set-theoretic results. Most 
of the presented material can be considered common knowledge, nevertheless for sake of 
completeness we decided to include a proof for all results, adapted to the special needs and 
notation of earlier sections. 

Definition A.l. Let I ■.=]0, lq] with lo £]0, 1] and {Mi^)^^j a family of subset o/M". 
(Mjtgx ii o-n ascending family of sets, when 2 for a < (3. 
(Mjtgx ii o-n descending family of sets, when C Mg for a < p. 
Furthermore we define the limes superior by 



limsupMt -=0 U 

and the limes inferior by 



liminf M, := (J Q M,. 

ti6Xte]o,ti] 

In general it holds that liminft_»o -^t Q liinsup^_^Q M^. When we have M := liminf(,_»o = 
limsup^^o -^t? (-^Jtex converges to the set M and use the notation 

lim,^oM, :=M. 

Lemma A. 2. If {M^)^^x is an ascending (resp. descending) family of sets, it converges 

\imM, = \jM„ 

resp. 

limM, = f|M,. 
tex 

Proof. If (Mjtex is an ascending family of sets, we have that IJte]o n] ~ Ui,e/ -^t 
-^ti ^ fltelOti] -^t' ^^^^ lim sup^_>Q M(, = Utex-^t — hm inf (,_»o -^t which implies the state- 
ment lim^^oMi, = [j^^jM^. The proof for the descending family is analogous. □ 

Lemma A. 3. Let (Mj^gx be an ascending family of sets such that M := lim^^o = 
Utex exists. If f : — ^ R is some function, then it follows that 



and 



sup f{x) = lim ( sup f{x)) = sup ( sup f(x)) 



inf fix) = lim ( inf fix)) = inf ( inf fix)). 
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Proof. Note that by Lemma [A2] it holds that M = \J,^jM,. 

We set := sup^jg^^ f{x) and since (Mj^gj is an ascending family we have that it is a 
monotone increasing net. As M = Utex have that < sup^g^/ /(x). It follows 

immediately that a := lim^^o ctt = supja^ | l G I}. For all x G M there exists some l £ I 
such that X G M^, thus f{x) < a,, which implies that supja^ | t G X} is an upper bound of 
the set f{M). As sup^.^^/ f{x) is the smallest upper bound of / we have a = sup^gj^,^ fi^)- 

We set /3t := inix^M, f{x) and since (Mj^gj is an ascending family we have that it is a 
monotone decreasing net. As M = Uigx-^t have that /?t > mix&M f{x). It follows 
immediately that (3 := lim^-^o A = inftgj (3^. For all x G M there exists some i G T such that 
X G Mt, thus /(x) > /3t which implies that inf^gj /^^ is an lower bound of the set f{M). As 
infjjgjv/ f{x) is the greatest lower bound of f{M) we have P = inf M. □ 

Lemma A. 4. Let {Mi_)i_^x be an descending family of sets in such that M := lim^^o-^t = 
ritex exists. If f : M" — > M is some function, then it follows that 

sup /(x) = lim ( sup /(x)) = inf ( sup /(x)). 

and 

inf f(x) = lim ( inf f(x)) = sup ( inf f(x)). 
Proo/. Note that by Lemma E2] it holds that M = f]^^jM,. 

We set := supM^ and since {M^)^^x is a descending family we have that it is a monotone 
decreasing net. As M = Cli^^jM^ we have that > supM. It follows immediately that 
a := lim^^o Q^t = infjot | l G I}. For all x G M we have that x G for all l € I, thus 
X < which implies that inf {a^ | z. G X} is an upper bound of M. As supM is the smallest 
upper bound of M we have a = supM. 

We set (3i := inf Af^ and since {M^)^^x is a descending family we have that it is a monotone 
increasing net. As M = Htex have that (3^ < inf M. It follows immediately that 

j3 := limt^o A = sup{/3t | i G X}. For all x G M we have that x G for all i G I, thus 
X > I3l which implies that sup{/3t | t G X} is an lower bound of M. As inf M is the greatest 
lower bound of M we have [3 = inf M. □ 

Lemma A. 5. Let f : 0, ^ M be a function, then for any bounded set C C there exists a 
sequence {xk)keN in C such that 

sup/(y) = limsup/(xfe). 

y&C k^oo 

Proof. Denote M := supy^Q f{y). 

If M = oo, we can pick find Xfc G C such that f{xk) > for all A; G N. It is now trivial that 
limsupfc^oo/(a^fc) = 

In the case M < oo: Let x^ G C, such that /(x^) < M, then we can find some x^+i such 
that M - /(xfc+i) < 1/2(M - /(xfc)). 

The last step followed by contradiction: Assuming M — /(x) > l/2(Af — f{xk)) for all 
X G C holds, implies /(x) < 1/2(M + f{xk)). The right-hand side is obviously an upper 
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bound, smaller than M, contradicting the fact that M is the smallest upper bound of the set 
{fix) I X G C}. 

Starting with any xi G C, we iteratively obtain a sequence with M - f{xk) < 2-^, thus 

sup f{y) - limsup f{xk) = 0. 

□ 

Lemma A. 6. Let f : Vti ^ 0.2 he a continuous map. If {M^)i_^x be an descending family of 
sets, such that linii^^Q = Cli^^jM^ is a compact subset ofi}i, then it follows that 

Iim/(MJ = /(limM,). 
t— >o t— >o 

Proof. First of all we note that /(H^gjMj C f]^^j f(M^). We proove by contradiction: 
Assume there exists some y E Clit^i fi^i.) such that y f{M). As f{M) is a compact subset 
of ^2 we have that there exists some e > 0, with i?e(/(M)) P|{y} = 0. The continuity of / 
implies that there exists some (5 > such that f{Bs{M)) C B^{f{M)). Now there exist some 
lq el such that M^g C B^{M). Assuming the contrary, that Bs{M) C for all t G T, would 
lead to the contradiction Bs{M) C lim,_^oM,. Finally we obtain f]^^jf{M,) C f{M,g) C 
fiBsiM)) 

y B,{f{M)) D f{Bs{M)) D f{Mj D f] f{M,), 

which contradicts the inital choice of y, prooving the statement. □ 

Lemma A. 7. Let h : M" M. be a upper semi- continuous function and K (s M", then the 
supremum sup^g^ h{z) is attained, i.e. there exists a zq G K such that sup^g^ h{z) = h{zo) < 
oo. 

Proof. Let (xfc)fcgN the sequence from the proof of Lemma lA. 51 By the upper semi-continuity 
we get that 

sup/(y) = limsup/i(xfe) < /i(limsupxfc) 

y€K fc— >oo k—*oo 

and since {xk)keN Q K we have that x := lim sup^^^^ a;;;, is a point in K. Furthermore the 
semi-continuity of h implies pointwise boundedness, thus h{x) < oo. □ 

Lemma A. 8. Let g : J x ^ M 6e a function such that 

(i) t ^ g{t, x) is Lebesgue measurable for all x € M, 
(a) X I— > g{t,x) is upper semi- continuous for almost all t £ J, and 

(Hi) there exists a positive function f3 G L\^^{J) such that < I3{t) for almost 

all t £ J, 

theng£Ll^{J;L^{W)). 
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Proof. It has to be shown that t i— > z)! is a Lebesgue measurable function. Let 

{Ki)i^fq be an ascending family of compact sets such that IJ/^i ~ ■'^"j then we have by 
Lemma IA.3I that 

/ \ 

sup h{t,z) = liminf sup \h{t, z)\ 

26M" \z&Ki J 

Since g satisfies property (ii), i.e. the upper semi-continuity of x ^ h{t,x) for almost all 
t G J, we conclude by Lemma fA. 71 that there exists a sequence {zi)i^fq with zi G Ki such that 

sup h{t, z) = liminf \h{t, zi)\ 

for almost all t £ J. Property (i) yields that {\h{t, zi)\)ii^n is a family of Lebesgue measurable 
functions. It follows by |21[ Theorem 11.12] that liminf;^oo \h{t,zi)\ is Lebesgue measurable. 
Since \h{t, zi)\ < P{t) for almost all t G J as assumed by (iii), we obtain due to [21^ Theorem 
12.24] 

/sup h{T, z) dr = / liminf sup |/i(r, z)! dr < liminf / sup |/i(r, z)| dr < / /3(r) dr 
M M ^ ' M ^ ' M 

for all M (s which implies g £ L^^iJ; L°°(IR")). □ 

Lemma A. 9. Let h : J x M" -^M be a function satisfying properties (i) - (ii) of Lemma \A.8l 
then for any Lebesgue measurable map ^ : J — > R"',t i— > ^(t), the composition 

t^h{t,^{t)) 

is a Lebesgue measurable function from J ^ M. 



Proof. First we proove the simpler case where x i— > h{t, x) is continuous for almost t £ J 
instead of property (ii). We approximate ^ by simple functions Snit) := l^fc=i an,/clyl„ ^(0 
(see [21 Theorem 11.35]), such that is the pointwise limit of {sn)nm- Note An^k are Lebesgue 
measurable subsets of J. We have 

n 

/(*,5n(t)) = J^F(t,a„,fc)lA„.,(t), 
fc=l 

SO t 1-^ f{t,Sn{t)) is Lebesgue measurable as the finite sum of measurable functions. From 
the continuity of x i— > f{t,x) it follows that lim„^oo /(^) ^^(t)) = f{t,^{t)) holds pointwise 
for almost all t £ J, thus t i-^ /(t, ^(t)) is Lebesgue measurable as the limit of measurable 
functions (see [211 Theorem 11.12]). 

The case where x i— > f{t,x) is only upper semi-continuous is much more difficult. The idea 
is to approximate this function by functions f{t,x) = inf^gfj (/^^(t, x), where {ipk)k&'N satisfy 
(i) and x i— > ipk{t,x) is continuous. First of all we set hk{t,x) := supyg^^^^^^.) f{t,y), which is 
again upper semi-continuous and fulfills hk{t,y) > f{t,x) for all y £ Bi/f.{x). By the upper 
semi-continuity of x i-^ f{t, x) we have that for all / G there exists a constant 6i of x such 
that f{t,y) < f{t,x) + 1/1 for all y £ Bsi{x). Choose /cq S N such that Biik^^{x) C Bsi{X), 
then /ifc(i, x) < /(t, x) + l// for all k > k^. It immediately follows that inffcgpj hk{t^ x) = f{t, x). 
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Let p e C'^{W) be positive with supp(/9) C ^^(O) and / p(x) dx = 1. Then we define 
iPk{t,x) := /^^(Q) h2k{t,x-{2k)-^y)p{y)dy by convolution, so ipk € Lj^^(J; C°°(M")). Observe 
that \x- {x- {2k)-^y)\ < {2k)-^\y\, thus x - {2k)-^y £ Bi/2kix) and h2k{t,x- {2k)-^y) > 
f{t,x) which imphes 

/ /(MW.)d. = /(M) 

Bi(0) Bi(0) 
for all {t, rr) G J X and /c G N. 

Setting Zk = {2k)~^y and w G Bi/2kix — Zk) we obtain \w — {x — Zk)\ < l/{2k), so — 
■fi"! — \zk\\ ^ l/(2fc) which implies |a; — < 1/k for all w G Bi/2kix — Zk). Conclude that 
h2k{t,x- Zk) := sup^gsi/(2fe){x-zfe) /(i,^^) < sup^gBi/fc(x) = hk{t,x). It follows that 

f{t, x) = inf x) > inf 99A;(t, > f{t, x), 
fceN fceN 

thus f{t,x) = inffcgN x). It follows by the first part of the proof that t i-^ ^k{ti(,{t)) is 
Lebesgue measurable which yields that t fit,(,it)) is Lebegue measurable as the infimum 
of Lebesgue measurable functions (see [21, Theorem 11.12]). □ 
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B A generalized Stationary Phase theorem 



Appendix B contains a straight-forward generalization of the stationary phase theorem for 
nets of smooth functions. We essentially follow the proof of the classical stationary phase 
theorem as presented in [23\ Theorem 7.7.5], but with a pedantic book-keeping on the e- 
dependence of the estimates. Due to the technical nature of this result, we decided to put it 
in the appendix. 

Definition B.l. Suppose that g £: and let {ge)e some representative, then we introduce 
the notation 

\ge\k ■■= XI l^^s'^l 

\a\=k 

and 

IJ'K,k,e{9) ■= SUp|5e|fc 
K 

t^KAeia) ■= max^i^,,,e(5r). 

Lemma B.2. Suppose that g 0, with an non-negative representative {ge)ej such that there 
exists e' > 0, with ge{x) > for e < e' and all x E il. Then it follows that for any d 
there exists some constant C such that 

X \d^ge(.x)\ < Cy^)^fi4^) 
\a\=l 

holds for all X ^ K and M is some compact set with K d M° . 

Proof. Let K be some compact set and M d Q such that K M° . Then we can find some 
(5 > such that Bs{K) C M holds. Choose some xq G K. By Taylor's formula we obtain for 
any x = xq + ae^ with k = 1, n and a G [—(5, 5\ (note that xq + ae^ G B^{xq) C M) that 

< ge{x) < ge{xo) + dkge{xo)a + \ae{k)a^ 

holds with 

ms{k) := max (sup \dlge{x)\,2ge{xo)6~^). 
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B A generalized Stationary Phase theorem 



Now we distinguish two cases: If 



meik) < 4(/£(xo), then we put a = ±5 and get 



5\dkge{xo)\ < ge{xo) + -^rne{k) 




< a/ {geixof + 8'^ge{xQ)me{v) + 5^me{k)ge{xo) 



= V ge{xo)Vge{xo) + 25'^me{k). 



In the case where S'^rrieik) > 4gs{xo) we set a = ±^/2gs{xo)/ms{k) < S to obtain 



Idkgeixo)] < \/me{k)/2^yg^{xo) + 




= ^/2me{k)^/ge{xo) 
We can finahy conclude that the estimate 

\dkge{xo)\'' < 5-^ge{xo)igsixo) + 26^m,{k)) 

< 6'^geixo){sup geix) + 26^me{k)) < max {6'^ ,2, 26~'^)ge{xo)nl[ 2 eig) 
xeK 

holds for all xq e K where S and A*m2£(5) independent of xq. The statement fol- 
lows by adding up the estimates Yl\a\=i\^'^ 3^(^)1 ~ Ylk=i\^kg£{x)\ and putting C = n ■ 
Vmax ((5-2,2,2(5-4). □ 

Theorem B.3 (Stationary phase theorem). Let u G Qc{^) with support K ^ Q and (f)^ € 
£{^) with the property that there exists an Eq > and m G N such that 



(1) 



inf (|(^;(x)|) > X, for aU e < Eq 



holds, with Ae some net tending to zero. Then we have that 





e 



is a Colomheau function in the u variable and it is bounded by 



(3) 




\a\<k 



where 



Lk,e := Cfcmax{l,//^^fc^^((/)£) ' } 



and M is some compact set with K C M° . 
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Proof. It is obvious that ([2]) is a well-defined Colombeau function. The proof follows closely 
the proof of classical stationary phase theorem in [23^ Theorem 7.7.1]. By Ni;{x) := \(p'^{x)\'^ 
we denote the square of the norm of the gradient of the phase function. Let 



and since 



iio-;—^ exp (iiocj)^) = dy exp {iuj4>£) 



we obtain after an integration by parts 

(using the notation introduced in Definition IB. ip . We prove by induction: For A; = we have 
the obvious bound 



/lie exp (iu;0e)(ix| < C sup 
K 



Assume that the bound ([3]) holds for power A: — 1, then we have that 
uj^\ J Us exp {iu;(j)£)dx\ ='^'^""^1 J {duUu,e)^^P {iuj<j)£)da 



k-l / n N 



m=0 ^ \v=l 



holds. In the next step we are going to show that 

m 

(4) m I < Mm,e I 



r — m 
2 

e 



V r=0 

holds. Again we are prooving by induction: For m = we have 



(5) ivi^|u,,,| = J^iV-^ 



2 \-:r—\ < n\uF 



dxy 



and let us now assume that ^ holds up to m — 1. Let a be any multi-index with \a\ = m. 
We apply 9" on 

NeUu^s = Us- 



dx,, 

and obtain 



f3<a 
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B A generalized Stationary Phase theorem 



It follows that 



(6) |iV.9"<| 



f3^0,l3<a ^' ^ /3<a 
- ^ ( ^We\l\Uu,e\m-l +^\(t>e\m-l+l\Ue\l ] < C {\Ne\i\Ui,^s\m-l + \(pe\l\Ue 



vi=l 1=0 
m-2 



m-2 



+ max {sup \Ns\m-l, sup \(j)e\rn-l+l) ^{\Uu,e\l + \Us\l) + 



U£\m-l\(Pe\2) ■ 



Now we can apply Lemma lB.2l to bound |A^e|i. This yields 



H=i 

\^e\l < C20V^ 



and we can verify that 



a-P+ei 



i=l 



i=l P<a 



< 7^2 E \4>e{x)\\a\-k+l\4>e{x)\k+l 
k<\a\ 



holds. This yields 



SUp|iVe|/ < K^y^ Snp\(t)e{x)\l-k+lS\ip\4)e{x)\k+l 



k<l 
2 



We introduce 



max (sup I iVe I / , sup I I i+i ) 



< max {Ki^l'K,m+lM^)'^^^^K,m+lM<^)) 

in order to simplify the notation. Using these bounds we can estimate ([6]) by 



t'm-2 



, 1=1 



m-2 



m—1 



z=i 



1=1 
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and the induction hypothesis ([!]) for m — 1 gives 

/ m— 1 r— m+i m—2 I 



2 

e 



r=0 1=1 r=0 

m , , m—2 l+l 

r — m + l 



V r=0 /=0 ^=0 

/m-1 \ i+1 , , 



2 



< CscT^^ E l^-l-^- ' ^ ^ecT^) maxmax{l, M;,J Y 

\ 1=0 / r=0 r=0 

It follows immediately that 

m 

(7) Y. W'^,e\m < Mm,s Y. 1""!^^ 



r — m 
2 



f r=0 

holds, where Mm,£ is defined by 

M„,, := 4'^) max max {1, Mi, J = ^^"'^n^^^ max{cj«, 1} 
is a generalized number. We are finally able to estimate (jB]) from above by 

k-l / m+l 



r — m — l m + l ^ 

e 



Lk-l,e Y, Yj \'^e\rNs ^ iVf 

m=0 ^ \ r=0 

< Lk-i,sMk,s{k - 1) X^sup (in^l^A^I" 

r=0 ^ 

The generalized constant L^^g is recursivly defined by 

k 

Lk,e := Lfc_i,,Mfc,,(fc-l) = (fc-l)! J]Mz 

1=1 

k l-l 

= C{k-l)lllaPl{m^x{aP,l} 



1=1 i=l 

(k \ k l~l 

J]max{l,^l^,^,(^,)2} J] J]max{l,^^,(0,)2} 
1=1 / 1=1 i=l 



Note that for < r < A; we have that 



Ni < min{l,inf < 



holds. 
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Abstract 



As the title "Generalized regularity and solution concepts for differential equations" suggests, 
the main topic of my thesis is the investigation of generalized solution concepts for differential 
equations, in particular first order hyperbolic partial differential equations with real-valued, 
non-smooth coefficients and their characteristic system of ordinary differential equations. 

In Colombeau theory there have been developed existence results that yield solutions for ordi- 
nary and partial differential equations beyond the scope of classical approaches. Nevertheless 
this comes at the price of sacrificing regularity (in general a Colombeau solution may even 
lack a distributional shadow). It is prevailing in the Colombeau setting that the question of 
mere existence of solutions is much easier to answer than to determine their regularity prop- 
erties (i.e. if a distributional shadow exists and how regular it is). In order order to address 
these regularity question and encouraged by the fact that the solution of a (homogeneous) 
first order partial differential equation can be written as a pullback of the initial condition by 
the characteristic backward flow, a main topic of my thesis deals with the microlocal analysis 
of puUbacks of c-bounded Colombeau generalized functions. 

Another topic is the comparsion of Colombeau techniques for solving ordinary and partial 
differential equations to other generalized solution concepts, which has led to a joint article 
with Giinther Hormann. A useful tool for this purpose is the concept of a generalized graph, 
which has been developed in the thesis. 

Wie der Titel "Verallgemeinerte Regularitat und Losungskonzepte fiir Differentialgleichun- 
gen" andeutet, ist das Hauptthema die Untersuchung von verallgemeinerten Losungskonzcptcn 
fiir Differentialgleichungen, insbesondere fiir hyperbolische, partielle Differentialgleichungen 
erster Ordnung mit reellen, nicht-glatten KoefHzienten und deren charakteristisches System 
von gewohnlichen Differentialgleichungen. 

In der Colombeau Theorie gibt es Resultate die Losungen fiir gewohnliche und partielle Differ- 
entialgleichungen liefern, die ausserhalb der Reichweite gangiger klassischer Losungskonzepte 
liegen. Nichts destotrotz werden dicsc Resultate auf Kosten der Regularitat solcher Losungen 
erzielt (im AUgemeinen folgt nicht, dass Colombeau Losungen einen distributionellen Schat- 
ten besitzen miissen) . Es ist meist einfacher im Rahmen der Colombeau Theorie die Existenz 
einer Losung zu zeigen, als der deren Regularitatseigenschaften zu bestimmen (d.h. ob ein 
distributioncllcr Schatten existiert und wie regular dieser ist). Um solche Regularitatsfragen 
zu behandeln und durch die Tatsache bestarkt, dass die Losung einer (homogenen) par- 
tiellen Differentialgleichung erster Ordnung als Pullback der Anfangsbedingung durch den 
charakteristischen Riickwartsfluss bestimmt werden kann, ist ein Teil meiner Dissertation mit 
der mikrolokalen Analyse von Fullbacks kompakt beschrankter Colombeau verallgemeinerter 
Funktionen befasst. 

Ein weiteres Thema ist der Vergleich der Colombeau Techniken zum Loscn gcwohnlicher 
und partieller Differentialgleichungen mit anderen verallgemeinerten Losungskonzepten. Ein 
niitzliches Werkzeug fiir diese Zwecke ist das Konzept eines verallgemeinerten Graphen, das 
im Rahmen dieser Arbeit ebenfalls entwickelt worden ist. 
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